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Abstract 
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The  powering  requirement  of  a  ship  is  one  of  the  most  important  aspeets  of  naval 
arehitecture.  Traditionally,  ships  have  been  tested  for  hull  resistance  using  hydrodynamic  tank 
testing.  Tank  testing  evaluates  models  of  ships  by  measuring  their  resistance  in  a  tow  tank.  This 
has  proved  to  be  useful  and  accurate,  but  is  very  time  consuming,  expensive,  and  has  inherent 
scaling  errors.  Because  of  these  reasons,  today  many  vessels  are  sold  on  the  market  without  any 
model  testing.  Another  set  of  design  tools  is  parametric  predictions.  Parametric  predictions 
contain  acquired  data  for  a  specific  family  of  hull  forms  and  use  key  hull  parameters  to  evaluate 
a  particular  design.  Parametric  predictions  run  quickly,  but  cannot  be  used  to  evaluate  any  hull 
outside  of  the  limits  of  previously  tested  hulls.  Computational  Fluid  Dynamics  (CFD)  codes  are 
a  much  newer  method  of  determining  the  powering  requirement  for  ships.  CFD  uses  numerical 
modeling  to  simulate  the  forces  which  act  on  the  ship.  Unlimited  testing  takes  a  fraction  of  the 
time  and  effort  as  compared  to  tank  testing.  However,  it  is  not  yet  wholly  proven  in  its  accuracy. 

An  investigation  has  been  made  using  tank  testing  and  CFD  to  predict  the  new  David 
Pedrick  Mk  II  Navy  44  Sail  Training  Craft's  (STC)  performance.  The  performance  prediction 
process  proved  most  accurate  when  using  tank  data  to  aid  the  CFD  calculations.  Using  this 
combined  method  along  with  an  aerodynamic  calculation  for  the  Mk  II  Navy  44  STC,  a  Velocity 
Prediction  Program  (VPP)  was  constructed.  Finally,  the  new  VPP  as  well  as  previous  parametric 
VPPs  were  used  to  predict  a  full  set  of  performance  figures.  These  tools  were  also  used  to 
improve  the  rudder  design  of  the  Mk  II  Navy  44  STC. 
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Performance  prediction  has  always  been  an  important  parameter  in  ship  design.  From 
cruising  boats  to  America's  Cup  yachts  to  naval  surface  ships,  the  speed  of  the  vessel  will  greatly 
influence  its  mission  effectiveness.  Today  there  are  three  methods  of  making  performance 
predictions  of  vessels: 

•  Parametric  predictions 

•  Hydrodynamic  tank  testing 

•  Computational  fluid  dynamics 

What  this  study  entailed  was  an  examination  of  these  methods  of  performance  prediction.  In 
order  to  research  these  tools  extensively,  a  sailing  craft  was  evaluated  because  of  the  complexity 
in  its  hydrodynamic  characteristics.  For  example,  a  naval  vessel  or  a  cargo  ship  is  typically  only 
evaluated  in  its  forward  motion.  A  sailing  vessel,  however,  must  be  evaluated  not  only  in 
forward  motion,  but  also  for  varying  amounts  of  heel  angle,  yaw  angle,  rudder  angle,  and  trim 
angle.  Furthermore,  the  lifting  surfaces  on  a  sailing  craft  such  as  the  keel  and  rudder  have 
tremendous  effect  on  the  performance  of  the  vessel  and  so  must  also  be  evaluated. 

For  sailing  craft,  parametric  predictions  often  take  the  form  of  Velocity  Prediction  Programs 
(VPPs).  These  parametric  VPPs  are  the  result  of  the  testing  of  many  hull  forms  in  the  tow  tank. 
VPPs  are  computer  programs  which  predict  the  hydrodynamic  as  well  as  aerodynamic 
configurations  of  sailing  craft.  Parametric  VPPs  take  key  measurements  of  the  hull  and  sailing 
rig  and  then  interpolate  through  supplied  data  to  make  predictions  for  any  vessel.  However, 
parametric  VPPs  are  restricted  in  their  ability  to  extrapolate  beyond  the  limits  of  their 
hydrodynamic  data.  For  each  wind  condition,  the  program  will  iterate  boat  speed,  heel  angle, 
and  sail  conditions  so  that  the  sum  of  the  aerodynamic  and  hydrodynamic  forces  and  moments 
on  the  vessel  is  zero.  The  speed  of  VPPs  is  what  makes  these  tools  so  valuable.  Today,  VPPs 
can  be  run  even  through  personal  computers  in  a  matter  of  minutes.  Parametric  VPPs  are  used 
mostly  for  preliminary  yacht  designs.  They  will  give  as  output  the  estimated  speed  of  the  boat 
for  any  wind  condition  as  well  as  the  angle  of  heel,  side  force,  and  drag. 

Over  many  years,  naval  architecture  has  refined  the  method  of  hydrodynamic  tank  testing. 
Simply  defined,  tank  testing  consists  of  measuring  the  forces  on  a  geometrically  scaled  model 
which  is  being  towed  through  a  tank  of  water.  The  measured  forces  from  the  tank  tests  can  then 


9 

be  extrapolated  to  predict  the  forces  on  the  full-sized  vessel.  For  sailing  craft,  an  elaborate 
testing  program  is  necessary.  Due  to  aerodynamic  forces,  the  towing  rig  must  be  equipped  to 
change  the  vessel's  yaw,  heel,  trim,  and  rudder  angles.  Although  drag  is  the  primary  output  from 
the  tests,  lift  and  yawing  moment  are  also  important  to  a  sailboat.  Besides  the  complexity  of  the 
testing  rig,  tank  testing  has  further  inherent  difficulties  in  scale  effects.  Though  the  model  can  be 
geometrically  scaled  down,  the  water  molecules  cannot,  leading  to  inherent  flow  differences 
between  the  ship  and  the  model.  For  this  reason,  there  are  many  numerical  adjustments  and 
engineering  assumptions  which  add  to  the  extrapolation  of  the  data. 

Computer  codes  known  as  Computational  Fluid  Dynamics  (CFD)  use  fundamental  fluid 
dynamics  equations  to  predict  the  resistance  of  the  ship  moving  through  water.  Using  a  digitized 
hull,  a  computer  puts  the  ship  into  a  virtual  tow  tank  and  calculates  the  lift,  drag,  yawing 
moment,  and  trim  on  the  body.  However,  the  accuracy  of  the  simulation  is  primarily  a  function 
of  the  quality  of  the  code  and  the  quality  of  the  digitized  hull.  Furthermore,  there  are 
complexities  due  to  viscous  effects  and  the  free-surface. 
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The  Performance  Prediction  Process 

The  intention  of  this  study  was  to  compare  these  different  performance  prediction  tools  and 
create  a  process  which  will  accurately  evaluate  the  Mk  II  Navy  44  Sail  Training  Craft  (STC). 
This  study  compared  the  parametric  predictions  of  the  International  Measurement  System  (IMS) 
and  PCSail  VPPs  to  small  model  tow  tank  testing  and  to  the  FKS  and  SPLASH  CFD  codes.  The 
IMS  VPP  results  were  supplied  by  the  Mk  II  Navy  44  designer,  David  Pedrick.  The  PCSail  VPP 
code  was  developed  by  David  Martin  and  Robert  Beck  from  the  University  of  Michigan  and  was 
distributed  at  the  15*’^  Chesapeake  Sailing  Yacht  Symposium.  The  tank  testing  was  conducted  in 
the  120'  tow  tank  at  the  United  States  Naval  Academy.  The  FKS  code  was  developed  at  the 
David  Taylor  Model  Basin,  Naval  Surface  Warfare  Center  -  Carderock,  by  Dr.  F.  Noblesse,  in 
cooperation  with  C.  Yang  of  George  Mason  University.  The  SPLASH  code  was  developed  by 
Bruce  Rosen  and  Joe  Laiosa  of  South  Bay  Simulations. 

One  important  aspect  of  this  study  was  to  compare  the  influence  of  turbulent  and  laminar 
flow  on  the  lifting  appendages  of  the  model.  The  model  has  a  waterline  length  of  45"  with  the 
keel  having  a  span  of  less  than  6"  and  the  rudder's  span  of  less  than  2".  The  model's  length  was  a 
function  of  the  size  of  the  tooling  equipment  in  the  model  shop  and  the  size  of  the  tow  tank. 
Furthermore,  both  CFD  codes  were  inviscid  -  the  frictional  component  of  the  fluid  was  not 
calculated  in  the  initial  drag  prediction.  Both  codes  were  chosen  based  on  their  speed  of 
calculation  and  their  availability. 

The  testing  program  made  performance  predictions  of  the  sailboat  at  four  wind  speeds  and 
five  wind  angles  as  well  as  the  upright  condition  which  was  equivalent  to  the  vessel  under 
mechanical  propulsion.  The  IMS  VPP,  Tank  Testing,  FKS,  and  SPLASH  evaluated  the  boat  in 
each  performance  condition.  The  PCSail  VPP  evaluated  the  boat’s  initial  sailing  characteristics 
in  each  condition. 
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An  initial  performance  prediction  process  shown  in  Figure  (1)  was  developed  through  this 
study.  The  hull  lines  and  rig  diagram  of  the  Mk  II  Navy  44  STC  were  the  initial  data  necessary 
for  all  calculations.  From  these  diagrams  of  the  Navy  44,  the  PCSail  VPP  was  able  to  compute 
the  sailing  attitudes  and  velocities  of  the  Navy  44.  These  positions  were  used  to  determine 
which  conditions  would  be  tested  in  the  tank.  The  tank  testing  gave  both  upright  and  sailing 
hydrodynamic  forces  as  output. 

As  in  PCSail,  the  diagrams  of  the  Navy  44  were  the  only  necessity  to  run  the  FKS  CFD 
code.  The  primary  output  from  FKS  was  upright  hydrodynamic  forces.  Sailing  hydrodynamic 
conditions  were  evaluated,  but  the  accuracy  of  FKS  in  evaluating  sailing  conditions  was  poor. 

In  order  for  the  SPLASH  CFD  to  be  more  accurate,  some  of  the  tank  testing  data  was  used 
as  input  to  SPLASH.  SPLASH  made  predictions  as  a  virtual  tow  tank,  being  able  to  predict  both 
upright  and  sailing  hydrodynamic  forces.  Finally,  all  of  these  hydrodynamic  forces  were  used  by 
a  custom-made  VPP  for  the  Mk  II  Navy  44.  The  output  of  this  custom  VPP  was  the  optimum 
sailing  attitudes  and  velocities  of  the  Mk  II  Navy  44. 
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PCSail  Velocity  Prediction  Program 

VPPs  operate  by  the  principle  that  a  sailing  vessel  is  in  equilibrium.  Therefore,  the 
summation  of  forces  and  moments  on  the  vessel  must  be  zero.  In  this  way,  a  VPP  is  divided  into 
two  parts:  hydrodynamic  and  aerodynamic.  The  hydrodynamic  forces  are  a  hull  form’s  reaction 
to  the  vessel  being  driven  through  the  water.  These  forces  include  hydrodynamic  drag,  side 
force,  yawing  moment,  and  righting  moment  (Figures  (2),  (3)).  The  hydrodynamic  forces  resist 
the  aerodynamic  forces  on  the  boat  due  to  sails.  These  forces  are  aerodynamic  drive,  side  force, 
yawing  moment,  and  heeling  moment.  When  each  of  the  hydrodynamic  forces  equals  the 
opposite  aerodynamic  forces,  the  boat  is  in  equilibrium  and  is  considered  a  possible  sailing 
condition.*  After  a  VPP  solves  for  possible  sailing  conditions,  the  VPP  finds  the  sailing 
conditions  which  are  fastest  around  a  course  and  considers  them  the  optimum  conditions. 


Hydrodynamic  Side  Force 


Figure  2  -  Z-plane  force  and  moments  diagram  Figure  3  -  X-plane  force  and  moments  diagram 

PCSail  was  supplied  as  freeware  at  the  15*  annual  Chesapeake  Sailing  Yacht  Symposium. 
PCSail  is  a  spreadsheet  for  Microsoft  Excel®  which  predicts  the  forces  on  the  boat  and  the 
optimum  sailing  conditions. 

PCSail  predicts  the  hydrodynamic  forces  on  a  vessel  from  the  Delft  series  of  sailing  yacht 
model  tests.  The  Delft  series  is  a  set  of  hydrodynamic  data  taken  from  a  series  of  different 


'  Longitudinal  trimming  moments  are  typieally  ignored  in  a  VPP  since  its  calculations  have  limited  effects  on 
velocity.  In  the  same  way,  Sinkage  and  buoyant  forces  are  also  assumed  to  have  negligible  effects  on  velocity. 
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sailing  hull  forms.  By  giving  key  hull  parameters,  PCSail  interpolates  through  the  Delft  series  to 
predict  the  hydrodynamic  qualities  of  a  given  hull  form. 

The  Hazen  method  documented  by  Larsson  and  Eliasson  calculates  the  aerodynamic  forces 
generated  by  the  sailing  rig.  The  characteristic  lift  and  drag  of  sails  for  PCSail  were  found  via 
experimental  data  as  well  as  theoretical  calculations.  After  giving  sail  shape  inputs,  PCSail 
predicts  the  lift  and  drag  on  the  sails  specific  to  a  vessel.  Finally,  lift  and  drag  in  the  sail-axis  are 
corrected  to  aerodynamic  drive  and  lift  in  the  vessel-axis. 

PCSail  only  solves  to  equate  aerodynamic  drive  to  hydrodynamic  drag  and  heeling  moment 
to  righting  moment.  Forces  which  are  ignored  included  hydrodynamic  side  force  and  yawing 
moments.  Aerodynamic  side  force  is  considered  only  for  calculating  heeling  moment.'^^^ 


PCSail  Setup 

Designer  David  Pedrick  supplied  the  IMS  VPP  predictions  for  the  Mk  II  Navy  44'  STC. 
The  IMS  VPP  serves  as  the  standard  of  performance  prediction  and  a  rating  tool  for  racing 
yachts.  However,  since  the  IMS  VPP  output  did  not  give  sailing  characteristics  such  as  sail 
forces  and  heel  angles,  the  PCSail  VPP  was  used.  The  boat  characteristics  from  the  IMS  VPP 
were  used  as  inputs  to  the  PCSail  VPP. 

The  initial  inputs  into  the  PCSail  program  were  very  straightforward  as  in  Figure  (4).  Some 
modifications  had  to  be  made  to  the  program  so  that  the  IMS  VPP  and  PCSail  VPPs  agreed  in 
their  output. 


Figure  4  -  PCSail  sample  hull  characteristics  input 


14 


Modifications  to  PCSail 

One  modification  was  that  of  GM,  the  metacentric  height.  The  metacenter  is  a  calculation 
which  helps  define  initial  stability  in  a  vessel.  PCSail  computed  GM  as: 


GM= 


(1) 


While  acceptable  for  trend  studies,  this  approximation  did  not  provide  the  righting  moment 
accuracy  within  the  tolerance  needed  for  this  study.  The  hydrostatics  program  FASTSHIP® 
evaluated  the  stability  of  the  Mk  II  Navy  44  STC.  GM  was  calculated  for  use  in  the  VPP,  and 
4.94  ft  was  used  instead  of  the  PCSail  calculated  GM  of  4.27  ft.  This  indicated  that  the  vessel 
was  more  stable  than  predicted  by  the  trend  from  PCSail. 


Secondly,  yaw  angle  and  rudder  angle  determinations  were  added  to  the  VPP.  PCSail  had 
already  calculated  the  total  aerodynamic  lift  from  each  sail.  From  the  supplied  rig  diagrams,  the 
longitudinal  center  of  effort  of  each  sail  was  approximated  at  the  geometric  center  of  the  sail. 
From  these  two  values,  the  yawing  moment  on  the  boat  due  to  aerodynamic  forces  was  found. 
From  the  lines  plan,  the  center  of  effort  on  each  appendage  was  approximated  so  that  the  lift 
force  required  on  each  appendage  could  be  found.  From  these  values,  the  necessary  angle  of 
attack  of  each  appendage  was  found. 

Another  small  modification  was  that  the  VPP  used  the  average  freeboard  height  as  the 
height  of  the  base  of  the  mast.  PCSail  initially  calculated  that  the  total  height  of  the  rig  was 
approximately  6  inches  less  than  the  true  height.  This  made  the  righting  moments  and  wind 
speeds  at  the  height  of  the  sails  lower  than  they  should  have  been.  The  height  of  the  base  of  the 
mast  was  modified  to  reflect  the  true  base  mast  height. 


VPP  Results 

The  VPP  was  run  at  wind  speeds  ranging  from  6  to  24  knots.  The  PCSail  output 
corresponded  extremely  well  with  most  of  the  IMS  VPP  output.  There  was  a  fairly  consistent 
discrepancy  at  low-wind  speeds.  This  was  expected  since  the  PCSail  aerodynamic  algorithm  is 
not  as  rigorous  as  the  IMS  algorithm  at  these  low  wind  speeds.  Figure  (5)  shows  this 
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discrepancy  in  light  air  especially  for  heel  angles  at  reaching  angles  (120-170  degrees)^.  This 
sailing  condition  is  very  sensitive  to  apparent  wind  effects  as  the  sails  are  often  stalled  and 
requires  more  complex  wind  calculations  for  more  accurate  results.  However,  most  of  the  range 
had  to  be  checked  against  their  respective  velocities-made-good  (VMG),  which  is  the  speed  of 
the  vessel  on  the  wind-axis.  For  instance,  if  leeway  is  neglected  at  a  true  wind  direction  of  90 
degrees,  the  VMG  is  zero.  By  choosing  the  maximum  absolute  value  of  velocity-made-good,  the 
wind  angle  at  which  a  boat  can  most  optimally  travel  upwind  or  downwind  can  be  found.  Most 
of  the  deep-reaching  angles  were  at  angles  past  the  best  VMG,  and  therefore  were  out  of  the 
range  of  the  optimum  sailing  angle. 


Figure  5  -  PCSail  VPP  compared  to  IMS  VPP  at  6-knots  wind  speed 


The  PCSail  VPP  seemed  to  disagree  with  the  IMS  VPP  for  heel  around  a  true  wind  angle  of 
90  degrees.  This  disagreement  was  based  on  the  two  algorithm’s  different  interpretations  on 
when  to  raise  and  lower  the  jib  and  spinnakers.  Therefore,  the  heel  angle  was  not  a  continuous 
function  of  true  wind  angle  since  a  spinnaker  and  jib  have  much  different  effects  on  heel. 
Velocity  for  jib  and  spinnaker,  however,  was  nearly  equal  at  the  point  which  the  VPP  chose  to 
switch  sails.  Velocity  was  considered  continuous  as  a  function  of  true  wind  angle  even  while 
performing  sail  changes. 


^  Wind  angles  are  reported  as  the  true  wind  angle  aft  of  the  bow.  The  bow  is  an  angle  of  0  degrees,  wind  from  the 
beam  is  at  90  degrees,  and  wind  from  the  stem  is  at  180  degrees. 
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For  medium  to  high  wind  speeds  as  in  Figure  (6),  the  PCSail's  data  came  very  near  the  IMS 


VPP's  data.  Even  heel  angle  showed  much  better  correlation  between  the  two  VPPs. 


Velocity  as  a  Function  of  True  Wind  Angle  at 
12-knots  True  Wind  Speed 


Heel  Angle  as  a  Function  of  True  Wind  Angle  at 
12-knots  True  Wind  Speed 


0  30  60  90  120  150  180 

True  Wind  Angle  [deg] 


Figure  6  -  PCSail  VPP  compared  to  IMS  VPP  at  12-knots  wind  speed 
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Tank  Testing 

Tank  testing  is  the  traditional  method  of  predicting  the  performance  of  vessels.  Since  it  is 
far  too  expensive  to  build  full-size  ships  to  test  different  hull  forms,  models  are  built  for  testing 
instead. 


Performance  Prediction  Theory 

Performance  prediction  tank  testing  accomplishes  one  goal:  determining  the  forces  on  the 
hull  at  different  attitudes  and  speeds.  To  accurately  predict  forces,  the  model  and  ship  need 
dynamic  similarity.  There  are  three  parts  to  dynamic  similarity:  scaling  of  model  sizes,  model 
speeds,  and  model  forces. 

Geometric  similarity  is  the  linear  scaling  of  the  model  to  the  ship  (Eqn  (2)).  This  is  met 
when  all  characteristic  lengths  of  the  model  (Lmodei)  are  made  proportional  to  the  ship  (Lship)  by 
the  scaling  factor 


L. 


'^ship  _  2^ 


(2) 


'model 


Kinematic  similarity  is  the  scaling  of  speeds  (Eqn  (3)).  Kinematic  similarity'^'^^  requires  that 
the  ratio  of  velocities  of  the  model  (Vmodei)  and  ship  (Vship)  are  the  same.'^'"^^ 


^ship  _  ^ship  _ 


model 
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'model 


(3) 


Dynamic  similarity  is  the  scaling  of  forces  (Eqn  (4)).  In  the  tank,  hull  resistance  is  the  force 


of  the  water  acting  on  the  ship.  All  forces  (F)  can  be  decomposed  into  their  non-dimensional 
coefficient  (C).^'^^ 


C  = 


F 

YipSV^ 


(4) 


Besides  force,  a  coefficient  is  dependent  on  the  mass  density  of  the  fluid  (p),  the  velocity 


(V)  of  the  fluid  relative  to  the  characteristic  length  (L)  of  the  body,  and  the  characteristic  wetted 
surface  area  of  the  hull  (S). 


Total  resistance  coefficient  (Ct)  is  the  summation  of  viscous,  wavemaking,  air,  and  added 
wave  resistances  (Eqn  (5)).  The  primary  resistances  are  viscous  (Cy)  and  wavemaking  (Cw),  and 
the  others  are  omitted  for  simplicity. 
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c  =c  +c 


The  viscous  resistance  coefficient  (Cy)  is  a  function  of  Reynolds  number  (Re)  and  the 
frictional  resistance  coefficient  (Cf).  Reynolds  number  defines  the  speed  of  the  vessel  relative  to 
the  condition  of  fluid  flow.  It  is  the  ratio  of  inertia  forces  to  viscous  forces  (Eqn  (6)).'^"^^ 

Re=—  (6) 

V 

Reynolds  number  is  dependent  on  the  velocity  of  the  fluid,  the  characteristic  distance  of 
fluid  flow  over  the  body  (L),  and  the  kinematic  viscosity  of  the  water  (v). 

There  is  one  modification  to  the  Reynolds  number  calculation  which  is  exclusive  to  sailing 
craft.  For  normal  non-sailing  vessels,  the  flow  over  the  body  (L)  is  typically  the  length  of  the 
waterline.  For  sailing  craft,  the  keel  and  rudder  make  a  large  difference  in  frictional  resistance. 
Moreover,  the  local  Reynolds  number  for  both  keel  and  rudder  are  much  less  than  the  Reynolds 
number  across  the  entire  hull.  Therefore,  the  flow  conditions  across  much  of  the  wetted  surface 
area  for  sailing  craft  is  at  a  lesser  Reynolds  number  than  if  the  whole  body  was  the  characteristic 
flow  length.  The  result  is  an  approximation  to  Reynolds  number  (Rn)  for  sailboats:  a  factor  of 
0.7  is  added  to  the  formula  for  Reynolds  number  (Eqn  (7)). 

R,  =0.7—  (7) 

V 

The  frictional  resistance  coefficient  (Cf)  can  be  closely  approximated  without  modeling  by 
using  the  standard  1957  International  Towing  Tank  Conference  (ITTC)  ship-model  correlation 
equation  (Eqn  (8)).  The  frictional  resistance  equation  is  a  function  created  through  the  curve  fits 
of  many  years  of  experimental  data.'^'"^^ 

0.075 

The  viscous  resistance  equation  includes  the  skin  friction  term  and  a  form  factor  (k).  The 
form  factor  accounts  for  added  resistance  from  the  shape  of  the  vessel.  This  reduces  to  an 
equation  for  the  viscous  resistance  coefficient  (Eqn  (9)).'^'"^^ 


C^=(l  +  A:)-C^ 


(9) 
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The  wavemaking  resistance  coefficient  (Cw)  is  a  function  of  the  non-dimensional  Froude 
number  (Fn).  Froude  number  defines  the  speed  of  the  vessel  relative  to  the  type  of  wave  systems 
created.  It  is  the  ratio  of  inertial  forces  to  gravity  forces  (Eqn  (10)). 


F  = 


V 


(10) 


Froude  number  is  dependent  on  the  velocity  of  the  fluid,  the  acceleration  due  to  gravity  (g), 
and  the  characteristic  distance  of  fluid  flow  over  the  body. 


There  is  no  equation  for  the  wavemaking  resistance  coefficient,  and  it  remains  a  function  of 
only  Froude  number  (Eqn  (1 1)).'^'^^ 


C„=/(F.) 


(11) 


Since  it  is  assumed  that  only  viscous  and  wavemaking  resistances  are  significant,  and  there 
exists  an  equation  for  the  viscous  resistance  coefficient,  an  equation  of  resistance  can  be  defined 
in  Eqn  (12).  Wavemaking  becomes  the  only  unknown  variable  in  the  resistance  equation  and  is 
found  experimentally  through  tank  testing  or  CFD. 


Rt=-pSV^ 


(l  +  A:) 


0.075 


_(loglO(i?J-2)^ 


+  C„ 


(12) 


Hydrodynamic  Fluid  Flow  Theory 

The  flow  condition  around  a  ship  is  composed  of  three  basic  parts:  turbulent  flow,  laminar 
flow,  and  separated  flow. 
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The  aerodynamicist,  Ludwig  Prandtl,  developed  the  eoncept  of  the  boundary  layer  for  fluid 
flow  (Figure  (7)).  Flow  in  the  region  surrounding  a  body  can  be  divided  into  the  boundary  layer 
and  anywhere  outside  of  the  boundary  layer.  Outside  of  the  boundary  layer,  the  viscosity  on  the 
body  is  negligible,  and  the  flow  may  be  treated  as  inviscid.  Therefore,  only  flow  inside  the 
boundary  layer  is  important  for  determining  frictional  resistance. 

The  laminar  boundary  layer  can  be  assumed  to  have  a  parabolic  velocity  profde  “U”  as 
shown  in  Figure  (7).  At  the  body,  the  velocity  of  the  fluid  is  zero.  Reynolds  number,  pressure 
gradient,  and  surface  roughness  determine  the  length  of  the  laminar  boundary  layer  and  disturbed 
flow.  For  a  flat,  smooth  plate  in  undisturbed  flow,  the  flow  typically  moves  from  laminar  flow  to 
the  transition  zone  around  Rn  =  500,000.'^'^  A  description  of  laminar  flow  is  that  it  is  unmixable, 
meaning  that  if  the  flow  were  to  oscillate  between  positive  and  negative  flow  directions,  a  fluid 
particle  would  always  remain  on  the  same  streamline.'^^^  The  coefficient  of  skin  friction  on  a  flat 
plate  in  laminar  flow  is: 


_  0.730 


(13) 


Fluctuating  velocities  and  eddying  motion  characterize  the  turbulent  boundary  layer. 
Compared  to  the  laminar  boundary  layer,  the  turbulent  boundary  layer  is  thicker,  and  the  velocity 
profde  of  the  flow  increases  much  faster  away  from  the  surface.  However,  before  the  flow  enters 
into  the  turbulent  region,  it  passes  through  the  transition  zone  (Figure  (7)).  The  transition  zone 
contains  sudden  spots  of  turbulence  generation  and  a  general  breakdown  of  laminar  flow. 
Introduction  into  the  transition  zone  is  primarily  a  function  of  Reynolds  number,  but  surface 
roughness  and  disturbed  flow  can  also  trip  the  condition  from  laminar  to  turbulent.  The 
coefficient  of  skin  friction  on  a  flat  plate  in  turbulent  flow  is: 


_  0.0576 


(14) 


A  special  region  inside  the  turbulent  boundary  layer  is  called  the  viscous  sublayer  (Figure 
(8)).  The  viscous  sublayer  is  mainly  laminar  flow,  but  for  a  typical  sailboat,  is  only  0.1  mm 
thick.  The  importance  of  the  layer’s  size  is  that  it  is  affected  by  surface  roughness,  and  the  layer 
is  liable  to  trip  back  to  the  turbulent  boundary  layer  in  enough  disturbance. When  laminar,  the 
viscous  sublayer  significantly  decreases  the  skin  friction  on  the  vessel. 
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Turbulent  Boundary  Layer 


Figure  8  -  Exaggerated  flow  conditions  on  a  typical  sailing  vessel^^^ 


The  separated  zone  is  typically  found  near  the  stem  of  a  vessel  and  is  a  function  of  the 
pressure  gradient.  For  a  body  that  is  not  a  flat  plate,  the  pressure  gradient  will  remain  positive  as 
the  form  of  the  body  increases  (from  the  bow  to  after  amidships  on  a  typical  sailboat).  However, 
as  the  form  of  the  body  bends  inward,  the  velocity  profde  will  still  try  to  cling  to  the  surface. 
The  resultant  is  that  the  flow  will  invert  under  this  negative  pressure  gradient.  The  inverted  flow 
will  cause  large  eddies  much  larger  than  the  eddies  found  in  normal  turbulent  flow.  Separated 
flow  drag  is  significantly  higher  than  turbulent  flow  drag. 

For  tank  testing,  the  flow  condition  is  a  function  of  Reynolds  number,  the  roughness  of  the 
surface  of  the  hull,  and  the  turbulence  in  the  surrounding  water.  For  the  full-size  Mk  II  Navy  44 
STC,  the  flow  condition  along  most  of  the  hull  is  turbulent  as  shown  in  Figure  (8).  The  flow 
condition  around  most  of  the  model  would  be  laminar  since  the  model  is  small,  the  tank  water  is 
generally  stable,  and  the  model  surface  is  very  clean.  Although  a  model  can  be  scaled 
geometrically  from  a  ship,  the  water  and  therefore  the  flow  condition  cannot.  If  resistance  data 
with  a  model  in  mostly  laminar  flow  were  scaled  to  full-scale,  the  predicted  resistance  would  be 
significantly  lower  than  actual. 

Keeping  Reynolds  and  Froude  numbers  constant  defines  dynamic  similarity  in  terms  of 
viscous  and  wavemaking  resistance  respectively.'^"^^  Therefore,  if  both  Reynolds  and  Froude 
numbers  of  the  model  matched  those  of  the  ship,  then  the  model’s  behavior  would  be  scaled 
identically  to  the  ship’s  behavior.  However,  because  the  acceleration  due  to  gravity  and  the 
kinematic  viscosity  of  water  cannot  be  scaled  at  the  same  time,  perfect  modeling  conditions 
cannot  exist.  In  some  cases,  Reynolds  number  similarity  is  used  in  a  wind  tunnel  for  testing 
submerged  appendages.  For  most  hydrodynamic  tank  testing,  however,  Reynolds  number 
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similarity  is  neglected  because  it  would  require  the  model  to  travel  at  speeds  much  greater  than  a 
tow  tank  could  provide.  Froude  number  similarity  is  used  because  the  model’s  speed  is  slower 
than  the  ship’s  speed.  Furthermore,  viscous  resistance  is  fairly  well  predicted  using  skin  based 
methods  (such  as  the  given  ITXC  equation)  while  the  wavemaking  resistance  is  generally  more 
sensitive  to  the  details  of  the  hull  form  and  thus  must  be  determined  for  each  new  design.  Using 
the  Froude  hypothesis,  the  forces  on  the  model  can  be  predicted  for  the  full-scale  ship.'^*"^^ 

The  Froude  hypothesis  states  that  the  wavemaking  coefficient  is  the  only  variable  which  can 
be  directly  equated  between  the  ship  and  model.  In  accordance,  for  each  test  performed  in  the 
tank,  the  Froude  number  of  the  extrapolated  ship’s  speed  is  equal  to  the  Froude  number  of  the 
model.  Extrapolation  is  begun  by  decomposing  the  tank  test  data  into  the  coefficient  of 
wavemaking.  Thereafter,  the  Reynolds  number  as  well  as  the  coefficient  of  viscous  resistance  is 
recalculated  for  the  ship.  The  new  viscous  resistance  term,  the  wavemaking  term,  and  a 
correlation  allowance  are  summed  to  produce  a  new  coefficient  of  total  resistance  for  the  ship. 
The  correlation  allowance  mainly  takes  into  effect  the  discrepancy  between  the  smooth-painted 
hull  of  the  model  and  the  rough,  biologically  coated  hull  of  the  ship.  Typically,  the  correlation 
allowance  is  0.0004,  although  this  varies  in  accordance  with  the  assumed  roughness  of  the  ship’s 
hull. 

Once  the  coefficient  of  total  resistance  for  the  ship  is  calculated,  the  total  resistance  force  is 
calculated,  and  the  ship  prediction  is  completed.  According  to  dimensional  analysis,  geometric 
force  scaling  can  be  performed  by  multiplying  the  model  force  by  the  cubed  scaling  factor.  If 
this  is  done  for  resistance  tests,  the  ship  prediction  will  be  over-estimated.  The  reason  for  this  is 
that  although  there  is  a  greater  amount  of  viscous  resistance  for  the  ship  due  to  its  larger  wetted 
surface  area,  the  coefficient  of  viscous  resistance  is  smaller.  Turbulent  flow  theory  predicts 
smaller  coefficients  of  viscous  resistance  for  larger  Reynolds  numbers.  Since  the  ship  has  a 
much  larger  length  and  travels  at  larger  velocities,  the  Reynolds  number  for  the  ship  is  larger 
than  for  the  model.  The  ITTC  resistance  equation  which  is  used  in  this  extrapolation  reflects  the 
differences  in  the  viscous  coefficient. 


Model  Preparation 
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The  tank  used  in  testing  the 
Mk  II  Navy  44  STC  was  the  120’ 
tow  tank  loeated  in  USNA’s 
Hydromechanics  Laboratory.  The 
towing  rig  on  the  120’  tank 
required  extensive  preparation. 
The  attitude  of  the  model  needed 
to  be  set  in  almost  any  condition 
of  heel,  yaw  angle,  trim,  heave, 
and  rudder  angle  (Figure  (9)). 


The  requirements  on  the  rig  were  therefore: 

1)  Infinitely  adjustable  in  yaw  from  0  to  10  degrees.  Yaw  is  the  angle  at  which  the  bow 
points  when  traveling  through  the  water. 

2)  Infinitely  adjustable  in  heel  from  0  to  35  degrees  (or  deck  submersion).  Heel  is  the  angle  at 
which  the  boat  rolls  around  the  longitudinal  axis  of  the  boat. 

3)  Infinitely  adjustable  in  rudder  angle  from  0  to  10  degrees. 

4)  Infinitely  adjustable  in  initial  trim  until  deck  submersion.  Trim  is  an  angular  measurement 
of  how  either  the  bow  or  stem  is  submerged  relative  to  a  lateral  axis  through  midships. 

5)  Free  to  trim. 

6)  Free  to  heave.  Heave  is  the  vertical  movement  of  the  entire  boat  into  or  out  of  the  water. 

Because  of  the  forces  and  attitudes  specific  to  sailing  craft,  more  instmmentation  than  just 
drag  was  needed.  For  the  Navy  44,  the  recorded  data  included: 

1)  Drag  parallel  to  motion  of  travel 

2)  Lift  perpendicular  to  motion  of  travel 

3)  Yawing  moment 

4)  Trae  speed 

5)  degrees  of  change  in  trim 

The  model  was  received  from  the  USNA  Model  Shop  built  to  the  lines  supplied  by  Pedrick 
Yacht  Design.  The  model  was  measured  and  put  to  a  position  of  0  degrees  heel  and  trim  to  find 
midships  and  the  centerline.  A  small  aluminum  plate  was  mounted  to  be  flat  at  this  baseline 
upright  condition.  Lines  of  centerline  and  midships  were  etched  into  the  plate. 


Figure  9  -  Degrees  of  freedom  in  model  testing 
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A  bulkhead  was  installed  in  the 
model  a  few  inches  forward  of 
midships.  Attached  to  the  aft  side  of 
the  bulkhead  were  two  aluminum 
plates  which  could  be  adjusted  so  that 
the  model  was  set  in  a  position  of 
heel.  Forward  of  the  aft  plate  was  a 
cylinder  of  bearings  which  housed 
the  heave  post  and  left  the  model  free 
to  trim  (Figure  (10)). 


A  pin  was  attached  to  the 
transom  on  centerline.  This 
connected  to  a  stiff  carbon  fiber  rod 
which  could  be  extended  and  was 
itself  connected  to  the  tow  rig.  The 
rod  allowed  the  model  to  be  set  in  a 
yawed  position  (Figure  (11)). 

Figure  11  -  Yaw  adjuster 

A  4"  tiller  was  attached  to  the 
rudderpost.  At  the  end  of  the  rudder 
was  a  pin  which  connected  to  an 
electronic  caliper.  The  other  end  of 
the  caliper  was  pinned  to  the  hull. 

The  caliper  allowed  for  extremely 
precise  measurements  in  rudder  angle 
(Figure  (12)). 


Figure  12  -  Rudder  adjuster 
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Next,  three  force  blocks  and  two 
inclinometers  were  calibrated. 
Force  blocks  measured  resistance  in 
the  direction  of  travel  (drag),  lift 
perpendicular  to  motion  at  the  heave 
post  (fwd-lift),  and  lift 
perpendicular  to  motion  at  the  aft 
pin  (aft-lift).  The  directions  of  these  forces  are  shown  in  Figure  (13).  Inclinometers  would 
measure  change  in  trim  and  heel.^  The  aft-lift  block  would  be  used  to  measure  the  longitudinal 
(yawing)  moment.  All  force  blocks  were  calibrated  using  standard  weights  up  to  twice  the 
expected  loads.  The  inclinometers  were  calibrated  using  a  stand-alone  electronic  inclinometer 
on  an  adjustable  surface.  Both  inclinometers  were  calibrated  to  their  rated  maximum  incline. 

The  drag  and  lift  force  blocks  were  both  rated  to  a  25  lb  capacity.  This  was  to  minimize  the 
effect  of  cross-moments  interfering  with  the  force  acquisition.  The  aft-lift  block  did  not  have 
any  interfering  moments  and  was  rated  at  5  lbs.  In  the  calibration  for  all  force  blocks,  the 
maximum  linear  deviation  was  0. 1 16%  on  the  5  lb  aft- lift  block.  This  corresponded  to  a  standard 
deviation  of  0.0037  lbs,  which  was  considered  very  acceptable. 

The  heave  post  was  attached  to  the  drag  block,  mounted  underneath  the  lift  block,  and 
finally  secured  to  a  square  aluminum  beam  on  the  carriage.  The  aft-lift  block  was  mounted  to  an 
"L"-shaped  beam,  which  was  clamped  onto  the  after  end  of  the  carriage.  The  inclinometers  were 
mounted  to  the  aluminum  plate  inside  the  model.  The  speed  of  the  carriage  had  been  previously 
calibrated  and  was  measured  at  the  motor's  gears. 

After  the  inclinometers  were  added,  the  model  was  weighed  and  trimming  weights  were 
added  in  order  to  scale  the  model  geometrically.  Placing  the  model  in  the  water,  the  weights 
were  moved  so  that  the  model  had  no  trim.  The  position  of  these  weights  marked  0-degrees  of 
trim,  and  a  ruler  was  drawn  on  the  centerline  to  measure  the  trimming  arm  induced  by  the 
weights.  Figure  (14)  shows  the  final  model  setup: 


Figure  13  -  Schematic  indicating  measured  forces 


^  The  heel  of  the  vessel  was  set  through  thumb  screws  mounted  to  the  bulkhead.  An  inclinometer  was  set  to 
measure  any  change  in  the  heel  of  the  vessel  -  which  would  not  be  desired.  The  acquired  data  from  the  tests  showed 
that  the  heel  did  not  change. 
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Figure  14  —  Model  ready  for  testing 


All  force  blocks  and 

inclinometers  had  wires  which  ran 
to  two  amplifiers  mounted  on  the 
tow  rig  (Figure  (15)).  The  data 
then  ran  from  the  amplifiers  to  an 
analog-digital  converter  and 

finally  through  a  USB  connection 
to  a  PC.  Before  gage  calibrations, 
the  amplifier  was  set  to  return 
values  in  the  linear  region  of 
output.  Figure  15  -  Amplifiers  mounted  on  the  tow  rig 
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The  final  stage  in  testing  was  to  find  the  position  of  0  degrees  of  yaw.  The  model  was 
attached  to  the  rig  without  the  rudder  and  with  0-degrees  heel  and  trim.  The  model  was  then  run 
down  the  tank  a  minimum  of  four  times  where  each  trial  varied  by  yaw.  Model  speed  was  held 
at  a  constant  3.2  ft/s.  This  speed  was  to  ensure  typical  flow  conditions  across  the  appendages, 
but  with  minimal  wavemaking  interference."^  Plots  were  constructed  of  inches  of  yaw  adjustment 
versus  pounds  of  side  force.  The  position  of  0  degrees  of  yaw  was  marked  at  the  point  of  no  side 
force. 

Next  the  rudder  was  added,  and  the  process  of  finding  0  degrees  of  rudder  was  identical  to 
the  above  process,  except  the  rudder  angle  was  adjusted  vice  yaw.  Again,  the  point  of  no  side 
force  marked  the  position  of  0  degrees  of  rudder.  This  position  of  0-rudder,  0-yaw,  0-trim,  and 
0-heel  was  considered  the  standard  upright  condition. 

The  typical  standard  deviation  in  calibrating  the  yaw  and  rudder  angles  was  0.019  lbs  as 
measured  by  both  lift  blocks.  These  calibrations  were  checked  every  time  the  model  was  taken 
out  of  the  water,  which  was  at  least  after  every  day  of  testing. 


At  higher  speeds,  the  wavemaking  effects  would  have  produced  large  waves,  causing  the  model  to  heave  and  trim 
substantially.  This  would  disrupt  the  measurements  taken  in  calibration. 
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Tank  Tests 

There  were  three  parts  of  tank  testing: 

1)  Unstimulated  Upright  Conditions 

2)  Stimulated  Upright  Conditions 

3)  Stimulated  Sailing  Conditions 

The  goal  of  unstimulated  testing  was  to  find  a  baseline  for  the  evaluation  of  the  later-added 
turbulence  stimulators. 

Stimulators  were  added  to  the  model  to  try  to  create  turbulent  flow  condition  around  the 
model.  These  stimulators  were  placed  at  a  position  to  trip  the  flow  where  turbulent  flow  would 
begin  on  the  ship.  Details  of  these  stimulators  are  discussed  in  the  “Stimulated  Upright 
Condition”  section.  Additional  stimulators  were  added  after  all  other  testing  was  completed  in 
order  to  factor  out  the  added  resistance  due  to  the  stimulators. 

The  stimulated  upright  condition  provided  results  which  represent  the  motoring  condition  of 
the  vessel.  The  upright  resistance  of  a  vessel  is  also  typically  used  to  judge  its  overall  resistance 
against  other  vessels.  In  the  upright  condition,  heel  angle,  yaw  angle,  and  rudder  angle  were  all 
set  at  zero  degrees.  The  vessel  was  trimmed  to  its  designated  waterline  and  was  left  free  to  trim 
and  heave  while  at  speed  due  to  hydrodynamic  forces. 

Stimulated  sailing  conditions  allowed  for  calculations  of  the  vessel  at  different  attitudes. 
For  sailing  conditions,  the  vessel  had  a  set  heel  angle,  yaw  angle,  rudder  angle,  and  trim  angle. 
The  heel,  yaw,  and  rudder  angles  were  set  initially  from  a  chosen  matrix.  Trim  angle  was  set 
through  trim  weights  which  approximate  the  aerodynamic  trimming  moment.  This  trimming 
moment  was  determined  experimentally  through  an  initial  run  for  each  condition  and  equaled  the 
countering  measured  hydrodynamic  trimming  moment.  The  model  was  free  to  heave  in  the 
sailing  condition. 
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Upright  Unstimulated  Tests 

The  unstimulated  testing  was  eondueted  at  the  standard  upright  eondition.  The  eoeffieient 
of  total  resistanee  is  shown  in  Figure  (16).  The  generated  eurve  maintained  eohereney  above 
Rn=400,000  signaling  turbulent  flow.  Below  Rn=400,000,  there  was  signifieant  disorder  in  the 
resistanee  eurve,  showing  that  the  flow  eondition  was  moving  from  laminar  through  the 
transition  zone. 
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Figure  16  -  Drag  in  upright  unstimulated  flow 
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Upright  Stimulated  Tests 

After  the  unstimulated  flow  testing  was  eompleted,  sand  strips  were  applied  to  the  hull  to  aet 
as  turbulent  stimulators.  Sand  strips  were  plaeed  on  the  hull  1"  aft  of  the  stem  and  on  the  keel 
and  rudder  at  25%  of  their  ehords’  girth.  A  grain  size  of  0.5mm  to  1.0mm  was  used  and  at  a 
width  of  a  0.25".  The  grain  size  was  ehosen  to  be  as  large  as  the  laminar  boundary  layer  in  order 
to  disrupt  the  flow  suffieiently.  This  disturbed  boundary  layer  should  ehange  from  laminar  to 
turbulent  flow.  The  upright  stimulated  testing  was  eondueted  at  the  standard  eondition. 
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Figure  17  -  Drag  in  upright  conditions 

The  results  of  both  unstimulated  and  stimulated  conditions  are  shown  in  Figure  (17).  The 
stimulated  condition  showed  a  uniform  increase  in  total  resistance.  The  added  resistance  was 
from  the  friction  caused  by  the  sand  strips  and  also  indicated  that  the  flow  condition  did  not 
substantially  change  at  Rn  >  400,000.  The  meaning  of  this  was  that  it  was  considered  safe  to 
assume  ship-like  turbulent  flow  above  Rn  =  400,000.  Stimulated  tests  performed  below 
Rn=400,000  showed  that  the  fluid  flow  was  still  in  the  transition  zone  on  at  least  some  part  of  the 
vessel. 
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Additional  sand  strips  were  added  to  the  model  to  ealeulate  the  added  frietion  from  the 
strips.  The  first  set  of  additional  strips  doubled  the  width  of  the  previous  strips  and  was  loeated 
aft  of  and  adjaeent  to  the  previous  strips.  The  seeond  set  of  strips  brought  the  width  up  to  three 
times  the  width  of  the  original  strips.  The  third  set  of  strips  was  aft  of  and  adjaeent  to  the 
seeondary  strips.  Testing  of  the  strips  populated  only  the  lower  Reynolds  number  eonditions 
sinee  the  strips  were  used  only  for  Prohaska  analysis  of  the  data.^ 
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Figure  18  -  Drag  in  the  upright  conditions  with  additional  sand  strips 

Again,  the  additional  sand  strips  show  a  general  increase  in  resistance  (Figure  (18)).  The 
trend  in  the  resistance  curve  was  maintained  after  Rn  =  300,000.  Below  Rn  =  300,000, 
separation  between  the  second  and  third  sand  strip  increased.  This  showed  the  presence  of 
laminar  flow  even  with  the  extra  turbulence  stimulators.  Because  of  this,  results  below  Rn  = 
300,000  with  additional  sand  strips  were  questionable. 
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^  Prohaska  analysis  will  be  discussed  in  the  “Prohaska  Theory”  section. 
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Stimulated  Sailing  Tests 

The  sailing  conditions  were  estimated  from  attitudes  and  conditions  predicted  by  the  PCSail 
VPP.  This  matrix  varied  by  speed,  heel  angle,  yaw  angle,  and  rudder  angle. 

PCSail  was  run  at  true  wind  speeds  of  6,  12,  18,  and  24  knots  and  at  40,  60,  90,  120,  and  170 
degrees  to  the  wind.  The  VPP  determined  the  values  of  boat  speed  and  heel  angle  for  each  of 
these  20  conditions.  The  values  of  speed  and  heel  are  shown  in  Figure  (19). 
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Figure  19  -  Sailing  tests  matrix 

The  points  picked  for  positions  of  speed  and  heel  were  near  to  the  VPP’s  prediction.  These 
points  were  picked  so  that  the  true  sailing  conditions  could  be  interpolated  from  the  matrix  once 
the  testing  had  concluded.  Furthermore,  the  speed  of  the  carriage  could  be  set  only 
approximately  at  its  input  dial.  Therefore,  the  desired  speed  of  each  test  differed  slightly  from 
the  resulting  speed  of  the  carriage. 

Each  sailing  test  had  a  predicted  yaw  and  rudder  angle  from  calculations  added  to  the  PCSail 
VPP.  However,  as  the  sailing  testing  was  conducted  at  first  in  unstimulated  flow,  it  was  found 
that  the  aft-lift  block  was  recording  negative  forces  on  the  model  which  was  opposite  of  the 
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Figure  19  -  Sailing  tests  matrix 
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predicted  aerodynamic  force.  The  modifications  to  the  PCSail  VPP  had  estimated  the  angle  of 
attack  on  the  keel  and  rudder  from  their  calculated  coefficients  of  lift,  while  the  lifting  force  due 
to  the  canoe  body  was  assumed  zero.  The  PCSail  modifications  to  find  rudder  angle  and  keel 
angle  proved  to  be  incorrect.  This  was  due  either  to  the  lift  caused  by  the  canoe  body  itself  or 
from  unusual  flow  patterns  at  the  rudder  possibly  due  to  laminar  separation. 

The  result  of  this  finding  was  two-fold.  First,  the  unstimulated  testing  was  adjusted  to 
include  only  the  upright  condition.  Secondly,  the  test-matrix  was  reevaluated  so  that  true  sailing 
conditions  could  be  found. 

A  test  was  performed  to  find  how  yaw  affected  the  boat's  characteristics.  From  this  test  of 
yaw  shown  in  Figure  (20),  it  was  determined  that  for  small  angles  of  changes  in  yaw  the  lift 
would  change  linearly.  In  the  new  test  matrix,  instead  of  trying  to  predict  the  yaw  and  rudder 
angles,  set  angles  were  chosen.  Yaw  angles  of  0  and  4  degrees  and  rudder  angles  of  0,  3,  and  6 
degrees  were  used  for  each  sailing  condition.  Ultimately,  for  the  sailing  matrix  with  20 
variations  of  speed  and  heel,  there  were  a  total  of  120  tests  which  included  rudder  and  yaw  angle. 
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Figure  20  -  Lift  linearity  for  small  changes  in  yaw 
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Prohaska  Theory 

The  total  resistance  is  assumed  to  be  composed  of  viscous  and  wavemaking  factors  (Eqn 
(12)).  Viscous  resistance  is  composed  of  frictional  resistance  multiplied  by  a  form  factor  (Eqn 
(9)).  Because  of  Bernoulli’s  principle,  an  incompressible  fluid  has  to  change  velocity  while 
moving  across  a  form  which  is  not  flat.  The  increased  velocity  of  the  fluid  creates  increased 
frictional  resistance  across  the  form,  and  that  component  is  calculated  in  the  viscous  resistance 
term  as  the  form  factor, 


To  calculate  the  form  factor,  a  process  similar  to  the  Prohaska  method  is  used.  Prohaska 
suggested  that  the  form  factor  could  be  reduced  from  the  equation: 


(15) 


The  International  Tank  Testing  Committee  in  1978  recommended  modifying  Prohaska’ s 
method  to  use  the  equation: 


where  n  is  a  power  of  Froude  number  being:  4  <  «  <  6 . 


(16) 


Typically,  a  plot  is  constructed  where  Cj  jCp  are  ordinal  values  and  F„"/c^  are  abscissa 

values.  The  form  factor  is  derived  as  an  ordinal-intercept  of  the  best-fit  linear  regression. 
Values  used  for  computation  should  be  at  as  slow  a  speed  as  possible  without  reaching  into  the 
scale-effects  region. 
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Prohaska  Analysis 

The  unstimulated  and  stimulated  upright  tests  were  analyzed  using  Prohaska  plots  to  find  the 
form  factor  of  the  model. 


Figure  21  -  Prohaska  plot  of  unstimulated  data 

The  data  which  showed  linearity  without  scale  effects  was  used  for  the  Prohaska  plots  in 
Figure  (21).  The  best-fit  lines  produced  ordinal-intercepts  of  1.009,  1.035,  and  1.054  for  n 

values  of  4,  5,  and  6,  respectively.  The  jCp  Prohaska  plot  produced  the  best-fit  line,  so  the 


value  of  1+k  for  the  form  factor  in  unstimulated  flow  was  1.009. 
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The  same  method  was  used  for  stimulated  flow. 
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Figure  22  -  Prohaska  plot  of  stimulated  data 

The  data  was  sorted  to  find  good  points  for  the  Prohaska  plots  in  Figure  (22).  The  best-fit 
lines  produced  ordinal-intercepts  of  1.168,  1.181,  and  1.209  for  n  values  of  4,  5,  and  6, 

respectively.  Again,  the  jCp  Prohaska  plot  produced  the  best- fit  line,  so  the  value  of  1+k  for 

the  form  factor  in  stimulated  flow  was  0.168.  The  value  of  0.168  included  both  form  factor  and 
the  effect  of  one  strip  of  sand.  This  form  factor  was  used  to  compute  the  viscous  friction  from 
the  model  to  the  coefficient. 
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Finally,  the  sand-strip  stimulation  data  was  analyzed. 


Figure  23  -  Prohaska  plot  of  additional  sand-strips  data 

The  data  for  both  extra  sand  strips  were  analyzed  for  n  values  of  4,  5,  and  6  (Figure  (23)). 
The  jCp  Prohaska  plot  fit  the  best  line  through  both  sets  of  data.  The  ordinal-intercept  of  the 
two  sand  strips  plot  was  1.290  and  of  the  three  sand  strips  plot  was  1.356. 

The  difference  in  form  factors  between  the  additional  sand  strips  gave  the  added  form  factor 
due  to  the  sand  strips.  The  difference  between  these  values  was  0.065.  Assuming  that  the  testing 
with  only  one  sand  strip  indicated  the  true  flow  conditions  around  a  full-size  vessel,  the 
difference  from  the  sand  strips  was  subtracted  from  the  stimulated  form  factor.  Therefore,  the 
adjusted  form  factor  for  the  hull  was  0.103.  The  form  factor  of  0.103  was  used  to  scale  up  the 
coefficient  data  to  the  ship  scale. 

An  interesting  look  into  flow  affects  was  seen  by  superimposing  the  viscous  coefficient  over 
the  total  resistance  coefficient  for  each  sand  strip  condition  (Figure  (24)).  The  points  to  notice 
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were  at  the  intersections  of  each  total  resistance  curve  and  viscous  resistance  curve.  Most  of 
these  interactions  occurred  around  Rn  =  375,000.  Data  acquired  at  Reynolds  numbers  lower  than 
375,000  had  significantly  lower  resistance  than  at  points  where  Rn  >  375,000.  The  data  having 
lower  resistance  indicated  a  large  amount  of  laminar  flow,  and  did  not  properly  model  the  flow 
conditions  found  on  the  full-size  vessel. 
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Figure  24  -  Total  and  viscous  resistance  for  all  upright  data 

The  exception  to  the  rule  of  a  decreasing  coefficient  of  total  resistance  with  decreasing 
velocity  is  below  Rn  =  250,000.  These  very-slow  test  points  had  entirely  laminar  flow,  and  the 
laminar  friction  coefficient  dominated  resistance.  Eventually,  these  slow  test  points  would  drive 
the  total  resistance  coefficient  to  a  limit  of  infinity  at  low  velocities  (Eqns  (8),  (13)). 

Since  the  form  factor  was  found,  the  viscous  component  of  resistance  was  subtracted  from 
the  total  resistance  to  leave  wavemaking  resistance  (Figure  (25)).  Notice  that  below  a  Froude 
number  of  0.18  (Rn<400,000),  there  was  a  considerable  amount  of  data  which  had  a  negative 
wavemaking  coefficient.  The  negative  wavemaking  coefficient  indicated  that  there  was  a 
considerable  amount  of  laminar  flow  across  the  model  at  those  points. 
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Figure  25  -  Wavemaking  resistance  in  the  upright  condition 

Froude  proposed  that  ship  predictions  could  be  made  by  scaling  only  the  wavemaking 
resistance  from  model  tests.  The  viscous  component  of  resistance  for  the  ship  would  be 
calculated  using  the  form  factor  which  was  found  for  the  model. 

Powering  required  to  pull  a  ship  through  the  water  is  known  as  effective  horsepower.  Power 
is  the  product  of  force  over  a  velocity,  therefore,  for  English  units  power  is: 

where  veloeity  is  expressed  in  knots,  and  resistanee  is  in  pounds. 

From  the  aequired  upright  data,  both  the  resistanee  and  effeetive  horsepower  were  ealeulated 
for  the  full-size  Mk  II  Navy  44  STC  (Figure  (26)). 
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Figure  25  -  Wavemaking  resistance  in  the  upright  condition 
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Figure  26  -  Upright  resistance  and  horsepower 


Difficulties  in  Tank  Testing 

From  the  beginning  of  tank  testing,  there  were  two  major  problems.  The  first  problem 
consisted  in  realignment.  After  the  initial  alignment,  the  rig  was  assumed  straight.  However, 
after  checking  the  alignment  again  after  some  of  the  initial  unstimulated  flow  cases,  it  was  found 
that  the  square  beam  was  not  rigidly  attached  to  the  carriage.  To  solve  this,  both  the  square  beam 
and  the  "L"  beam  were  screwed  into  the  carriage  -  not  just  clamped. 

Another  major  difficulty  throughout  the  project  was  the  presence  of  transient  noises  in  the 
system  (Figure  (27)).  The  result  of  this  noise  was  to  multiply  the  raw  forces  to  many  times  the 
actual  forces.  Many  Fourier-transforms  were  performed  on  the  data  to  try  to  identify  the  source 
of  the  errors.  The  Fourier-transforms  could  not  find  a  frequency  of  the  noise  which  was 
distinguishable  from  the  standard  data  frequencies. 

Ultimately,  it  was  shown  from  an  oscilloscope  that  the  noise  had  a  beat  signal.  The  source 
of  the  beat  could  not  be  found.  However,  the  data  acquisition  system  was  modified  to  truncate 
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the  data  near  the  zero-erossings  of  the  beats.  This  method  adjusted  the  starting  and  ending  points 
of  the  data  to  find  the  lowest  standard  deviation.  In  between  the  truneation  limits,  the  data  was 
averaged  to  find  the  steady-state  foree  without  the  influenee  of  noise. 


Another  work-around  to  reduee  the  error  assoeiated  with  the  transient  noise  was  to 
drastically  increase  the  rate  of  sampling.  Using  this  method,  a  greater  percentage  of  the  noise 
averaged  itself  out  as  sampling  was  increased  from  25  Hz  to  143  Hz.  This  method  required 
moving  to  a  laptop  instead  of  a  PC.  This  was  because  the  PC's  data  acquisition  was  connected 
by  a  USB  connection  which  could  not  handle  the  high  sampling  rate.  The  laptop  used  a  parallel 
connection  for  data  acquisition  which  could  receive  high  data  rates.  On  the  higher  data  sampling 
system,  the  beat  signal  was  visible  in  the  data  (Figure  (27)). 
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FKS 

FKS  is  a  CFD  code  developed  primarily  from  the  work  of  Dr.  Noblesse  from  NSWC-CD. 
FKS  is  called  a  “free-surfaee  code”  which  describes  its  goal  of  predieting  wave  making 
resistance.  The  free-surface  is  the  boundary  between  two  fluids.  For  ships,  this  boundary  is 
between  the  seawater  and  the  air  and  mathematically  is  modeled  as  a  continuous  surface.  In  real 
life  at  the  free-surface,  non-eontinuum  effects  can  take  place.  For  example,  a  bow  wave  from  a 
ship  ean  often  roll  over  on  itself,  or  spray  can  form  and  break  off  free  from  the  rest  of  the 
surface.  Furthermore,  free-surfaee  eodes  are  diffieult  to  eonstruct  because  little  is  known  about 
the  viscous  relationship  of  wavemaking.  FKS  is  an  inviseid  code,  and  this  viscous  interference  is 
ignored,  leading  to  some  error  in  the  solution.  Determining  this  error  was  a  major  foeus  of  the 
testing. 


FKS  Theory 

As  waves  are  formed  near  the  surface  of  the  hull,  they  radiate  towards  infinity.  However,  as 
they  radiate,  the  waves  interfere  with  each  other.  As  an  example  of  this  interferenee  at  Fn=0.40, 
the  bow  wave  has  construetive  interference  with  the  stern  wave  which  creates  larger  waves. 
Sinee  waves  are  created  from  the  energy  of  the  boat  driving  through  the  water,  it  requires  more 
force  for  the  boat  to  create  the  larger  bow  and  stem  waves.  Therefore,  total  resistanee  inereases 
at  Froude  number  0.40.  For  ships,  Fn=0.40  is  eonsidered  the  “hull  speed.”  Without  extra 
powering,  most  ships  will  not  be  able  to  travel  at  faster  Froude  numbers. 

Farfield  waves  are  the  summation  of  waves  at  a  great  distanee  away  from  the  ship. 
Individual  waves  created  by  the  hull  merge  a  short  distance  away  from  the  ship,  ereating  a 
uniform  wave  pattern  that  propagates  outward.  In  the  absence  of  viseosity,  these  would  travel  to 
infinity.  By  inspecting  the  farfield  waves,  the  total  resistanee  from  all  of  the  wave  interferenee 
ean  be  ealculated.'^'*^^ 


The  Haveloek  formula  (Eqn  (18))  ealeulates  the  resistance  eaused  by  waves  through  a 
calculation  of  the  farfield  wave-speetmm  (S).  In  Eqn  (18),  the  wave-speetram  is  broken  into  its 
real  (Sr)  and  imaginary  parts  (Si). 


71  ^  K-V 


(18) 
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The  Havelock  formula  is  a  function  of  Froude  number,  wavenumber  (k),  the  Fourier  variable  (P), 
and  the  wave-spectrum.  Froude  number  is  incorporated  into  the  equation  via  Eqn  (19)  which 
acts  as  a  description  for  velocity. 


Wavenumber  is  a  description  of  the  wavelength  (X)  of  the  radiated  waves  (Eqn  (20))  and 
therefore  acts  as  a  function  of  wave  celerity  (c),  or  the  speed  of  the  waves  (Eqn  (21)). 


In 

K  = - 

X 


(20) 


Wave  celerity  (Eqn  (21))  is  also  a  function  of  the  frequency  of  the  waves  (co).  The  wave 


frequency,  along  with  wavelength  can  be  used  to  calculate  the  energy  in  each  wave. 


CO 

c  =  — 

K 


Moreover,  wavenumber  is  a  function  of  both  velocity  and  the  Fourier  variable  in  Eqn  (22). ' 


:(^)  —  V  +  -V P‘ 


The  farfield  wave-spectrum  can  be  calculated  as  a  function  of  disturbance  velocity,  or  the 
velocity  of  water  particles  due  to  the  ship’s  velocity,  along  the  hull.  This  calculation,  the 
Fourier-Kochin  representation  of  farfield  waves,  was  used  by  Dr.  Noblesse  to  develop  FKS.'^"^ 


FKS  calculates  the  disturbance  velocity  at  the  hull  using  a  slender-ship  approximation.  This 
approach  treats  a  hull  as  a  set  of  cylinders  without  end  conditions.  The  benefit  of  this  approach 
is  that  calculations  can  occur  on  limited  resources.  The  negative  aspect  of  this  method  is  that  the 
approximation  breaks  down  once  the  ship’s  beam  is  too  large  for  its  length.  Slender-ship 
theory’s  assumptions  will  fail  when  the  beam  is  much  greater  than  the  total  wavelength. 


FKS  calculates  the  velocity  at  the  hull  of  the  ship  via  the  slender-ship  approximation.  Using 
these  velocities,  the  farfield  waves  can  be  calculated  using  the  Fourier-Kochin  method.  Finally, 
the  Havelock  formula  calculates  the  wave  resistance  from  the  farfield  waves. 
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FKS  Process 

In  order  to  calculate  the  disturbance  velocity  at  the  hull,  FKS  needed  to  be  able  to  make 
calculations  along  a  numerical  representation  of  the  hull.  FKS  used  a  three-dimensional  surface 
model  divided  into  panels  to  form  a  structured  grid,  or  mesh  (Figure  (28)).  The  model  itself  had 
to  be  scaled  to  a  non-dimensional  length  of  one  unit. 


Figure  28  -  Nondimensional  Mk  II  Navy  44  STC  with  trapezoidal  mesh 

Care  had  to  be  included  in  the  formation  of  the  meshed  hull.  ^  FKS  used  only  the 
underwater-body  of  the  hull  as  the  input  mesh.  Therefore,  the  mesh  was  constructed  of  the  body 
below  the  waterplane,  or  z=0.  Since  the  initial  free-surface  at  z=0  moves  to  both  above  (+z)  and 
below  (-z)  the  meshed  waterplane,  FKS  calculates  the  hull  at  z>0  as  a  linear  continuation  of  the 
hull  from  below  the  surface  of  the  water.  If  the  meshed  hull  had  protrusions  at  the  waterplane 
which  were  not  characteristic  of  the  continued  hull,  the  calculated  free-surface  would  be 
inaccurate.  An  example  would  be  if  the  boat  was  at  a  maximum  heel  where  the  deck  touched  the 
free-surface.  At  this  point,  FKS  would  perform  its  calculations  as  if  there  was  still  more  hull  for 
the  free-surface  to  attach. 


®  Before  the  beginning  of  this  Trident  Project,  the  author  went  to  “Gridgen  training”  by  the  Gridgen  Corporation  in 
Ft  Worth,  TX.  This  training  taught  how  to  form  usable  meshes  using  the  Gridgen  program  which  were  used 
extensively  in  this  project. 
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The  CFD  meshing  program,  “Gridgen”  was  used 
to  create  a  grid  of  panels  over  the  underwater  surface 
of  the  hull. 

The  first  step  in  using  Gridgen  was  to  create 
items  called  “connectors.”  These  connectors  were 
placed  over  the  boundaries  of  the  hull.  Connectors 
act  as  limits  to  the  area  being  analyzed.  They  also 
determine  the  sizes  and  amounts  of  panels  used  for 

analysis.  Sizes  of  panels  were  determined  through 

Figure  29  -  Sample  x-y  surface  in  Gridgen  .  ,  ,  „  , 

putting  nodes  on  the  connectors  m  a  process  called 

“dimensionalizing.” 

Finally,  an  area  called  a  “domain”  was 
constructed  by  arranging  adjacent  connectors  into  a 
four-sided  surface.  The  order  in  which  the  connectors 
were  picked  was  important  since  the  normal  of  a 
domain  determines  if  the  domain  is  an  inside  or 
outside  piece.  A  domain  with  an  inverted  normal 
would  tell  FKS  that  water  was  flowing  through  the 
inside  of  the  hull,  not  over  the  outside. 

Figure  30  -Dimensionalized  connectors 

After  exporting  the  domains  which  represent  the 
hull  from  Gridgen,  a  preprocessor  was  used  to  cut  the 
trapezoidal  panels  in  half  to  create  triangular  panels. 
The  preprocessor,  included  in  FKS,  called  “g2FKS” 
created  the  final  model  ready  for  computations. 

Besides  the  hull  itself,  the  only  other  input  to 
FKS  was  the  Froude  number.  By  running  a  set  of 
Froude  numbers,  an  upright  resistance  curve  could  be 
completed. 

Fignre  31  -  Gridgen  with  domain  over 


surface 
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FKS  Tests 

To  validate  FKS,  a  sample  hull  known  as  the  Wigley  wedge-hull  was  used.  The  Wigley 
wedge-hull  is  a  mathematical  form  which  has  wavemaking  data  published  in  Principles  of  Naval 
Architecture. 


The  Wigley  wedge-hull  was  created  in  Gridgen,  and  a  range  of  speeds  were  tested.  The 
result  from  FKS  was  the  variation  of  the  wavemaking  coefficient  as  a  function  of  Froude 
number.  Table  (1)  indicates  that  FKS  performed  well  compared  to  the  published  calculated 
coefficient  of  wavemaking. 

_ Values  of  Fn _ 

Wigley  Minima  Cw  0.23  —  0.35  — 

Calculations  Maxima  Cw  — _ 0.27 _ _ 0.48 

FKS  Minima  C*  0.25  —  0.33  — 

Calculations  Maxima  C*  —  0.29  —  0.48 

Table  1  -  FKS  comparison  of  the  Wigley  wedge-hull  to  published  data'*'*' 

The  published  data  was  of  a  simple  closed-form  calculation  of  the  coefficient  of 
wavemaking.  Like  FKS,  the  published  solution  was  also  inviscid.  As  a  consideration  of  the 
published  data,  the  Wigley  calculations  were  from  the  wave  patterns  from  only  four  points  on  the 
hull.  However,  at  Fn=0.48,  the  bow  wave  and  the  stem  wave  dominated  the  resistance  curve. 
That  FKS  exactly  matched  the  closed-form  wavemaking  solution  at  the  higher  Froude  number 
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showed  that  FKS  can  properly  evaluate  a  simple  hull.  It  was  considered  that  since  FKS  was 
more  complex  and  calculated  wave  patterns  from  all  parts  of  the  hull,  the  lower  Froude  numbers 
were  more  accurate  in  the  FKS  model  than  the  published  solution.  At  these  lower  Froude 
numbers,  the  resistance  curve  was  not  dominated  by  just  one  point,  but  considered  the 
wavemaking  of  the  entire  surface. 

For  the  Mk  II  Navy  44,  a  script^  was  written  in  the  TCL  language  to  facilitate  generating 
grids  in  Gridgen.  This  script  received  as  input  any  heel,  yaw,  rudder,  or  trim  angle  on  the  Navy 
44.  Furthermore,  heave  and  the  number  of  panels  on  the  hull  could  be  set  internally  in  the  script. 
Two  subsequent  scripts  were  written  in  the  BASH  shell  which  would  automatically  run  FKS 
over  a  set  range  of  speeds. 

Two  sets  of  testing  were  performed  on  the  Mk  II  Navy  44:  upright  and  sailing  conditions. 
For  the  upright  condition,  all  attitude  angles  were  set  at  zero  degrees.  The  range  of  speeds  was 
set  from  Fn  =  0.05  to  0.70.  Each  condition  took  approximately  20  seconds  to  analyze  on  a 
Pentium  4  laptop  running  Linux.  Sixty-six  conditions  were  analyzed. 

Sailing  conditions  were  composed  of  a  square  test  matrix  which  varied  speed,  heel  angle, 
yaw  angle,  and  rudder  angle.  Speed  ranged  between  Fn  =  0. 1  and  0.5.  Heel  angle  varied  between 
0  and  27.5  degrees.  Yaw  angle  and  Rudder  angle  was  set  as  if  in  the  tow  tank  where  yaw  was  set 
at  0  and  4  degrees.  Rudder  angle  was  similarly  set  at  0,  3,  and  6  degrees.  Five  hundred  forty 
sailing  conditions  were  analyzed. 

FKS  Upright  Analysis 

To  post-process  FKS,  several  corrections  had  to  be  applied.  The  main  problem  with  the  raw 
FKS  data  was  that  the  characteristic  length  used  in  FKS’s  Froude  number  calculations  was  not 
representative  of  the  true  length  of  the  fluid  in  question. 

A  variable  called  “virtual  length”  was  created  to  revise  the  Froude  number  scaling.  Virtual 
length  was  created  to  help  visualize  the  speed  of  the  boat  relative  to  the  virtual  tow  tank  in  FKS’s 
code. 


’  See  Appendix  B  for  the  TCL/CL  Code 
*  See  Appendix  B  for  the  BASH  Codes 
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^fluid  body 


V'hull 


Figure  33  -  FKS  corrections 


The  first  step  in  correcting  the  Fronde  number  was  to  use  the  entire  model  length  as  the 
fluid  body  (Figure  (33).  This  allowed  the  body  of  water  to  travel  as  fast  as  the  entire  hull  (Lhuii). 
The  second  correction  was  to  divide  the  nondimensional  length  (Lwifks)  by  the  length  of  body 
used  in  FKS  (Lfks)-  However,  the  second  correction  is  not  reducible  to  a  value  of  one.  One 
more  correction  was  needed  to  reduce  the  terms  to  one.  The  nondimensional  length  was  a 
waterline  length,  while  the  term  Lfks  was  a  total  length.  Therefore,  the  preceding  term  had  to  be 
multiplied  by  the  value  of  the  total  length  (Ltot)  over  the  waterline  length  (Lwi). 


(23) 


The  upright  resistance  curve  could  now  be  evaluated  and  compared  to  the  tank  testing  curve. 
The  total  resistance  curve  (Figure  (34))  for  FKS  included  the  wavemaking  coefficient  and  a 
viscous  coefficient  based  off  the  ITTC  viscous  equation  (Eqns  (8),  (9)).  In  the  viscous 
formulation,  the  form  factor  used  was  k=0.103. 


FKS  did  not  use  the  typical  hydrodynamic  wavemaking  coefficients  as  in  Eqn  (4).  Instead, 
FKS  uses  the  Eqn  (24)  where  wavemaking  resistance  is  defined  as  (Rw)  and  the  velocity  of  the 
fluid  is(U).["^ 


(24) 
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Figure  34  -  Upright  resistance  from  FKS  and  tank  testing 

FKS  compared  favorably  to  a  complex  hull  form  like  the  Mk  II  Navy  44.  However,  as 
Froude  number  increased  after  Fn=0.35,  spray  would  develop  as  the  assumptions  from 
continuum  mechanics  break  down  in  that  the  free-surface  would  roll  over  on  itself.  This  effect 
became  prominent  after  Fn=0.5.  In  Figure  (34),  FKS  helped  to  show  where  laminar  flow  existed 
on  the  tank  model  in  areas  mostly  below  Fn  =  0.2. 

Furthermore,  the  local  velocities  at  the  hull  calculated  using  slender-ship  theory  did  not 
account  for  the  interaction  of  the  waves  generated  upstream  (near  the  bow).  While  slender-ship 
theory  is  a  very  quick  method  of  determining  the  local  velocities,  its  approximations  can  break 
down  for  a  beamy  hull  like  the  Mk  II  Navy  44.  At  high  speeds,  the  hull’s  wider  entry  angle 
would  cause  local  velocity  disturbances  which  are  not  accounted  for  in  slender  ship  theory.  The 
presence  of  this  effect  contributed  to  the  error  seen  at  the  high  Froude  numbers  (Fn  >  0.5). 

Finally,  the  humps  and  hollows  of  the  coefficient  of  total  resistance  curve  also  existed  in  the 
coefficient  of  wavemaking  curve.  Because  FKS  performed  a  pure  inviscid  wavemaking 
calculation,  the  output  was  thought  of  as  a  numerical  approximation  to  a  closed  form  inviscid 
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wavemaking  solution.  In  this  solution,  the  humps  and  hollows  would  be  magnified  since 
viscosity  would  normally  decrease  wavemaking  effects. 

Further  explanations  of  FKS  results  are  located  in  the  conclusions. 
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SPLASH 

SPLASH  is  a  free-surface  CFD  code  which  was  developed  especially  for  sailboats.  It  has 
been  used  extensively  in  the  America’s  Cup  and  is  considered  one  of  the  leading  CFD  codes  for 
sail  craft.  SPLASH  is  a  panel-based  code,  meaning  that  it  calculates  potential  flow,  and  much 
revision  has  gone  into  the  code  to  evaluate  the  effectiveness  of  lifting  surfaces. 

SPLASH  Theory 

SPLASH  stands  for  Small  Perturbation  Linearized  Analysis  of  free-Surface 
Hydrodynamics.  SPLASH  is  inviscid,  meaning  that  shear  stresses  from  the  fluid  are  not 
calculated.  Secondly,  SPLASH  calculates  the  flow  of  incompressible  fluids.  Incompressibility 
combined  with  inviscidness  creates  the  conditions  for  an  ideal  fluid  meaning  that  the  potential 
flow  of  the  fluid  can  more  easily  be  calculated. 

Potential  flow  can  be  calculated  from 
singularities  called  sources  and  sinks.  To 
visualize  this  concept,  sources  can  be 
thought  of  as  points  which  radiate  fluid 
symmetrically,  while  sinks  suck  fluid 
(Figure  (35)).  In  reality,  sources  and  sinks 
simulate  the  flow  around  a  body  by  acting 
as  either  receptors  or  transmitters  of 
velocity  potential  (O).  The  underlying 
property  of  potential  flow  is  found  in  Laplace’s  equation  for  valid  fluid  flow  where  the  velocity 
potential  (O)  must  reach  equilibrium: 


Figure  35  -  Potential  flow'*®' 


=  0  (25) 

From  the  velocity  potential  (O),  the  nondimensional  velocity  vector  (V)  can  be  calculated: 

[12] 


F  =  VO  =  (!),/  + 0^7 +  0  A 


(26) 
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SPLASH  Testing 

SPLASH  was  compiled  on  an  IRIX  SGI  machine  with  eight  processors.  Because  of  the 
many  differing  versions  of  POSIX  based  systems  and  the  complexity  of  the  Mk  II  Navy  44,  the 
setup  of  SPLASH  exclusively  took  almost  two  months. 

As  in  FKS,  the  first  step  to  running  the  code  was  to  generate  a  mesh  around  the  body.  The 
current  version  of  SPLASH  was  ideally  suited  to  handle  International  America’s  Cup  Class 
(lACC)  Yachts.  SPLASH  normally  automatically  creates  the  mesh  required  to  set  up  a  hull  for 
testing.  However,  the  Mk  II  Navy  44  differed  quite  a  bit  from  an  lACC  boat.  Changes  first  had 
to  be  made  to  the  Fortran  setup  files  in  SPLASH  to  eliminate  the  calculations  for  a  keel-bulb  and 
winglets. 

Elements  known  as  “iges  surfaces”  defined  the  hull  shape  for  the  Mk  II  Navy  44  three- 
dimensionally.  SPLASH  normally  had  the  ability  to  determine  how  to  place  meshes  onto  the 
iges  surfaces  of  the  hull.  However,  the  iges  surface  of  the  Mk  II  Navy  44  hull  was  of  an  irregular 
format,  and  automatic  meshing  failed.  Instead,  Gridgen  was  again  used  to  generate  an  initial 
mesh  for  SPLASH. 

SPLASH  was  made  more  accurate  by  incorporating  sail  forces  into  the  setup  files.  SPLASH 
is  a  valuable  tool  for  sailing  vessels  because  in  its  internal  calculations,  SPLASH  will  change  the 
attitude  of  the  boat  based  on  an  equilibrium  found  with  hydrodynamic  and  aerodynamic  forces. 
Sail  forces  could  be  determined  either  from  a  VPP  or  from  tank  testing  data.  In  the  tank, 
aerodynamic  forces  are  equal  and  opposite  of  the  hydrodynamic  forces.  Therefore,  the  tank  data 
of  the  Mk  II  Navy  44  was  modified  to  become  aerodynamic  data  and  was  inserted  into 
splash’s  sail  forces  file. 

Finally,  the  actual  conditions  to  be  analyzed  in  SPLASH  had  to  be  coded  into  FORTRAN. 
Conditions  were  chosen  based  on  a  method  that  SPLASH  needed  to  post-process  the  data.  For 
the  upright  condition,  40  points  were  picked.  For  the  sailing  condition,  heel  angle  was  set  at  0, 
10,  15,  20,  and  25  degrees  while  speed  ranged  from  0  to  10.5  knots.  Yaw  and  rudder  angle 
varied  as  in  FKS  and  tank  testing  -  using  six  tests  for  each  heel  and  speed  combination. 

SPLASH  Upright  Analysis 

As  in  both  tank  testing  and  FKS,  the  output  from  SPLASH  consisted  of  upright  and  sailing 
data.  SPLASH  did  not  require  any  scaling  of  the  output  since  forces  were  computed  at  the  full- 
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scale.  Furthermore,  it  was  not  necessary  to  calculate  the  viscous  drag  of  the  ship  in  the  post¬ 
processing.  This  was  performed  as  a  part  of  SPLASH’s  computations. 

Figure  (36)  shows  how  SPLASH  compared  the  upright  resistance  of  the  Mk  II  Navy  44  to 
tank  testing  and  FKS  predictions.  At  high  velocities,  SPLASH  predicted  slightly  greater 
resistance  for  the  ship  compared  to  both  tank  testing  and  FKS.  One  explanation  for  this  is  in  the 
viscous  drag  calculation,  called  “viscous-stripping.”  In  tank  testing  and  FKS,  viscous  drag  was 
computed  using  a  standard  wetted  surface  area.  However,  since  SPLASH  calculated  the  changed 
wetted  surface  area  for  each  sailing  condition,  the  updated  area  was  used  in  the  calculation.  The 
result  of  this  is  that  SPLASH  more  accurately  predicted  viscous  drag  than  the  other  two  methods. 


Figure  36  -  Resistance  of  all  upright  data 

Interestingly,  at  a  speed  of  10.5  knots,  SPLASH  predicted  a  significant  drop  in  resistance. 
At  this  speed,  SPLASH  also  predicted  that  the  pitch  angle  of  the  boat  would  jump  as  well.  As 
boat  speed  would  accelerate  through  10.5  knots,  the  vessel  would  actually  begin  to  ride  over  the 
bow  wave  (a  similar  principle  to  breaking  the  sound  barrier  in  aerodynamics).  Moving  over  the 
bow  wave  would  create  trim  by  the  stem  of  about  5  degrees.  The  effect  of  this  is  that  the  bow 
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would  rise  out  of  the  water,  decreasing  the  waterline.  At  this  velocity,  the  bow  wave  would 
begin  to  disintegrate,  leaving  only  a  partial  bow  wave  and  a  stem  wave.  Without  the  massive 
bow  wave,  the  coefficient  of  wavemaking  would  drop.  Furthermore,  with  a  decreased  wetted 
surface  area,  viscous  resistance  would  decrease  as  well.  This  speed,  known  as  semi-planing  is 
characterized  by  a  decrease  in  the  coefficient  of  total  resistance  and  facilitates  high-powered 
boats  to  reach  speeds  much  larger  than  their  own  hull  speeds.  That  the  tank  test  data  did  not 
clearly  indicate  such  a  condition  of  semi-planing,  was  not  cause  for  questioning  SPLASH ’s 
predictions.  On  the  contrary,  there  were  not  enough  upright  high-speed  tests  performed  in  tank 
testing  to  sufficiently  populate  this  region. 


Figure  37  -  Total  resistance  coefficient  of  all  npright  data 

Noticeably,  SPLASH  did  not  have  the  minor  humps  and  hollows  in  the  coefficient  curves 
that  were  present  from  FKS  (Figure  (37)).  Since  SPLASH  calculated  only  the  potential  flow, 
some  of  the  wave  interference  effects  were  minimized  when  compared  to  FKS.  However,  many 
of  the  humps  and  hollows  found  in  FKS  would  not  be  as  great  if  viscous  interference  were  to 
occur.  Therefore,  the  actual  coefficient  of  wavemaking  curve  was  estimated  to  be  between  the 
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curves  of  SPLASH  and  FKS.  The  use  of  a  viscous  CFD  code  could  help  to  identify  where 
discrepancies  exist. 

Furthermore,  it  was  noticed  that  at  0.1  <  Fn  <  0.2  there  was  a  difference  in  the  coefficient  of 
total  resistance  between  FKS  and  SPLASH.  At  this  speed,  viscous  resistance  dominated  the 
resistance  on  the  boat  and  this  again  indicated  that  the  viscous  resistance  calculation  differed 
between  both  methods. 
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Sailing  Data  Analysis 

Tank  testing  provided  120  data  points  for  sailing  conditions.  Each  point  varied  by  boat 
speed,  heel  angle,  rudder  angle,  and  yaw  angle.  FKS  also  gave  force  outputs  for  a  much  larger 
matrix  of  over  500  data  points.  SPLASH  generated  a  matrix  of  150  data  points.  The  difficult 
part  of  analyzing  the  sailing  conditions  was  that  each  method  generated  results  for  slightly 
different  sailing  conditions. 

Many  attempts  were  made  at  making  use  of  the  sailing  data  before  a  sure  method  was 
created.  The  main  problem  with  the  sailing  data  was  that  the  sailing  conditions  tested  in  the 
tank,  FKS,  and  SPLASH,  did  not  exactly  match  the  predicated  conditions  from  the  PCSail  VPP. 
In  the  tank,  this  was  because  the  velocity  input  for  the  tank  could  only  be  set  approximately  to 
the  desired  velocity.  Furthermore,  the  internal  data  analysis  in  SPLASH  confined  tests  from 
being  singular  test  conditions.  In  SPLASH,  the  PCSail  VPP  predicted  conditions  could  not  be 
individually  tested  efficiently. 

The  first  attempt  of  calculations  with  the  tank  data  was  to  perform  a  regression  of  all  of  the 
data  points  to  find  functions  for  drag,  lift,  and  the  aft-lift  terms.  These  functions  would  vary  by 
velocity,  heel  angle,  yaw  angle,  and  rudder  angle.  Several  multivariate  regressions  were  made, 
but  ultimately,  a  polynomial  expression  for  the  data  with  acceptable  accuracy  could  not  be  found. 

A  second  approach  used  MATLAB  to  interpolate  between  the  acquired  test  points.^  The 
interpolation  technique  took  the  sparsely  populated  data  matrix  and  transformed  it  into  a  square 
matrix  using  splines  to  calculate  intermediate  points.  Since  the  transformed  matrix  was  both 
square  and  dense,  linear  interpolation  through  this  matrix  could  be  used  to  find  the  required 
values  for  any  velocity  and  heel  angle.  Each  yaw  angle  and  rudder  angle  condition  created  six 
matrices  each  for  the  drag,  fwd-lift,  and  aft-lift  variables. 

Making  square  matrices  of  the  tank  data  provided  not  only  a  method  to  perform 
interpolation,  but  also  visualization  and  analysis  of  the  acquired  sailing  data.  Figure  (38)  shows 
the  collected  data  laid  onto  the  matrix  used  in  calculations. 


®  See  Appendix  C  for  the  MATLAB  script 
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Figure  38  -  Tank  testing  drag  with  0  degrees  rudder  and  0  degrees  yaw 

Figure  (38)  indicates  that  for  a  given  heel  angle,  drag  did  not  change  dramatically.  Upon 
further  inspection  of  the  data,  two  things  became  apparent.  First,  values  chosen  for  calculations 
were  required  to  be  inside  the  range  of  tested  values.  Furthermore,  between  the  testing  matrix 
and  velocity  =  0,  the  splines  created  a  very  slight  negative  drag  in  order  to  smooth  their  curves. 
At  speeds  and  heels  outside  the  range  of  the  testing  matrix,  there  was  no  data  and  extrapolation 
would  be  inappropriate.  This  applied  for  all  data  collected  (drag,  fwd-lift,  and  aft-lift).  Notice 
that  the  square  matrix  in  Figure  (39)  modified  for  calculations  reported  null  answers  for  any  data 
which  was  not  in  the  range  of  interpolation. 

Secondly,  drag  did  not  change  dramatically  between  heel  angles,  but  Figure  (39)  shows  that 
there  was  a  significant  amount  of  change  in  drag  across  a  range  of  heel  angles  with  constant 
model  speed.  This  shows  where  tank  testing  and  CFD  have  great  resolution.  The  change  in  drag 
across  a  range  of  heels  is  a  function  of  only  the  shape  of  any  particular  hull.  A  VPP  would  not 
be  able  to  detect  these  fine  changes  in  drag. 
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4.5  r 


Figure  39  -  Closer  inspection  of  tank  testing  drag  with  0  degrees  rudder  and  0  degrees  yaw 

For  tank  testing,  CFD,  and  FKS,  six  drag  matrices  were  made  for  all  rudder  angles  and  yaw 
angles.  One  of  the  disadvantages  of  the  method  of  making  square  matrices  for  the  sailing 
conditions  had  to  deal  again  with  resolution.  Individual  points,  such  as  the  point  at  model  speed 
4.2  ^s  and  heel  of  10  degrees  (Figure  (38)  had  a  great  effect  on  the  shape  of  the  entire  matrix. 
Notice  that  above  4.2  Q)s,  around  10  degrees  of  heel,  a  dip  was  present  (Figure  (39).  The 
precision  of  the  test  at  this  point  could  be  drawn  into  question,  especially  since  it  had  far 
reaching  affects  on  the  rest  of  the  matrix. 
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Figure  40  -  FKS  drag  with  0  degrees  rudder  and  0  degrees  yaw 

FKS  did  not  show  a  dip  in  drag  at  10  degrees  of  heel  at  all  (Figure  (40)).  However,  FKS 
showed  interesting  humps  right  underneath  0.5  pounds  of  drag.  These,  unfortunately  were  in 
error,  and  the  sailing  data  for  FKS  will  later  be  shown  to  be  unrealistic. 
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3.0 


Heel  angle  [deg] 


Figure  41  -  SPLASH  drag  with  0  degrees  rudder  and  0  degrees  yaw 

Since  SPLASH  results  were  given  in  the  full-size  condition,  the  results  presented  in  Figure 
(41)  were  calculated  to  model  scale  to  facilitate  the  comparison  between  the  three  experimental 
methods.  First,  total  model  drag  in  Figure  (41)  did  not  have  as  high  a  limit  in  SPLASH  as  in 
tank  testing  or  FKS.  The  speed  of  the  tests  in  SPLASH  were  limited  to  5  Q)s,  unlike  tank  testing 
which  proceeded  to  5.6  Q)s  for  the  sailing  conditions. 

Secondly,  SPLASH  did  not  show  as  much  variation  in  drag  as  a  function  of  heel  angle  as 
when  compared  to  tank  testing.  This  further  draws  into  question  the  precision'^  of  individual 
points  in  tank  testing.  A  quality  of  CFD  which  was  exposed  is  that  its  precision,  unlike  tank 

In  engineering  terminology,  “precision”  is  defined  as  repeatability  of  tests  while  “accuracy”  is  defined  as  absolute 
error  of  the  tests  from  a  datum. 
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testing,  can  always  be  100%.  This  made  matrices  such  as  these  very  able  to  obtain  reasonable 
results.  However,  in  CFD,  if  the  meshing  of  a  model  at  certain  conditions  was  incorrect,  the 
accuracy  of  the  entire  grid  would  be  very  low. 

The  questionable  accuracy  was  the  case  for  FKS.  When  forming  the  meshes  for  the  sailing 
conditions,  FKS  suffered  from  two  items.  First,  as  the  model  was  heeled  in  Gridgen,  it  was  not 
heaved  at  corresponding  amounts  to  control  the  displacement.  This  created  a  mesh  which  was 
not  an  accurate  representation  of  the  hull.  Secondly,  FKS  did  not  iteratively  remesh  the  hull 
during  its  calculations.  This  tended  to  leave  the  hull  in  a  position  which  was  also  not  accurate. 
Finally,  FKS  was  not  programmed  to  make  lift  circulation  calculations.  The  meaning  of  this  is 
that  FKS  could  not  accurately  model  lifting  forces  or  induced  drag  due  to  lifting  surfaces. 
Therefore,  the  use  of  FKS  for  resistance  calculations  was  only  appropriate  for  the  upright 
condition  with  no  heel  angle. 

SPLASH  also  had  accuracy  problems.  First,  although  the  hull  was  remeshed  during  the 
calculations,  the  data  which  provided  the  aerodynamic  force  balance  was  calculated  from  tank 
testing  data.  The  aerodynamic  force  balance  provided  only  a  small  input  into  the  trim  of  the 
vessel,  but  the  realization  is  that  the  accuracy  of  the  remeshing  was  only  as  high  as  the  precision 
and  accuracy  of  the  tank  testing. 

Interesting  plots  were  seen  as  a  series  through  different  rudder  and  yaw  angles.  Figure  (42) 
shows  sailing  data  recorded  from  the  aft-lift  force  block  from  the  tank.  These  forces,  when 
multiplied  by  the  distance  between  the  aft-pin  and  the  heave  post  gave  the  yawing  moment. 
From  this  data,  the  effect  of  yaw  angle,  and  rudder  angle,  and  heel  angle  on  the  yawing  moment 
were  investigated. 

Similarly,  Figure  (43)  is  the  yaw  moment  measured  in  SPLASH.  The  moment  was  around 
the  heave  post  at  its  placement  on  the  model.  Although  Figure  (42)  is  force  data  and  Figure  (43) 
is  moment  data,  both  plots  are  proportional. 


Aft  Uft  (lbs)  Aft  Uft  (lbs) 
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Aft  Uft  with  0  degrees  of  yaw  and  0  degrees  of  rudder 


Speed  [ft/s] 


Aft  Uft  with  4  degrees  of  yaw  and  6  degrees  of  rudder 


Speed  [ft/s] 


Figure  42  -  Data  from  the  aft-lift  force  block  in  tank  testing 


YawMoment  [ft-lbs] 
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YawMoment  with  0  degrees  of  yaw  and  0  degrees  of  rudder 


Angle  [d 


Speed  [ft/s] 

YawMoment  with  4  degrees  of  yaw  and  6  degrees  of  rudder 


Angle  [d 


Speed  [ft/s] 


Figure  43  -  Full-scale  yaw  moment  around  the  heave  post  in  SPLASH 
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Both  tank  test  data  and  SPLASH  data  showed  very  similar  results  with  respect  to  yawing 
moment.  The  similarity  in  both  acted  to  validate  each  other  by  arriving  at  the  same  answer  from 
two  very  different  methods.  The  yawing  moment  was  a  difficult  measurement  to  make  in  the 
tank  because  of  scaling  errors.  Because  the  size  of  the  rudder  was  so  small,  it  had  a  high  chance 
of  laminar  flow  thereby  not  creating  the  necessary  lift  across  the  rudder  or  yaw  moment  for  the 
entire  hull.  The  yawing  moment  was  also  a  difficult  computation  to  make  in  CFD  but  for 
different  reasons.  On  a  sailboat,  there  is  a  strong  interaction  of  flow  between  the  keel,  the  hull, 
and  the  rudder.  If  some  part  of  the  mesh  around  an  appendage,  or  between  an  appendage  and  the 
hull,  was  flawed,  then  the  entire  lift  circulation  calculated  by  the  appendage  would  have 
tremendous  error. 


However,  there  were  differences.  For  instance,  in  the  tank,  especially  at  the  condition  of  0 
degrees  yaw,  heel,  and  rudder,  there  was  a  measurable  yaw  moment  at  high  velocities.  This 
moment  was  attributed  to  large  amounts  of  heave  and  trim  in  the  model  at  high  speeds.  The 
effect  of  this  movement  was  that  the  distance  between  the  aft-pin  and  the  force  block  tried  to 
lengthen  and  therefore  put  tension  on  the  aft-lift  force  block.  The  presence  of  this  effect  is 
visible  in  Figure  (44). 


Figure  44  -  High  speed  0  degrees  yaw,  heel,  rudder  tank  test 
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The  conclusion  is  that  although  both  methods  had  good  abilities  at  measuring  very  specific 
forces  on  the  hull,  SPLASH  was  far  more  able  to  accurately  measure  these  minute  forces  without 
any  induced  testing  error. 

From  a  sailing-analysis  perspective,  both  methods  showed  interesting  results  for  yaw- 
moment  in  that  heel  moment  created  a  tremendous  amount  of  inverse  yaw.  The  meaning  of  this 
is  that  the  more  the  Mk  II  Navy  44  heeled  over,  the  more  it  wanted  to  put  its  bow  into  the  wind. 
Although  this  was  already  known  to  sailors  who  have  experienced  this  effect  of  “rounding  up,” 
the  quantification  of  this  principle  showed  that  it  occurred  most  at  high  boat  speeds.  Notably, 
both  SPLASH  and  tank  testing  predicted  the  same  speed  for  the  maximum  rounding  up  effect:  7 
to  7.5  knots. 

One  of  the  tremendously  useful  abilities  of  SPLASH  was  in  its  visualization  of  each  test. 
Figure  (45)  shows  the  result  of  the  flow  calculations  visually  which  helped  analysis  in  two  ways. 
First,  visualization  helped  ensure  that  there  were  no  meshing  errors  associated  with  the 
calculations.  Secondly,  Figure  (45)  and  other  angles  of  the  image  show  the  contour  plot  of  the 
pressure  distribution  across  the  entire  hull  and  free-surface.  Labeled  Cp  for  coefficient  of 
pressure,  this  pressure  diagram  could  be  integrated  to  calculate  the  center  of  lateral  resistance  for 
each  appendage  or  for  the  entire  hull.  Not  only  did  Cp  show  the  result  of  the  calculations,  but  Cp 
could  aid  with  the  development  of  the  hull  form  or  even  the  design  of  the  sail  arrangements. 
Stream  vectors  along  with  Cp  also  showed  areas  where  a  design  may  not  be  fair.  For  instance,  at 
the  rudder  tip,  large  vectors  existed  indicating  non-continuous  flow  calculations.  Although  some 
of  these  vectors  may  be  the  result  of  induced  drag  vortices,  it  also  indicated  that  the  tip  of  the 
rudder  was  not  a  smooth  surface." 


"  After  analyzing  the  grid  near  the  rudder,  Joe  Laiosa  and  Bruce  Rosen,  the  creators  of  SPLASH,  reported  that  the 
tip  of  the  rudder  created  from  the  iges  surface  itself  was  not  fair 
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MarkllNavy445TC 


Speed  =  7.1  knots 
Heel  =  15  deg 
Yaw  =  4  deg 
Rudder  =  6  deg  ^ 


Figure  45  -  Potential  flow  visualization  in  SPLASH 
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MAPLE  was  used  to  process  data  between  the  yaw  and  rudder  angles.'^  Since  all  variables 
were  considered  to  behave  linearly  in  yaw,  an  equation  for  yaw  angle  was  formed  which  could 
interpolate  between  each  matrix.  Similarly,  there  were  three  rudder  angles  tested  for  each  point. 
A  quadratic  regression  was  used  to  solve  between  different  yaw  angles.  Ultimately,  a  function 
was  created  which  treated  yaw  and  rudder  as  two  axes  of  a  three-dimensional  surface.  The 
surface  (Q  was  defined  in  terms  of  rudder  angle  (y)  and  yaw  angle  (X)  by  the  function: 

r) = *  /  ■  ^  V(4,0)  -  i  r  a  ■  f(0,0) + i  /  a .  4-(o,3) 

-mo)-}/ -cm 

+irC(4,3)+ir-A-cm)-^rA-a4,6)-ir<(0,0) 

+ f  r-  «o,3)  -  i  r  V(0,6)  -  4  ^  ■  ao.o) + i  A  ■  a4,0) + «0,0) 

An  example  of  this  interpolation  between  matriees  is  shown  in  Figure  (46)  where  surfaee  (Q 
is  model  drag  at  2.5  feet  per  seeond  and  10  degrees  heel  with  the  results  from  tank  data. 
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Figure  46  -  Model  drag  from  tank  data  defined  by  rudder  angle  and  yaw  angle 


Now  that  an  interpolation  method  was  found,  the  tank  data  was  compared  to  PCSail.  For 
each  sailing  condition,  an  interpolated  hydrodynamic  tank  condition  was  found.  Each  condition 

See  Appendix  C  for  the  MAPLE  script 
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had  a  specified  heel  angle  and  boat  speed  which  could  be  matched  as  a  tank  condition.  Each 
PCSail  condition  also  included  the  aerodynamic  forces  acting  on  the  sails.  Therefore,  the 
aerodynamic  drive  and  side  force  were  matched  with  the  hydrodynamic  data  found  in  the  tank. 

Unfortunately,  there  were  two  problems  with  this  comparison  method.  First,  the  tank  data 
either  over  or  under-predicted  the  PCSail  forces  by  up  to  20%  for  each  condition.  This 
percentage  of  error  was  unacceptable  in  making  comparisons  between  tank  and  CFD  prediction 
methods.  The  source  of  this  error  was  due  to  the  simplified  algorithms  in  the  hydrodynamic 
code  of  PCSail.  The  predicted  velocity  and  heel  parametrically  determined  from  PCSail  for  a 
sailing  condition  was  not  hydrodynamically  the  same  as  the  similar  condition  measured  in  either 
CFD  or  the  tank. 

Secondly,  yaw  and  rudder  angle  could  not  be  solved  using  an  interpolative  comparison  to 
the  PCSail  data.  Since  the  forces  found  in  the  tank  included  drag,  fwd-lift,  and  aft-lift,  there 
were  too  many  variables  to  solve  for  both  yaw  and  rudder  angles. 

The  solution  was  to  create  a  new,  customized  Velocity  Prediction  Program  which  was 
focused  only  on  the  Mk  II  Navy  44  STC. 
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Navy  44  Velocity  Prediction  Program 

Having  the  hydrodynamic  data  acquired  from  multiple  sources,  there  were  two  components 
from  a  velocity  prediction  program  which  were  missing:  aerodynamic  data  and  a  solver. 
Because  of  its  ability  to  show  internal  results  and  generate  plots  useful  in  evaluating  the  code, 
Excel  was  chosen  as  the  platform  to  handle  the  code.  Acquired  data  from  tank  testing,  FKS,  and 

SPLASH  were  stored  in  a  worksheet,  and  most  of  the  code  was  written  in  Visual  Basic  as  well  as 

1  ^ 

cell  functions  in  the  spreadsheet. 

Aerodynamic  Theory 

The  aerodynamic  data  was  calculated  using  the  Hazen  method,  as  in  PCSail.  As  wind 
passes  over  a  sail  from  a  certain  direction,  lift  and  drag  are  created  as  with  a  wing.  The  lift  is 
perpendicular  to  the  angle  of  attack  of  the  wind,  while  the  drag  is  parallel  to  the  angle  of  attack 
of  the  wind.  Lift  and  drag  are  calculated  in  terms  of  coefficients  (Eqn  (28)). 

Faero=lCp^,AV^  (28) 

There  is  no  difference  between  this  calculation  of  aerodynamic  force  and  of  hydrodynamic 
force.  Designation  changes  are  that  the  density  of  the  fluid  (pair)  is  multiplied  by  the  nominal 
area  of  the  sail  (A)  and  the  velocity  of  the  surrounding  air  (V). 

The  coefficient  of  lift  (Cl)  generated  by  a  sail  is  assumed  a  function  only  of  apparent  wind 
angle.  There  are  different  coefficients  of  lift  for  each  sail:  mainsail,  jib,  and  spinnaker. 

The  coefficient  of  drag  (Cd)  is  the  summation  of  three  variables:  parasitic  drag,  induced 
drag,  and  windage.  Parasitic  drag  (Cdp)  is  the  friction  associated  with  wind  passing  over  a  sail. 
Like  lift,  it  is  assumed  also  to  be  only  a  function  of  apparent  wind  angle. 

Induced  drag  (Cdi)  is  the  result  of  vortices  created  as  the  airflow  is  sucked  around  the 
windward  portion  of  the  sail  to  the  leeward  side  (Eqn  (29)).  Induced  drag  is  a  result  of  the 
pressure  differential  on  the  sail  and  the  end  condition  at  the  top  of  the  sail.  Based  on  simple 
wing  theory,  it  is  a  function  of  the  square  of  the  lift  coefficient  for  a  sail.  Furthermore,  the  aspect 
ratio  (AR)  of  a  sail  plays  an  important  role  in  sail  efficiency.  A  high  aspect  ratio  sail  will  lessen 
the  induced  drag  by  effectively  delaying  the  amount  of  vortices  which  can  be  created. 


See  Appendix  D  for  the  Mk  II  Navy  44  yPP 
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The  aspect  ratio  used  in  calculating  induced  drag  is  not  of  any  particular  sail,  but  of  the 
entire  sailing  vessel.  For  upwind  sailing,  the  height  of  the  deck  above  the  water  (FBD)  is 
assumed  to  play  an  effect  in  the  induced  drag  of  the  vessel.  Therefore,  there  are  two  equations 
for  the  calculation  of  aspect  ratio.  Close-hauled,  aspect  ratio  is  a  function  of  the  freeboard 
(FBD),  the  effective  height  of  the  mast,  and  the  nominal  sail  area  (An). 

_{l.\(EHM  +  FBD)y 


For  any  other  sailing  condition  besides  upwind  (true  wind  angle  >  45  degrees),  the  freeboard 
term  is  simply  removed  from  the  aspect  ratio  equation. 

Windage  (Cdo)  is  a  function  of  the  characteristics  of  the  vessel,  and  is  a  crude  method  of 
determining  the  aerodynamic  drag  of  the  rigging  and  hull  (Eqn  (31)). 


C^^=1.13 


(b^„fbd)+(ehmmd) 


Windage  is  determined  from  the  maximum  beam  of  a  vessel  (Bmax),  the  average  freeboard, 
the  effective  height  of  the  mast  (EHM),  the  diameter  of  the  mast  (MD),  and  the  nominal  area  of 
the  sails  (An).  For  the  Mk  II  Navy  44,  all  of  these  variables  were  listed  in  the  IMS  certificate 
from  David  Pedrick. 


Compiled  from  experimental  sources,  the  IMS  recently  updated  and  published  its 
coefficients  of  lift  and  parasitic  drag  which  are  used  in  their  VPP.  Each  sail  has  its  own 
coefficients.  However,  interaction  between  the  headsail  and  the  mainsail  is  neglected.  For  the 
mainsail,  curve  fit  functions  were  created  to  turn  the  experimental  data  into  a  usable  form  for  the 
VPP  code, 
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Figure  47  -  Mainsail  lift  and  parasitic  drag  coefficients 

Figure  (47)  shows  how  both  lift  and  parasitic  drag  vary  as  a  function  of  apparent  wind  angle. 
Interestingly,  high  values  of  drag  at  wind  angles  aft  of  90  degrees  actually  increase  the  speed  of 
the  vessel.  At  these  wind  angles,  the  sails  act  more  to  catch  wind  rather  than  to  act  as  a  wing. 

The  coefficients  of  lift  and  parasitic  drag  are  presented  in  Figure  (48).  The  values  of  lift  and 
drag  for  the  jib  end  at  an  apparent  wind  angle  of  100  degrees.  After  this  angle,  the  use  of  the  jib 
would  be  of  limited  efficacy.  Assuming  that  the  rig  could  not  be  tom  down  for  any  wind  speed, 
the  use  of  a  spinnaker  at  these  high  apparent  wind  angles  would  always  produce  faster  boat 
speeds.  If  a  jib  was  to  be  used  at  high  wind  angles,  it  would  be  “poled-ouf  ’  and  used  on  the 
opposite  side  of  the  boat.  Experimental  data  of  this  effect  was  not  used.  Furthermore,  a 
spinnaker  would  still  be  more  efficient. 
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Figure  48  -  Headsail  lift  and  parasitic  drag  coefficients 


The  aerodynamic  lift  and  total 
drag  coefficients  are  recalculated  to 
a  new  axis.  Since  these  coefficients 
are  still  in  the  apparent  wind  axis, 
lift  and  total  drag  are  recalculated  as 
the  variables  drive  and  heel  force. 
Furthermore,  drag  and  heel  force  are 
calculated  not  being  towards  the 
bow,  but  in  the  direction  of  motion 
and  in  the  plane  of  heel  (O),  not  the 
waterplane  with  the  inclusion  of  yaw 
angle  (Figure  (49)). 


Figure  49  -  Sail  forces  in  the  heeled  xy-plane 
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The  side  force  on  the  boat  is  the  lateral  force  on  the  boat  in  the  xy-plane.  This  is  calculated 
as  the  heel  force  multiplied  by  the  cosine  of  the  heel  angle. 

To  calculate  the  total  aerodynamic  forces,  a  few  more  corrections  were  made.  First,  based 
on  the  Hazen  model  of  aerodynamics,  the  sails  could  be  trimmed  to  prevent  overpowering  the 
boat.  Overpowering  occurs  when  too  much  sail  area  is  exposed,  or  when  the  exposed  sail  area 
creates  too  much  lift.  The  typical  effect  of  overpowering  is  that  the  vessel  heels  over  too  far.  If 
heel  becomes  too  great,  deck-edge  immersion  occurs  causing  significant  drag  increases  and  a 
reduced  maximum  velocity.  Furthermore,  neither  CFD  or  tank  tests  went  so  far  as  deck-edge 
immersion,  so  the  boundary  of  the  experimental  data  would  be  reached. 

To  incorporate  sail  trim,  two  variables  are  used  in  the  calculation.  The  first  variable,  called 
“flat”  acts  as  if  the  camber  of  the  sail  was  being  reduced.  The  effect  of  this  is  that  coefficient  of 
lift  of  the  sail  is  reduced  proportionally.'^^^ 

The  second  variable,  called  “reef’  acts  as  if  the  exposed  sail  area  was  reduced.  The  reefing 
variable  affects  three  components.  First,  lift  is  reduced  proportionally  as  for  the  flat  variable. 
Secondly,  the  parasitic  drag  is  reduced  by  a  factor  of  reef  squared.  Finally,  as  if  truly  reefing  a 
sail,  the  vertical  center  of  effort  of  the  sail  is  reduced. 

Another  correction  to  calculate  the  actual  sail  forces  is  in  the  velocity  variable  (Eqn  (28)). 
From  the  solver  portion  of  a  VPP,  wind  speed  and  direction  are  given  as  true  wind  speed  and  true 
wind  angle.  However,  the  aerodynamic  functions  for  drag  and  lift  are  dependent  on  apparent 
wind  speed  and  apparent  wind  angle.  Furthermore,  the  velocity  of  the  wind  (Veff)  at  the  height 
of  the  center  of  effort  of  the  sails  over  the  water  (CEz)  is  not  the  same  as  the  true  wind  speed 
(Vtws)-  Because  of  friction  with  the  water,  at  low  altitudes,  the  true  wind  speed  is  reduced. 
True  wind  speed  can  be  found  using  Milgram’s  logarithmic  function  (Eqn  (32)). 

Veff  (CE.  )  =  (0.62 1 1 54  +  0. 1 0857(ln(CE  ^  )))  (32) 

Velocity  has  to  be  corrected  two  more  times  (Figure  (50)).  First,  the  true  wind  velocity  in  the  y- 
direction  must  be  factored  by  cosine  of  heel  in  order  to  find  the  velocity  in  the  heeled  plane. 


True  Wind  Angle 
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Figure  50  -  Velocity  diagram  in  the  xy-plane 

Finally,  the  vectors  from  boat  speed  and  wind  speed  are  added  to  create  an  apparent  wind  speed. 
The  boat  speed  vector  includes  both  speed  of  advance  and  leeway  angle.  Using  the  apparent 
wind  speed  in  Eqn(28),  the  total  aerodynamic  heelforce  and  drive  are  calculated. 

The  center  of  effort  of  each  sail  is  calculated  using  approximations.  For  the  mainsail,  the 
vertical  center  of  effort  (CEmz)  is  calculated  using  the  vertical  length  of  the  sail  (P),  the  height  of 
the  boom-above-deck  (BAD),  the  height  of  the  freeboard  at  the  mast  (FBD),  and  the  heel  angle 
(dJ).  The  0.39  value  is  an  approximation  that  the  center  of  effort  on  the  sail  is  39  percent  from 
the  bottom  of  the  sail  assuming  the  sail  is  semi-triangular  in  shape. 

CEm^  =  (0.39P  ■  reef  +  BAD  +  FBD)cos(0)  (33) 

The  center  of  effort  of  the  jib  (CEjz)  uses  the  same  calculation  as  that  for  the  mainsail, 
except  the  sail  is  assumed  to  start  at  the  deck,  not  at  a  certain  height  above  the  deck  (BAD). 
Furthermore,  the  vertical  length  of  the  jib  (I)  is  used  instead  of  the  luff  of  the  main  (P). 

CEj^  =  (0.39/  •  reef  +  FBD)cos(0)  (34) 

The  center  of  effort  of  the  spinnaker  (CEsz)  uses  a  calculation  also  similar  to  the  mainsail. 
Because  of  the  shape  of  the  spinnaker,  the  center  of  effort  is  assumed  to  be  50  percent  from  the 
bottom  of  the  sail.  However,  the  bottom  of  the  spinnaker  normally  begins  at  the  bottom  of  the 
mainsail  (BAD+FBD).  The  value  EHM  is  the  effective  height  of  the  mast  and  extends  from  the 
deck  to  the  top  of  the  mast  where  the  spinnaker  begins. 


CE^^  =  {0.5EHM  •  reef  +  BAD  +  FBD)cosiO) 


(35) 
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Furthermore,  the  longitudinal  centers  of  effort  of  each  sail  can  be  found  as  well.  For 
upwind  sailing,  the  angle  of  attack  on  the  mainsail  is  typically  around  15  degrees  while  the  jib  is 
at  12.5  degrees.  As  the  wind  moves  aft,  these  angles  increase  to  27  degrees  for  both  sails  when 
the  apparent  wind  is  at  90  degrees.  Finally,  both  sails  have  angles  of  attack  approximated  at  90 
degrees  when  the  wind  is  fully  aft  of  the  ship.  For  a  triangular  three  dimensional  foil  under 
typical  sailing  conditions,  the  longitudinal  center  of  effort  can  be  found  from  Figure  (51). 


Figure  51  -  Center  of  effort  for  a  three-dimensional  triangular  foil  with  moderate  camber 

Heel  force  and  drive,  however,  are  located  on  the  vertical  center  of  effort  of  the  sails 
(CEz)  above  the  datum  on  the  z-axis.  At  this  center  of  effort,  two  more  variables,  the  heeling 
moment  and  the  trimming  momentum  can  be  calculated  using  CEz  as  the  moment  arm.  The 
rolling  moment  is  the  heel  force,  renewed  at  an  axis  perpendicular  to  the  length  of  the  boat, 
multiplied  by  the  arm  CEz.  Heeling  moment  calculates  the  moment  which  heels  the  boat. 
Trimming  moment  is  renewed  at  an  axis  parallel  with  the  length  of  the  boat.  Trimming  moment 
calculates  the  moment  which  pushes  the  bow  of  the  boat  down. 
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By  multiplying  the  longitudinal  center  of  effort  moment  arm  of  each  sail  by  the  xy-planar 
forces  from  each  sail,  the  yawing  moment  due  to  aerodynamics  can  be  calculated. 

Hydrodynamic  and  Aerodynamic  Solution  -  Tank  Data 

Since  both  the  aerodynamic  and  hydrodynamic  forces  could  now  be  calculated,  a  solution 

had  to  be  found  to  equate  all  forces  in  order  to  maintain  equilibrium.  From  the  tank  data,  the 
drag  data  was  scaled  using  the  Froude  hypothesis  to  extrapolate  to  full-scale  resistance.  The 
resistance  could  be  equated  to  the  drive  variable. 

FroudeExtrapolation{Drag)  =  Drive  (36) 

Furthermore,  the  total  lift  measured  in  the  tank  was  the  summation  of  the  lift  measured  from 
the  forward  force  block  and  the  lift  at  the  aft  force  block.  This  summation  multiplied  by  the 
cubed  scaling  factor  (X)  could  be  equated  to  the  aerodynamic  sideforce. 

{Lift Aft  +  LiftFwd  )  •  =  Sideforce  (37) 

After  using  FASTSHIP  to  find  the  hydrostatic  righting  moment  of  the  ship,  it  was  coded  as 
an  automatic  visual  basic  macro  (RM(<1>))  which  provided  the  static  righting  moment  as  a 
function  of  heel.  However,  to  find  the  hydrodynamic  righting  moment,  three  approximations 
were  made.  First,  it  was  assumed  that  the  data  from  the  forward  lift  block  closely  approximated 
the  force  generated  by  the  keel.  Secondly,  the  data  from  the  aft  lift  block  approximated  the  force 
generated  by  the  rudder.  Third,  it  was  assumed  that  the  lift  of  the  keel  acted  on  a  point  at  Vi  the 
draft  of  the  keel,  and  the  lift  of  the  rudder  acted  on  a  point  at  Vi  the  draft  of  the  rudder.  These 
dynamic  functions  of  lift  acted  against  the  hydrostatic  righting  moment  of  the  vessel.  Therefore, 
in  summation,  the  righting  moment  could  be  equated  to  the  heeling  moment. 

- Liftfwd  •  Draft - Liftaft  •  ^ •  Draft =  Heelingmoment  (38) 

Finally,  the  yawing  moment  from  the  tank  data  was  the  force  from  the  aft  force  block 
multiplied  by  a  moment  arm  from  the  pin  location  to  the  heave  post.  This  moment  could  be 
equated  to  the  aerodynamic  yawing  moment. 


LiftAft  ■  Arm  ■  X^  =  Yawingmoment 


(39) 
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Hydrodynamic  and  Aerodynamic  Solution  -  SPLASH  Data 

In  general,  the  SPLASH  solution  was  identical  to  the  tank  solution.  Both  solutions  used  the 

same  aerodynamic  equations.  One  advantage  of  the  SPLASH  data  was  that  scaling  did  not  have 
to  occur.  However,  in  order  to  gain  better  correlation  to  the  tank  data,  the  same  viscous  stripping 
method  was  used.  The  viscous  component  of  SPLASH ’s  prediction  was  subtracted  from  its  total 
resistance  calculation.  The  viscous  prediction  using  the  ITTC  friction  equation  and  the  form 


factor  from  the  tank  were  used  instead.  The  equations  used  were: 

Resistance  =  Drive  (40) 

Lift  SPLASH  =  Sideforce  (41) 

-  Liftgpp^gu  •  y  •  Draft =  Heelingmoment  (42) 

Yawingmoment  ^pp^gfj  =  Yawingmoment  (43) 


VPP  Predictions 

Now  that  four  equations  were  made  to  satisfy  equilibrium,  a  prediction  code  was  necessary 
to  find  the  optimum  speed. 

Normally,  a  VPP  will  predict  conditions  which  create  the  greatest  boat  speed  for  a  given 
wind  angle  and  direction.  Some  modification  of  this  was  necessary  for  the  acquired  data  set. 
First,  the  Navy  44  VPP  included  yaw  angle  as  a  variable.  If  solving  an  upwind  condition  based 
upon  greatest  boat  speed  alone,  the  VPP  would  predict  extremely  high  yaw  angles  (which  would 
be  caused  by  large  amounts  of  sail  area  which  create  high  boat  speeds).  However,  the  high  yaw 
angle  would  not  mean  the  fastest  upwind  course,  or  highest  VMG.  Instead,  the  solver  was  made 
to  determine  the  greatest  VMG,  not  boat  speed  for  wind  angles  where  a  sailor  would  be  trying  to 
sail  upwind  or  downwind.  For  close  reaching  through  broad  reaching  wind  angles,  a  sailor 
would  be  content  with  the  maximum  boat  speed  since  VMG  is  not  the  goal.  Because  of  this, 
from  30  to  50  degrees  and  from  150  to  180  degrees  true  wind  angle,  the  solver  was  set  to  solve 
for  maximum  VMG,  not  boat  speed.  This  was  set  as  the  variable  “maxspeed.” 

In  solving  for  maxspeed,  only  two  variables  were  constant:  true  wind  speed  and  true  wind 
angle.  Initially,  a  method  of  predicting  the  optimum  condition  was  used  similar  to  PCSail’s 
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method.  This  involved  the  use  of  the  SOLVER  function  in  Excel.  By  holding  the  solver  to  the 
requirements  that  the  four  equilibrium  equations  must  not  exceed  0.05%  error,  one  variable, 
maxspeed,  could  be  maximized.  Items  changed  were  boat  speed,  yaw  angle,  heel  angle,  rudder 
angle,  reef,  and  flat. 

LFltimately  the  SOLVER  method  was  discontinued  from  use  in  the  code.  The  SOLVER 
code  made  predictions  which  were  unreasonable,  but  which  were  supported  by  the  data.  The 
problem  was  that  in  assembling  the  data  matrix,  allowances  were  not  made  on  which  headsail 
was  theoretically  flown.  In  tank  testing,  the  trim  weights  were  moved  which  simulated  the 
trimming  moment  as  a  result  of  either  the  jib  or  the  spinnaker.  In  Figure  (19),  these  conditions 
can  be  shown  as  clockwise  patterns  for  each  wind  condition.  The  jib  is  flown  first  at  the  lowest 
speed  (which  is  not  zero),  indicating  a  close-hauled  VMG  condition.  The  jib  is  flown  until  the 
third  following  point  where  it  was  predicted  the  spinnaker  would  be  flown  instead.  The 
spinnaker  continues  use  along  each  clockwise  pattern  until  a  heel  angle  of  zero  degrees  which  is 
assumed  to  be  a  downwind  VMG  course. 

The  solver  would  often  predict  heel  angles  for  upwind  conditions  at  downwind  heel  angles 
and  vice  versa.  Furthermore,  the  SOLVER  was  very  sensitive  to  any  dips  in  drag  for  a  given 
heel  angle  with  respect  to  constant  velocities.  The  SOLVER  would  often  predict  the  same  heel 
angle  for  a  large  range  of  speeds.  Finally,  with  four  sets  of  equations  as  well  as  large  matrices  of 
data  used  for  interpolation,  the  SOEVER  function  on  Excel  was  very  slow. 

A  new  method  was  created  which  used  circular  functions  in  Excel.  All  four  equilibrium 
conditions  produced  a  set  of  four  simultaneous  equations.  By  breaking  the  equations  into  their 
hydrodynamic  and  aerodynamic  components,  a  matrix  could  be  formed  which  represented  the 
equations.  The  four  tank  data  equilibrium  equations  in  matrix  form  are  shown  in  Eqn  (44). 

resistance  0  0  0  Drive 

0  -  Liftfwd  ■  X* Draft -  Liftaft  ■  X -jDraft^j  RM{^)HeelMoment 

0  Liftfwd  ■  X'  Liftaft  ■  X  0  Sideforce 

0  0  Liftaft  ■  Arm  ■  X  0  YawMoment 

By  taking  the  matrix  multiplication  of  the  inverse  hydrodynamic  matrix  with  the 
aerodynamic  matrix,  a  solution  matrix  can  be  formed.  This  solution  matrix  represents  the  factors 
by  which  resistance,  righting  moment,  fwd-lift,  and  aft-lift  should  change  in  order  to  come  to  an 
equilibrium.  The  SPLASH  data  matrix  was  identical  to  Eqn  (44)  with  corrections  made  using 
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Eqns  (40-43).  The  solution  for  the  SPLASH  data  represented  the  factors  by  which  resistance, 
righting  moment,  lift,  and  yawing  moment  should  change. 

Four  simplifications  were  made  in  the  iterative  approach.  First,  resistance  was  assumed  to 
be  a  function  mostly  of  velocity.  Secondly,  righting  moment  was  assumed  a  function  of  heel 
angle.  For  the  tank  data,  Fwd-lift  was  assumed  a  function  of  yaw,  and  aft-lift  was  assumed  a 
function  of  rudder  angle.  For  the  SPLASH  data,  lift  was  assumed  a  function  of  yaw,  and  yawing 
moment  a  function  of  rudder  angle.  By  modifying  these  four  variables  by  small  amounts  of  the 
solution  matrix,  the  spreadsheet  eventually  settled  on  a  solution  which  was  in  equilibrium.  This 
simple  finite  difference  method  was  written  directly  into  the  spreadsheet  cells  themselves.  Lach 
condition  took  at  most  5  seconds  to  solve. 

Consideration  in  the  finite  difference  solver  had  to  be  taken  in  the  movement  of  the 
variables.  If  an  initial  degree  of  freedom  (heel,  yaw,  rudder  angle,  or  velocity)  was  terribly 
unrealistic  for  wind  speed  and  wind  angle,  the  iterations  had  the  ability  to  self-destruct  and  some 
components  would  go  to  infinity.  This  would  stop  all  calculations.  Self-destruction  would  occur 
if  the  solution  matrix  steps  were  not  small  enough  to  fall  within  the  limits  of  the  data. 
Furthermore,  if  the  matrix  steps  were  too  large,  the  program  would  end  up  iterating  between  two 
conditions  forever,  never  falling  into  the  set  maximum  error  of  0.05%.  The  problem  with  setting 
the  steps  too  small,  however,  was  that  the  program  could  take  hours  to  run  a  complete  set  of 
conditions.  A  solution  was  found  in  that  the  step  size  would  decrease  per  number  of  iterations. 
If  the  code  were  to  self-destruct,  it  would  automatically  reset  and  start  again  at  a  new  initial 
positions  with  a  smaller  step  size  than  before. 

The  next  step  was  to  determine  how  the  sails  should  be  set  to  maximize  maxspeed.  By 
starting  at  a  predetermined  condition  of  reef  and  flat,  reef  and  fiat  would  be  changed  to  test  for 
greater  conditions  of  maxspeed.  For  wind  speeds  greater  than  18  knots,  it  was  found  that  reef 
and  flat  should  be  initially  set  to  0.8  for  both  variables.  At  mid-ranged  wind  speeds,  the  initial 
reef  and  flat  was  at  0.9.  For  wind  speeds  less  than  12  knots,  reef  and  flat  were  both  set  at  their 
maximum  values  of  1.0.  Finally,  the  spinnaker  was  always  flown  at  true  wind  angles  greater 
than  90  degrees. 

After  many  trials  with  the  new  VPP,  it  was  found  that  the  rudder  angle  was  not  being 
correctly  predicted.  This  was  due  fiilly  to  the  longitudinal  aerodynamic  center  of  effort 
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calculation.  While  this  calculation  approximated  the  center  of  effort  for  moderate  wind  speeds 
and  close-hauled  angles,  the  approximation  broke  down  as  reef  and  flat  changed  the  shape  and 
and  camber  of  the  sail.  Furthermore,  sailors  have  other  controls  to  balance  the  yaw  of  the  boat 
besides  rudder  angle.  Mast  rake  and  other  minor  sail  controls  can  significantly  change  the 
aerodynamic  yawing  moment  of  the  boat  without  much  loss  of  boat  speed.  The  solution  for  this 
VPP  was  that  this  fourth  equilibrium  condition  was  taken  out  of  the  solver.  Instead,  the  rudder 
position  was  set  at  0  degrees  while  reef  and  flat  iterated.  Once  a  final  sail  solution  was  found, 
the  rudder  angle  would  iterate  between  0  and  6  degrees  in  order  to  maximize  maxspeed. 
However,  the  yaw-moment  calculation  still  remained  in  the  code.  The  Yaw-moment  imbalance 
has  plagued  researchers  and  has  not  yet  been  adequately  solved.'"^  A  partial  solution  was  reached 
and  provided  the  preferred  mast  rake  position  which  is  valuable  information  for  mast  tuning. 


Discussion  with  Naval  Architect  Bill  Cook,  Oracle  Performance  Manager,  in  March  2003 


Performance  Prediction  Comparison 
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Polar  diagrams  are  plots  in  which  boat  speed  is  a  function  of  true  wind  angle  for  a  given 
wind  speed.  For  sailors,  polar  diagrams  help  to  tune  a  boat  and  to  find  the  quickest  course  to  a 
mark.  For  designers,  polars  help  provide  comparisons  between  boats  and  between  different 
configurations  of  the  same  boat.  Ultimately,  sets  of  polars  were  made  which  compared  the  IMS 
VPP,  PCSail,  tank  testing,  and  SPLASH  predictions.  Although  an  FKS  VPP  was  attempted,  its 
results  were  of  limited  use  since  lift  circulation  was  not  calculated  in  the  CFD  code.  With  four 
different  methods  at  four  different  wind  speeds,  it  became  difficult  to  interpret  all  predictions  on 
the  same  plot.  Therefore,  the  results  were  investigated  at  each  wind  speed. 
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Mk  II  Navy  44  SIC  Polar  Diagram 

6-knats  True  Wind  Speed 


■IMS 

■PCSail 

TANK 

SPLASH 
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The  6-knot  wind  speed  polar 
diagram  showed  the  most 
deviation  between  each 
prediction  method  for  a  given 
wind  speed.  Still,  the  average 
deviation  was  within  0.19  knots. 
In  the  6-knot  condition,  the  IMS 
VPP  was  assumed  to  have  the 
greatest  accuracy  because  of  its 
complex  aerodynamic  algorithm. 
Using  the  IMS  VPP  as  the 
baseline,  noticeably  lower  speed 
predictions  were  made  in  tank 
testing  particularly  at  deep 
reaching  angles.  This  was 
considered  error  due  to  scale 
effects  at  slow  model  speeds  and 
sub-turbulent  Reynolds  numbers. 

SPLASH,  however,  showed 
results  which  were  very 
consistent  with  the  IMS 
predictions  with  an  average 
deviation  of  less  than  0.08  knots. 
Even  though  the  aerodynamic 
algorithm  was  shared  by  the  tank 
VPP  and  the  SPLASH  VPP,  the 
interaction  of  the  hydrodynamics 
and  aerodynamics  in  the  solution 


Figure  52  -  Polar  diagram  for  6  knots  true  wind  speed 
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lead  to  a  chance  of  decreased  accuracy.  For  instance,  if  tank  testing  miscalculated  resistance  due 
to  scale  effects,  the  accuracy  of  the  aerodynamic  iteration  would  also  be  reduced. 

PCSail  showed  consistently  lower  velocity  predictions  for  this  windspeed.  Because  PCSail 
used  a  much  smaller  parametric  database  than  the  IMS  VPP,  these  predictions  deviated  both  due 
to  a  lack  of  sufficiently  populated  hydrodynamic  data  as  well  as  a  less  complex  aerodynamic 
functions.  Again,  the  iteration  of  the  hydrodynamic  and  aerodynamic  models  lead  to  errors  in 
this  low  wind  condition.  In  a  practical  sense  however,  the  lower  accuracy  of  the  prediction 
methods  may  not  be  too  significant  for  the  future  uses  of  the  Mk  II  Navy  44.  To  meet  Naval 
Academy  Sailing  schedules,  the  engine  is  usually  turned  on  when  boat  speeds  fall  below  five 
knots.  Furthermore,  races  are  rarely  started  until  wind  speeds  top  six-knots. 


Boat  Speed  [knots] 

CO  4:^  cn  O)  'sj  00  CD 
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Mk  II  Navy  44  SIC  Polar  Diagram 

12-knats  True  Wind  Speed 
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The  12-knot  wind  speed  data 
showed  excellent  correlation 
between  all  four  methods.  The 
average  deviation  between  all  of 
the  methods  was  less  than  0.08 
knots. 

The  results  at  this  wind 
speed  were  significant  in  that  12- 
knots  true  wind  speed  is  the 
normal  wind  speed  which  the  Mk 
11  Navy  44  STC  will  see  in 
operation.  In  both  the 

Chesapeake  Bay  and  near  the 
East  Coast  in  the  Atlantic  Ocean, 
the  Navy  44  will  often  find  itself 
sailing  in  12-knots  of  wind. 

The  12-knot  wind  condition 
showed  the  least  deviation  of  the 
experimental  data  because  the 
tank  testing  data  was  out  of  the 
scale  effects  region,  and  the 
SPLASH  data  was  unhindered 
by  any  discontinuities  in  the  free- 
surface  continuum. 


Figure  53  -  Polar  diagram  for  12  knots  true  wind  speed 
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II  Navy  44  STC  Polar  Diagram 

16-knats  True  Wind  Speed 
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At  16-knots  of  wind  speed, 
all  methods  predicted  very 
similar  results  upwind  (45-60 
degrees).  However,  as  the  boat 
turned  downwind,  SPLASH 
predicted  less  boat  speed.  At  the 
higher  speed  range,  SPLASH 
predicted  a  greater  resistance 
than  tank  testing  because  of  its 
different  method  of  viscous 
calculations.  The  tank  data 
mostly  fell  between  the  IMS  and 
the  SPLASH  predictions. 

Furthermore,  both  SPLASH 
and  tank  testing  predicted 
velocities  an  average  of  0.14 
knots  lower  than  the  IMS 
predictions.  Since  the  SPLASH 
and  tank  VPPs  solved  the  yaw 
balance,  this  error  was  found  to 
be  in  the  IMS  VPP  since  it  could 
not  calculate  the  specific  effects 
of  the  lifting  surfaces  (such  as 
vortices  and  wake  shed  from  of 
the  appendages). 

Overall,  at  16-knots  of  wind 


54  -  Polar  diagram  for  16  knots  true  wind  speed 


speed,  the  average  deviation  for 
the  methods  was  less  than  0.10 
knots. 
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Mk  II  Navy  44  SIC  Polar  Diagram 

20-kncts  True  Wind  Speed 
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At  20-knots  wind  speed,  the 
tank  and  SPLASH  VPPs 
continued  to  predict  boat  speeds 
generally  lower  than  the  IMS 
VPP.  The  lower  boat  speed  of 
SPLASH  as  compared  to  tank 
testing  near  beam-reaching 
angles  (85-105  degrees)  again 
reflects  the  differences  in  the 
viscous  calculations  used  by  both 
methods. 

While  PCSail  showed 
minimal  deviation  with  the  IMS 
prediction  at  downwind  courses, 
PCSail  overestimated  boat  speed 
in  close-reaching  conditions  (65- 
85  degrees).  Again  for  PCSail, 
this  error  was  due  to  not  having 
as  robust  a  hydrodynamic  dataset 
compared  to  the  IMS  VPP. 


Figure  55  -  Polar  diagram  for  20  knots  true  wind  speed 
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Performance  Prediction  Conclusions 

Ultimately,  writing  the  new  customized  Mk  II  Navy  44  VPP  was  the  only  way  to  compare 
different  methods  of  performance  prediction.  The  end  result  was  to  say  that  the  customized  VPP 
provided  very  good  correlation  to  both  the  PCSail  code  and  IMS  code  in  terms  of  relative 
performance  between  wind  angles  and  wind  speeds.  The  customized  VPP  also  made  rudder  and 
yaw  angle  predictions  which  neither  the  PCSail  or  IMS  VPP  calculated.  Therefore,  the 
customized  VPP  made  a  more  complete  performance  prediction  which  also  provided  more  input 
to  sailors  in  how  to  sail  the  Mk  II  Navy  44. 

In  comparing  the  hydrodynamic  predictions,  tank  testing  was  the  easiest  method  used  to 
evaluate  the  performance  of  the  boat.  Tank  testing  had  the  benefit  of  many  years  of  previous 
research  to  help  fine  tune  the  performance  prediction  process.  However,  tank  testing  required 
many  hours  spent  in  the  tow  tank.  The  hours  spent  testing  many  conditions  meant  that  most 
conditions  were  tested  only  once  or  twice.  There  was  no  possible  analysis  for  evaluating  the 
precision  of  each  run. 

Both  CFD  codes  had  the  advantage  of  having  100%  precision  in  their  tests.  FKS,  however, 
proved  that  having  100%  precision  was  meaningless  if  the  accuracy  was  insufficient.  FKS 
showed  itself  to  be  an  excellent  tool  for  predicting  the  upright  performance  of  vessels.  Once  the 
meshed  hull  was  built,  a  full  range  of  upright  tests  took  twenty  minutes.  This  compared  very 
favorably  to  a  week’s  worth  of  tank  testing. 

Because  of  the  far-field  calculations  used,  FKS  was  shown  to  be  very  sensitive  to  inviscid 
effects  and  assumptions  from  slender-ship  theory.  This  resulted  in  the  large  humps  and  hollows 
in  the  coefficient  of  resistance  (Figure  (34)).  FKS  was  not  written  to  calculate  the  effects  of 
lifting  surfaces,  therefore  only  upright  conditions  were  possible.  In  the  end,  FKS  was  shown  to 
be  an  excellent  numerical  approximation  of  a  closed-form  inviscid  resistance  calculation. 

SPLASH  showed  impressive  accuracy  compared  to  tank  testing  and  the  IMS  predictions. 
The  main  advantage  of  SPLASH  was  that  the  hull  was  free  to  move  and  pivot  in  any  direction 
based  on  aerodynamic  and  hydrodynamic  forces.  This  allowed  for  an  improved  mesh  over  the 
hull  and  therefore  more  accurate  results. 

The  disadvantage  to  SPLASH  was  from  a  user’s  perspective.  Being  such  a  complex  code, 
SPLASH  was  difficult  to  make  operational.  However,  once  working,  SPLASH  was  very  fast  to 
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run.  Whereas  the  full  sailing  matrix  took  weeks  to  perform  in  tank  testing,  a  similar  matrix  took 
nineteen  hours  to  calculate  with  SPLASH. 

Based  on  all  the  research,  the  performance  prediction  process  for  the  Mk  II  Navy  44  STC 
was  completed  by  making  a  final  set  of  polars. 


Polar  Diagram  for  the  Mk  II  Navy  44  STC 


180 


By  taking  input  from  each 
method,  a  best-fit  spline  was 
constructed  for  each  wind  speed. 
In  some  areas,  the  IMS  VPP  was 
deemed  the  most  accurate 
(especially  at  6-knots  wind  speed 
because  of  its  complex 
aerodynamic  code).  However, 
for  higher  boat  speeds,  tank 
testing  and  SPLASH  conditions 
were  always  incorporated 
because  of  their  highly  accurate 
hydrodynamic  data  (Figure  (56)). 


Figure  56  -  Polar  Diagram  for  Mk  II  Navy  44  STC 
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Given  the  inherent  errors  present  in  scale  effects  and  flow  simulation,  the  use  of  tank  testing 
to  setup  and  normalize  CFD  codes  would  produce  the  prediction  with  the  least  amount  of  error. 
For  sail  craft,  once  a  hydrodynamic  prediction  is  made,  its  VPP  should  be  used  alongside 
predictions  from  other  sources  for  a  full  estimate  of  the  performance  of  the  vessel. 

Once  a  prototype  Mk  II  Navy  44  STC  is  built,  full-scale  trials  could  validate  some  of  the 
results,  although  the  setup  and  testing  for  full-scale  tow  trials  would  be  difficult.  In  the  short 
term,  another  option  to  add  to  the  data  set  would  be  to  conduct  tank  testing  of  a  model  twice  or 
even  three  times  the  size  of  the  model  used  in  this  study.  While  taking  more  time  to  run  the  tests, 
this  would  improve  the  accuracy  of  the  tank  testing  data  and  would  reduce  the  scale  effects 
present  in  the  results. 

A  further  option  would  be  to  conduct  CFD  tests  using  a  viscous  code.  While  these  codes 
take  substantially  longer  to  run,  the  flow  modeling  on  the  ship  would  be  improved. 

Finally,  although  not  a  prediction,  an  accurate  GPS  along  with  months  of  sailing  on  the 
prototype  would  record  the  optimum  sailing  speeds  and  angles.  This  technique  would  validate 
the  predictions  made  in  the  design  phase. 


Rudder  Redesign 
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Since  SPLASH  had  shown  its  accuracy  in  the  previous  seetion  of  performance  predietion,  it 
was  used  to  help  perform  a  rudder  redesign.  Rudder  experimentation  is  an  area  in  yacht  design 
with  limited  development.  Indeed,  besides  wind  tunnel  testing  with  very  large  appendages  for 
lACC  yaehts,  most  sailboat  rudders  are  designed  with  little  experimental  data.  While  most  of 
the  hull  design  of  the  Mk  II  Navy  44  STC  is  nearly  finalized,  the  rudder  design  is  not  elose  to 
eompletion.  At  this  stage  of  the  design  of  the  Navy  44,  designer  David  Pedrick  has  simply  added 
a  generie  rudder  found  in  other  typical  sailboats.  By  completing  a  systematic  rudder  design  for 
the  Mk  II  Navy  44  STC,  the  knowledge  of  rudder  development  ean  be  improved. 

There  were  three  eomparisons  made  in  the  rudder  redesign: 

1)  Tank  testing  versus  SPLASH 

2)  Various  planform  shapes 

3)  Rudder  loeation  and  depth 

The  tank  testing 
eomparison  evaluated  two 
rudders  in  SPLASH  and  in  the 
tank.  The  first  rudder,  named 
the  “Pedriek”  rudder,  was  the 
rudder  which  was  supplied  by 
naval  arehitect  David  Pedriek 
with  the  original  lines  of  the 
Mk  II  Navy  44.  Bill  Beaver  of 
the  USNA  Laboratory  Support 
Staff  built  the  seeond  rudder, 
named  the  “Beaver”  rudder. 
Bill  Beaver’s  design  was  of  an 
elliptical  rudder  far  larger  than 
the  Pedrick  rudder. 

The  Beaver  rudder  was  positioned  at  the  same  rudder  stoek  location  as  the  Pedriek  rudder. 


Figure  57  -  Planforms  of  Pedrick  (left)  and  Beaver  (right)  rudders 
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The  Planform  testing  evaluated  SPLASH ’s 
capabilities  at  detecting  changes  in  the  planform  of  a 
rudder.  There  were  six  rudders  used  for  comparison. 
The  first  rudder  was  the  Pedrick  rudder.  The  second 
rudder  was  named  the  “Baseline”  rudder  as  it  was 
very  similar  to  the  Pedrick  rudder  with  slight 
modifications  to  the  tip.  The  tip  of  the  Baseline 
rudder  was  faired  and  the  base  of  the  tip  was  made 
flat  instead  of  rounded  as  in  the  Pedrick  rudder.  The 
Baseline  rudder  acted  as  the  starting  point  for  the  rest 
of  the  modified  rudders  of  these  series.  The  Baseline 
rudder  maintained  the  same  rudder  shaft  location  as 
the  Pedrick  rudder.  Furthermore,  for  the  entire  series 
of  rudders,  the  section  type  of  each  rudder  was  kept 
the  same  as  the  Pedrick  rudder. 


There  were  four 
modifications  of  the 
Baseline  rudder.  The 
“Tip”  rudder  was  the 
Baseline  rudder  with  the 
tip  extended  aft  to 
investigate  fluid  flow  at 
the  end  condition.  The 
“Bulge”  rudder  was  the 
Baseline  rudder  modified 
so  that  the  mid  portion  of 
the  rudder  was  1.3  times 
the  chord  of  both  the  root 
and  tip.  The  Bulge 
rudder  was  constructed  to  investigate  the  effects  of  fluid  flow  at  the  center  of  the  span. 


Figure  59  -  Planforms  of  the  Tip  (left)  and  Bulge  (right)  rudders 
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The  “Elliptical”  rudder 
was  the  Baseline  rudder 
which  was  faired  into  a  semi¬ 
ellipse.  Elliptical  planforms 
have  a  tradition  of  use  as 
wings  since  they  reduce 
induced  drag.  Furthermore, 
elliptical  rudders  maintain  a 
near  constant  coefficient  of 
lift  across  the  span.  This 
creates  even  stalls  across  the 
entire  rudder.  The  “Zoid” 
rudder  was  a  trapezoid  rudder 
of  the  same  dimensions  as  the 
elliptical  rudder.  Trapezoidal 
planforms  for  wings  have  similar  lift  and  drag  characteristics  as  elliptical  rudders  but  are  easier 
to  construct. 

The  location  and  depth  testing  included  four  rudders.  These  tests  were  conducted  to 
investigate  effects  of  aspect  ratio  as  well  as  location  of  the  rudder  relative  to  the  keel.  Often,  the 
vortices  from  the  keel  can  act  against  the  lift-producing  effects  of  the  rudder.  However,  if  a 
rudder  is  positioned  too  far  aft,  it  can  be  exposed  above  the  free  surface  and  lose  its  efficiency. 
Again,  the  Baseline  rudder  was  used  as  the  baseline  planform  for  the  depth  testing  as  well  as  the 
baseline  location.  The  “Maxdepth”  rudder  was  a  scaled  version  of  the  baseline  rudder  where  the 
tip  of  the  rudder  came  to  0.85  of  the  draft  of  the  keel.  This  is  the  maximum  draft  the  rudder 
could  be  with  respect  to  good  design  practice.  If  any  longer,  rudder  damage  due  to  grounding 
would  become  more  probable.  The  “Halffwd”  rudder  was  a  movement  of  the  Baseline  rudder  to 
a  position  one  half  of  a  span  length  forward  of  its  normal  position.  The  “Onefwd”  rudder  was  a 
movement  of  the  Baseline  rudder  to  a  position  one  full  span  length  forward  of  its  normal 
position.  Both  the  Halffwd  and  Onefwd  rudders  had  their  depths  held  at  the  same  depth  of  the 
Maxdepth  rudder.  Figure  (61)  shows  the  Baseline  rudder  interposed  in  the  pictures  of  the 
Halffwd  and  Onefwd  rudders. 


Figure  60  -  Planforms  of  the  Elliptical  and  Zoid  rudders 
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Figure  61  -  Rudder  planforms  for  Maxdepth,  Halffwd,  and  Onefwd,  respectively 


Tank  Comparison 

The  two  rudders  tested  both  in  SPLASH  and  in  the  tank  were  the  Pedrick  and  Beaver 
rudders.  For  both  methods,  three  yaw  angles  were  tested  of  0,  4,  and  6  degrees,  with  four  rudder 
angles  of  0,  2,  4,  and  8  degrees.  All  tests  were  condueted  with  zero  heel. 

Sinee  aetual  rudders  were  being  evaluated,  not  just  shapes  of  rudders,  the  requirement  was 
made  that  the  better  rudder  had  two  qualities:  a  higher  lift  to  drag  ratio  and  a  higher  yaw-moment 
to  drag  ratio.  The  higher  lift  to  drag  ratio  would  provide  faster  sailing  while  on  a  eonstant  course 
since  the  lift  of  the  rudder  acts  against  some  of  the  aerodynamic  sideforce.  The  higher  yaw- 
moment  to  drag  ratio  would  provide  for  faster  boat-handling  and  turning.  This  is  important  for 
inshore  races  where  maneuvering  is  very  important.  Likewise,  as  a  training  vessel,  the  ability  to 
turn  the  Mk  II  Navy  44  STC  would  add  to  the  safety  of  the  crew. 

In  the  original  assumptions  about  appendages,  yaw  was  assumed  to  have  a  linear  affect  on 
both  lift  and  drag,  while  rudder  angle  was  assumed  to  have  a  quadratic  affect.  If  plotting  Lift 
versus  Drag  in  which  both  yaw  and  rudder  change,  the  effect  would  be  cubic.  Figure  (62)  shows 
the  Pedrick  and  Beaver  rudders  tested  both  in  SPLASH  and  in  the  tank.  A  general  cubic  fit  of 
the  data  was  observed,  validating  these  assumptions. 
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SPLASH  and  tank  testing  analysis  showed  the  same  trend  for  both  rudders.  The  Pedrick 
rudder  produced  more  lift  for  less  drag  at  low  lift  conditions,  but  the  opposite  was  true  at  high  lift 
conditions.  This  can  be  further  described  as  that  the  Pedrick  rudder  was  more  efficient  at  low 
angles  of  attack  while  the  Beaver  rudder  was  more  efficient  at  high  angles  of  attack. 


Figure  62  -  Lift  versus  drag  comparing  the  Pedrick  and  Beaver  rudders 


This  data  was  applied  numerically  through  the  cubic  regression.  Two  new  curves  (Figures 
(63),  (64))  were  made  of  the  drag  of  the  Pedrick  rudder  subtracted  from  the  drag  of  the  Beaver 
rudder  over  a  range  of  lift  for  both  experimental  methods. 


Figure  63  -  Drag  delta  of  Beaver  to  Pedrick  in  SPLASH  Figure  64  -  Drag  delta  of  Beaver  to  Pedrick  in  tank 
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Tank  testing  predicted  a  greater  advantage  of  drag  in  low  lift  conditions  for  the  Pedrick 
rudder  than  in  SPLASH.  However,  some  error  in  the  low  lift  conditions  of  the  Beaver  rudder 
were  evident  from  the  data.  Furthermore,  the  inherent  scale  effect  errors  described  earlier  are 
amplified  when  using  small  appendages  leading  to  the  conclusion  that  SPLASH  produces  more 
reliable  lift  and  drag  predictions  especially  for  small  angles  of  attack.  Both  methods  still 
predicted  a  cross-over  point  where  the  Beaver  rudder  gains  advantage  over  the  Pedrick  rudder  in 
the  low  lift  condition. 

In  plotting  Yaw-moment  versus  Drag,  the  information  desired  was  mainly  how  the  rudder 
angle  affects  the  yawing  moment.  Therefore,  the  effect  of  rudder  on  the  boat  for  an  averaged 
yaw  was  quadratic. 


Figure  65  -  Yaw-moment  versus  drag  comparing  the  Pedrick  and  Beaver  rudders 


Figure  (65)  shows  that  the  Beaver  rudder  had  less  drag  for  the  same  amount  of  yaw  moment. 
SPLASH  and  tank  testing  showed  an  increase  of  the  Beaver  rudder’s  advantage  for  high  angles 
of  attack  of  the  rudder.  Therefore,  the  Beaver  rudder  would  have  reduced  drag  for  a  given  yaw- 
moment  against  the  Pedrick  rudder,  leading  to  faster  maneuvers.  Furthermore,  upwind  sailing 
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would  be  faster  with  the  Beaver  rudder  since  the  yaw  moment  of  the  rudder  is  applied  against  the 
aerodynamic  yaw-moment.  This  would  allow  for  more  power  to  be  produced  by  the  mainsail, 
leading  to  more  aerodynamic  drive,  and  therefore  more  boat  speed.  However,  a  note  of  caution 
is  that  all  rudder  tests  took  place  at  a  position  of  zero  degrees  heel.  To  fully  quantify  the  sailing 
performance  of  the  Beaver  rudder  a  full  sailing  spectrum  of  heel  and  yaw  in  SPLASH  should 
take  place. 

Planform  Comparison 

As  presented  in  the  tank  comparison  section,  tank  testing,  while  providing  general  trends  in 
appendage  performance,  does  not  accurately  predict  minute  changes  for  small  models. 
Therefore,  SPLASH  was  used  exclusively  to  determine  the  general  characteristics  of  different 
planforms  of  rudders.  Furthermore,  appendage  testing  in  CFD  required  roughly  one  half  hour  to 
build  each  rudder  in  a  computer-aided  design  program,  and  forty-five  minutes  to  complete  the 
testing  matrix  in  SPLASH.  For  model  design,  at  least  two  days  were  required  to  build  an 
appendage  followed  by  four  to  five  hours  of  tank  testing. 

The  object  of  the  planform  comparison  was  to  ascertain  different  qualities  of  the  shape  of 
each  rudder.  The  comparison  was  of  the  effectiveness  of  a  rudder  type,  not  just  the  best  rudder 
for  a  certain  boat.  Therefore,  a  desired  lift  or  yaw  moment  was  not  a  necessary  test.  The  rudder 
size  differences  were  negligible  as  long  as  they  were  all  within  limits  of  not  creating  excess 
wavemaking  drag.  Instead  of  lift,  drag,  or  yaw  moment,  nondimensionalized  coefficients  of  the 
respective  terms  were  used.  This  allowed  for  the  evaluation  of  different  sizes  easily. 

Figure  (66)  shows  the  results  of  the  efficiency  of  a  rudder  with  respect  to  lift.  SPLASH  was 
able  to  differentiate  between  all  of  the  shapes  even  though  the  average  performance  of  each 
rudder  (shown  by  the  cubic  fit)  were  almost  identical. 

Figure  (67)  shows  the  results  of  the  efficiency  of  a  rudder  with  respect  to  yaw  moment. 
There  was  more  visible  difference  in  performance  with  respect  to  turning  the  boat. 
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Figure  66  -  Lift  versus  drag  for  different  planforms 
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Figure  67  -  Yaw  moment  versus  drag  for  different  planforms 
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Improved  resolution  was  found  for  both  Figures  (66)  and  (67)  by  creating  plots  of  the 


difference  in  drag  coefficient  between  the  Baseline  rudder  and  the  other  rudders. 


-  Ellipse 

-  Zoid 

Figure  68  -Drag  delta  with  respect  to  lift  of  planforms  to  the  Baseline  rudder 

Figure  (68)  indicated  that  adding  a  tip  to  the  trailing  edge  of  the  rudder  decreased  the  drag. 
The  tested  geometry  of  the  Pedrick  rudder,  however,  clearly  indicated  decreased  performance 
from  the  Baseline  rudder.  The  Bulge  rudder  showed  advantage  for  small  amounts  of  lift  while 
the  opposite  was  true  for  the  Zoid  rudder.  From  these  results,  sailing  performance  would  be 
improved  downwind  with  the  Bulge  rudder  where  little  lift  is  needed,  but  would  be  improved 
upwind  using  the  trapezoidal  planform  in  high-lift  conditions.  No  general  trend  could  be  made 
of  the  elliptical  rudder,  except  that  in  only  a  very  limited  window  was  it  better  than  the  Baseline 
rudder.  The  Zoid  rudder  was  theorized  to  have  the  same  characteristics  of  the  elliptical  rudder. 
However,  the  elliptical  rudder  showed  a  very  wide  range  of  performance  while  the  Zoid  rudder 
stayed  consistently  close  to  the  baseline. 
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Figure  69  -  Drag  delta  with  respect  to  yaw  moment  of  planforms  to  the  Baseline  rudder 

Barely  visible  along  the  ordinate  axis  in  Figure  (69)  was  the  Zoid  rudder,  indicating  that  this 
rudder  was  almost  identical  to  the  Baseline  rudder  with  respect  to  the  efficiency  in  yawing  the 
boat.  Of  increasing  efficiency  in  yawing  the  boat  from  the  Baseline  there  were  two  rudders.  The 
Tip  rudder  was  the  next  most  efficient  followed  by  the  Bulge  rudder.  This  shows  that  increasing 
the  longitudinal  distribution  of  surface  area  increases  turning  performance.  The  Pedrick  rudder 
showed  more  efficiency  in  turning  the  boat  than  the  Baseline  rudder  for  small  amounts  of  yaw 
moment  (corresponding  to  small  rudder  angles).  However,  as  rudder  angle  increased,  the 
Pedrick  rudder  lost  efficiency  when  compared  to  the  Baseline.  Finally,  the  Elliptical  rudder 
showed  a  large  increase  in  efficiency  very  early,  but  also  lost  its  efficiency  as  rudder  angle 
increased,  much  like  the  Pedrick  rudder.  The  yaw  moment  calculations  for  the  Elliptical  rudder 
and  the  Pedrick  rudder  indicate  an  error  in  their  programmed  geometries  because  drag  rose 
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considerably  for  both  rudders  at  high  yaw  moments.  In  either  case,  high  yaw  moments  will 
accentuate  flow  patterns  around  the  leading  edge  of  the  rudder  and  in  the  induced  drag  vortex 
formed  at  the  tip.  However,  a  numerical  integration  of  forces  does  not  provide  a  method  for 
analyzing  the  actual  cause  of  the  problem. 

Explanations  for  the  differences  in  performance  could  only  be  found  in  plots  of  the  actual 
flow  patterns  across  each  rudder.  For  the  maximum  lift  condition  (6  degrees  yaw  and  8  degrees 
rudder),  where  the  differences  between  rudders  was  the  most  extreme,  SPLASH  plots  were 
constructed  of  both  the  back  and  face  surfaces  of  each  rudder.  The  face  surface  is  the  side  of  the 
rudder  which  contains  the  static  pressure  point  and  typically  has  positive  pressure  on  all  parts. 
The  back  surface  is  the  side  of  the  rudder  which  generates  most  of  the  suction  on  the  rudder  and 
contains  mostly  negative  pressures.  For  each  plot,  the  pressure  coefficient  as  well  as  velocity 
vectors  were  placed  on  top  of  the  rudder.  A  yellow  color  indicates  zero  pressure,  while  red 
indicates  positive  pressure.  Green  indicates  some  negative  pressure  while  blue  indicates  a  large 
amount  of  negative  pressure. 


Figure  70  -  Back  surface  of  Pedrick,  Baseline,  Tip,  Bulge,  and  Elliptical  Rudders  (direction  of  flow  is  left) 


The  vector  field  helped  to  indicate  flow  variations  between  the  rudders.  All  rudders  showed 
signs  of  induced  drag  vortices  at  the  tip.  However,  the  Tip  rudder  had  a  much  smaller  vector 
which  was  also  at  a  distance  from  the  rest  of  the  lifting  surface.  In  the  same  manner,  the  Bulge 
rudder  had  the  main  portion  of  its  lifting  surface  positioned  at  a  distance  away  from  the  vortex. 
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However,  the  Pedrick  rudder  showed  the  vortex  vector  moving  almost  vertically  up  the  span. 
The  Elliptical  rudder  showed  the  same  circumstance.  The  effect  of  having  a  vortex  at  a  distance 
away  from  the  section  increased  the  local  coefficient  of  lift.  Furthermore,  the  vector  fields  show 
a  flow  reversal  at  the  tip  of  the  Pedrick  rudder.  Because  of  the  flow  reversal,  the  tip  geometry  of 
the  Pedrick  rudder  was  seen  to  have  an  adverse  affect  on  the  efficiency  of  the  rudder  by  causing 
a  much  stronger  induced  drag  vortex.  A  closer  look  at  the  coded  grid  geometry  of  the  Pedrick 
rudder  showed  discontinuous  surface  gradients  at  the  tip,  which  proved  to  be  the  error. 

The  contour  plot  of  the  coefficient  of  pressure  also  indicated  the  effect  of  the  flow 
variations.  All  rudders  showed  large  low  pressure  near  the  leading  edge  at  the  root.  The  large 
low  pressure  signaled  separated  flow  and  indicated  that  the  flow  across  the  rudder  which  was 
next  to  the  hull  suffered  from  a  negative  pressure  gradient.  However,  both  the  Elliptical  rudder 
and  the  Pedrick  rudder  showed  separation  at  the  trailing  edge  near  the  tip.  This  loss  of  pressure 
was  also  due  to  a  negative  pressure  gradient  as  a  function  of  the  section  curve  at  the  tip.  As  a 
design  consideration,  this  negative  pressure  gradient  should  not  occur  at  the  tip  of  the  rudder  and 
was  due  to  improper  geometry.  Interestingly,  the  tip  of  the  Tip  rudder  showed  very  high 
pressure  as  a  result  of  the  strong  vortices  which  it  extended.  Finally,  along  the  leading  edge  of 
the  Elliptical  rudder,  separation  was  seen  at  Va  of  the  span  from  the  root.  This  separation 
indicated  that  the  leading  edge  of  the  rudder  at  this  point  was  too  thick  to  allow  for  a  normal, 
favorable  pressure  gradient. 
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Figure  71  -  Face  surface  of  Pedrick,  Baseline,  Tip,  Bulge,  and  Elliptical  Rudders  (direction  of  flow  is  right) 


The  face  surface  had  limited  use  in  qualifying  the  flow  effects  across  the  rudder  since  the 
flow  across  the  face  should  not  suffer  from  separated  flow  or  any  other  negative  flow 
characteristics.  However,  the  Elliptical  rudder  showed  vortices  which  had  an  unusual  bend 
across  the  face.  This  too  indicated  the  poor  construction  of  the  tip  geometry  of  the  Elliptical 
rudder. 


Location  and  Depth  Comparison 

The  location  and  depth  comparison  experimented  with  moving  and  extending  the  rudder  to 
different  positions.  The  Baseline  rudder  was  again  used  to  standardize  the  results  and  typical  lift 
and  yaw  moment  plots  were  constructed. 
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Figure  72  -  Lift  versus  drag  for  different  rudder  locations 
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Figure  73  -  Yaw  moment  versus  drag  for  different  rudder  locations 


The  lift  plot  resolution  was  improved  to  show  interesting  results  in  Figure  (72). 
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Figure  74  -Drag  delta  with  respect  to  lift  of  different  locations  to  the  Baseline 

Figure  (74)  shows  that  the  Baseline  rudder  had  good  placement  compared  to  the  Onefwd, 
Halffwd,  and  Maxdepth  rudders.  The  Onefwd  as  well  as  the  Maxdepth  rudders  clearly  increased 
drag.  The  Halffwd  rudder  showed  improvement  on  the  Baseline  rudder  for  small  amounts  of  lift, 
but  became  less  efficient  for  larger  lift  forces. 

All  of  the  rudders  suffered  from  flow  patterns  originating  at  the  keel.  As  the  keel  produced 
more  lift,  larger  vortices  formed  and  disrupted  the  flow  past  the  rudder.  The  forward  moving 
rudders  were  most  affected  at  these  high  lift  forces.  The  Maxdepth  rudder  suffered  from  this 
effect  as  well,  as  the  rudder  tip  was  deep  enough  to  by  influenced  by  the  keel’s  vortex. 
Furthermore,  the  Maxdepth  rudder  did  not  show  any  advantages  to  maximizing  the  rudder’s 
aspect  ratio.  The  additional  wetted  surface  area  also  increased  drag  and  would  hinder 
performance  downwind. 
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Finally,  the  Halffwd  rudder  showed  a  noticeable  improvement  in  performance  at  conditions 
of  low  lift.  At  these  positions,  the  Halffwd  rudder  was  more  submerged  than  the  Baseline  rudder 
which  resulted  in  increased  static  pressure  on  the  rudder  as  well  as  less  wavemaking  drag. 
However,  these  advantages  only  lasted  to  a  point  where  the  keel  vortex  gained  enough  strength 
to  disrupt  the  flow  over  the  Halffwd  rudder. 

Without  need  for  increased  resolution.  Figure  (73)  shows  the  efficiency  of  each  rudder  in 
regards  to  yawing  the  boat.  Because  of  a  decreased  moment  arm,  both  the  Onefwd  and  Halffwd 
rudders  had  increased  drag  in  turning  the  boat.  However,  the  Maxdepth  rudder  was  just  slightly 
more  efficient  than  the  Baseline  rudder  due  to  its  larger  and  more  submerged  area. 

Rudder  Redesign  Conclusions 

There  were  three  advantages  in  using  SPLASH  for  the  rudder  redesign  versus  tank  testing. 
First,  SPLASH  combined  with  the  FASTSHIP  CAD  program  took  less  than  two  hours  to 
complete  each  comparison  whereas  tank  testing  took  two  days.  Secondly,  SPLASH  did  not 
suffer  from  scale  effects  as  did  tank  testing  with  a  small  rudder.  Finally,  the  visualization  of  the 
flow  over  the  rudder  helped  to  analyze  each  design. 

The  planform  tests  showed  that  two  very  important  elements  of  rudder  design  evaluation  are 
having  a  quality  construction  of  the  geometry,  and  secondly,  locating  the  rudder  vortex  away 
from  the  main  lift  generating  surface.  Having  a  flattened  tip  of  the  Baseline  could  be  considered 
an  advantage  from  the  rounded  tip  of  the  Pedrick  rudder.  However,  the  Pedrick  rudder’s 
geometry  had  minor  programming  problems,  and  a  reliable  comparison  of  its  tip  shape  cannot  be 
made.  The  location  and  depth  test  showed  that  moving  the  rudder  away  from  the  keel  vortex 
also  improves  efficiency. 

In  regards  to  the  Mk  II  Navy  44’ s  rudder,  the  research  produced  two  recommendations. 
First,  the  rudder  should  be  designed  to  minimize  the  induced-drag  vortices’  interference  with  the 
lifting  surfaces  of  the  rudder.  Secondly,  the  rudder  should  remain  at  the  same  position,  and  its 
depth  should  not  be  increased. 

Therefore,  the  design  of  the  rudder  for  the  Mk  II  Navy  44  STC  should  be  either  similar  to 
the  Tip  or  Bulge  rudder,  or  even  a  combination  of  the  two.  These  rudders  provided  the  most 
efficiency  in  regards  to  both  lift  and  yawing  moment.  Full-scale  testing  using  different  rudders 
however  would  be  the  only  method  of  verifying  these  results. 
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Recommended  Performance-Based  Design  Approach 

From  this  study  of  performance  prediction  methods,  a  refined  design  process  was  created  for 
the  performance-based  design  of  sailboats.  In  a  full  ship  design,  the  performance  prediction 
phase  would  only  be  one  step  in  the  complete  design  spiral.  Using  the  tools  of  tow  tank  testing, 
a  simple  parametric  VPP  (PCSail),  a  simple  CFD  code  (FKS),  and  a  complex  CFD  code 
(SPLASH),  an  accurate  method  of  predicting  performance  can  be  made  through  iteration  in  the 
design. 

1 .  From  the  design  requirements,  obtain  the  initial  principal  dimensions  of  the  hull,  sail  plan, 

and  center  of  gravity. 

2.  Use  PCSail  to  find  the  target  boat  speed  and  hull  form  coefficients. 

3.  Develop  preliminary  hull  lines  from  the  hull  form  coefficients. 

4.  Run  upright  conditions  in  FKS  for  a  range  of  Froude  numbers  from  0.05  to  0.7  to  get  the 

effective  horsepower  (EHP).  From  the  EHP,  make  an  initial  engine  selection.  Also, 
determine  the  tow  tank  drag  and  lift  force  block  sizes  from  force  estimates.  Refine  the 
hull  lines  based  on  the  predicted  wave  profile  from  FKS. 

5.  Construct  a  scale  model  with  a  basic  keel  and  rudder. 

6.  Run  upright  tank  tests  with  and  without  turbulence  stimulators  for  the  likely  speeds  predicted 

by  PCSail,  emphasizing  speeds  corresponding  to  a  model  Reynolds  number  of  less  than 

1,000,000. 

7.  Determine  the  form  factor  using  the  sand  strip  correction  factor. 

8.  Run  tank  tests  at  6,  16,  and  24  knots  of  true  wind  speed  at  true  wind  angles  of  40,  90,  and  120 

degrees  using  heel  and  speed  estimates  from  PCSail.  Include  yaw  angle  tests  of  0  and  5 
degrees  and  rudder  angles  of  0,  3,  and  6  degrees.  Record  drag,  lift,  and  yaw  moment. 

9.  Standardize  SPLASH’ s  trimming  moment  estimations  with  the  tank  testing  data. 

10.  Run  SPLASH  sailing  predictions  until  the  hull  lines  provide  the  desired  performance 

A.  Utilize  SPLASH  for  the  full  sailing  matrix  at  wind  speeds  of  6,  12,  18,  and  24  knots 
and  wind  angles  of  40,  60,  90,  120,  and  170  degrees.  Include  yaw  angles  of  0,  3, 
and  5  degrees  and  rudder  angles  of  0,  3,  and  6  degrees. 


B.  Export  updated  lines  to  a  hydrostatic  program  to  obtain  the  righting  moment. 
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C.  Set  the  SPLASH  data  into  the  custom- VPP  and  ran  the  full  polar  prediction 

1 1 .  Evaluate  high  pressure  and  large  vortex  regions  on  the  appendages  in  SPLASH  and  create 

modified  appendages. 

12.  Run  SPLASH  appendage  predictions  until  satisfied  the  upwind  and  downwind  performance. 

A.  Utilize  SPLASH  for  testing  upwind  and  downwind  conditions  at  yaw  angles  of  0,  3, 

and  5  degrees  with  rudder  angles  of  0,  3,  and  6  degrees. 

B.  Run  VPP  predictions  from  updated  SPLASH  appendage  data. 

C.  Reevaluate  the  high  pressure  and  large  vortex  regions  on  the  appendages.  Modify  the 

appendages. 

13.  Obtain  the  final  EHP  and  final  engine  selection  from  the  SPLASH  upright  predictions. 

14.  Run  VPP  predictions  from  finalized  SPLASH  sailing  data  until  satisfied  by  the  predicted 

aerodynamic  performance. 

A.  Modify  the  Mainsail  area,  Jib  Area,  and  Spinnaker  areas  to  maximize  boat  speed. 

B.  Assemble  matrix  of  mast  rake  and  trim  predictions. 

15.  Provide  the  resistance  data,  rig  and  trim  predictions,  and  righting  moment  curve  to  a 

sailmaker  for  evaluation  in  an  aerodynamic  CFD  program. 

16.  Generate  the  final  hull  lines  and  assemble  the  predicted  aerodynamic  forces  for  the  structural 

and  construction  design  phase. 

By  using  this  performance-based  design  process  in  the  design  spiral,  the  predicted 
performance  of  the  boat  is  maximized  by  varying  the  hull,  the  appendages,  and  the  sails. 
Through  utilizing  the  strengths  of  each  component  of  the  performance  design  process,  the 
predicted  error  is  minimized.  Furthermore,  because  the  experimental  elements  of  the  process 
have  been  shown  to  correlate  well  to  other  prediction  methods,  validation  tests  of  the 
components  are  not  required. 

Ultimately,  while  performance  prediction  remains  a  process  of  extrapolation  to  the  full- 
scale,  the  methods  used  in  prediction  will  continue  to  change  as  theories  improve  and  technology 
evolves. 
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Appendix  A  -  Navy  44s  Performance  Prediction  Comparison 

Comparison  to  other  boats  is  an  important  part  of  performance  prediction.  Since  the  Mk  II 
Navy  44  is  a  replacement  of  the  Mk  I  Navy  44,  the  comparison  of  the  two  designs  is  very 
relevant. 

The  basic  performance  prediction  measure  for  any  vessel  is  the  upright  resistance  results. 
The  Mk  I  Navy  44  was  tested  using  the  same  procedure,  setup  and  analysis  as  the  Mk  II. 


Velocity  (knots) 

Figure  A1  -Upright  testing  of  Mk  I  Navy  44  and  Mk  II  Navy  44 

Figure  (Al)  shows  that  the  Mk  II  Navy  44  has  reduced  resistance  for  speeds  over  seven 
knots.  This  improvement  is  due  mostly  to  the  improved  stem  shape  of  the  Mk  II  Navy  44.  The 
waterline  near  the  stem  of  the  Mk  II  Navy  44  is  beamier  than  the  Mk  I  Navy  44,  and  does  not 
include  a  bustle.  At  these  higher  velocities,  boats  suffer  from  an  effect  called  “squatting.”  The 
squatting  of  a  boat  is  the  stem  sinking  into  the  water,  making  it  more  difficult  to  enter  the  semi¬ 
planning  condition  of  high  speeds.  When  the  waterplane  area  of  the  stem  is  increased,  the  effect 
of  squatting  is  reduced,  and  the  boat  can  more  easily  increase  its  speed. 


no 

Sailing  performance  is  the  second  measure  of  performance  prediction.  Although  tank 
testing  sailing  tests  were  completed  with  the  Mk  I  Navy  44,  the  data  was  not  populated 
sufficiently  to  make  a  full  prediction  using  the  customized  VPP.  Therefore,  the  IMS  VPP 
prediction  was  used  exclusively  to  compare  the  sailing  performance  of  the  two  designs. 


Polar  prediction 

-  -  -eMki 

-  -  ■12Ukl 

Mk  i  Navy  44 

-  -  -ISMkl 

-  -  ■20Mkl 

and  the 

5F-1kll 

Mk  ii  Navy  44 

- 12Mkll 

16Mkll 

standard  sea  state  with  dynamic  allowance 

- 20Mkll 

0 


Figure  (A2)  shows  the  IMS  VPP 
Polars  of  the  predicted  speeds  of  the  Navy 
44s.  Unlike  the  earlier  polars  these  include 
the  added  resistance  due  to  waves  for  the 
appropriate  sea  states.  The  polars  showed 
two  important  design  changes  between  the 
boats. 

First,  the  Mk  II  Navy  44  showed  an 
average  speed  improvement  of  0.3  knots. 
Not  only  was  this  a  result  of  improved 
stem  waterplane  area,  but  was  also  a 
function  of  a  longer  waterline  and  a 
smoothed  underwater  profile.  The  longer 
waterline  of  the  Mk  II  Navy  44  acted  to 
decrease  the  Froude  number  of  the  hull, 
which  decreased  the  wave  making 
resistance.  The  smoothed  underwater 
profile  improved  the  flow  near  the  stem 
by  reducing  flow  separation. 

Finally,  the  upwind  sailing  angle  of 
the  Mk  II  Navy  44  was  improved  by 
almost  3  degrees.  This  increased  the 
ability  of  the  Mk  II  Navy  44  to  sail 
quickly  around  a  racecourse  or  to  reach  a 


Figure  A2  -Polars  of  Mk  I  Navy  44  and  Mk  II  Navy  44 


Ill 

destination  when  sailing  into  the  wind  at  sea.  The  VPP  did  not  include  rudder  design 
characteristics  however.  This  would  favor  the  newer  design  as  an  earlier  midshipman  project 
identified  the  Mk  1  rudder  as  having  significantly  more  drag  than  the  newer  design.*^  From  these 
studies,  it  appears  that  the  newer  design  should  perform  better  in  most  if  not  all  conditions. 


Ref:  Miller,  Paul  H.  “Student  Research  Projects  for  the  New  Navy  44  Sail  Training  Craft.”  16“^  Chesapeake 
Sailing  Yacht  Symposium.  (2003):  135. 
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Grid  Setups 

The  first  scripts  generated  the  code  to  be  run  in  gridgen.  The  first  code,  “mkglf  exec” 

generated  another  code  to  be  run  in  gridgen.  An  input  of  heel,  yaw,  rudder,  and  pitch  would 

create  the  necessary  gridgen  script.  A  file  named  “mkglf  data”  was  called  in  this  file.  This  is  a 

TCL/CL  script  previously  generated  which  would  be  added  to  the  script  at  runtime.  An  example 

of  the  output  of  this  script  is  the  last  code  included  in  this  section. 

==>  mkglf. exec  <== 

# ! /bin/sh 

name=H$lY$2R$3 . gif 
namegrd=H$lY$2R$3 . grd 

echo  $name 
touch  $name 

echo  "#  Gridgen  Journal  File  VI  (Gridgen  14.02  REL  1)"  >  $name 

echo  "gg::varSet  name  $namegrd"  >>  $name 

echo  "gg::varSet  heelalpha  -$1"  >>  $name 

echo  "gg::varSet  yawalpha  $2"  >>  $name 

echo  "gg::varSet  rudderalpha  -$3"  >>  $name 

echo  "gg::varSet  pitchalpha  0"  >>  $name 

cat  $name  mkglf. data  >  mkglf. tmp 

mv  mkglf. tmp  $name 


The  second  script  is  a  setup  of  the  sailing  condition  which  were  run  in  FKS.  This  script  uses 
the  previous  “mkglf  exec”  script  to  generate  the  files  for  use  in  Gridgen. 


==>  sail2.sh  <== 
# ! /bin/sh 
. /mkglf . exec 
. /mkglf . exec 
. /mkglf . exec 
. / mkglf . exec 
. / mkglf . exec 
. / mkglf . exec 
. /mkglf . exec 
. /mkglf . exec 
. / mkglf . exec 
. / mkglf . exec 
. / mkglf . exec 
. /mkglf . exec 
. /mkglf . exec 
. / mkglf . exec 
. / mkglf . exec 
. / mkglf . exec 
. /mkglf . exec 
. /mkglf . exec 
. / mkglf . exec 
. / mkglf . exec 


0  0  0 

2.5  0  0 

5  0  0 

7.5  0  0 

10  0  0 

12.5  0  0 

15  0  0 

17.5  0  0 

20  0  0 

22.5  0  0 

25  0  0 

27.5  0  0 

0  0  3 

2.5  0  3 

5  0  3 

7.5  0  3 

10  0  3 

12.5  0  3 

15  0  3 

17.5  0  3 


0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 
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/mkglf . exec 

20 

0 

3 

0 

/mkglf . exec 

22.5 

0 

3 

0 

/mkglf . exec 

25 

0 

3 

0 

/mkglf . exec 

27.5 

0 

3 

0 

/ mkglf . exec 

0 

0 

6 

0 

/ mkglf . exec 

2.5 

0 

6 

0 

/mkglf . exec 

5 

0 

6 

0 

/mkglf . exec 

7.5 

0 

6 

0 

/mkglf . exec 

10 

0 

6 

0 

/ mkglf . exec 

12.5 

0 

6 

0 

/ mkglf . exec 

15 

0 

6 

0 

/mkglf . exec 

17.5 

0 

6 

0 

/mkglf . exec 

20 

0 

6 

0 

/mkglf . exec 

22.5 

0 

6 

0 

/ mkglf . exec 

25 

0 

6 

0 

/ mkglf . exec 

27.5 

0 

6 

0 

/mkglf . exec 

0 

2 

0 

0 

/mkglf . exec 

2.5 

2 

0 

0 

/mkglf . exec 

5 

2 

0 

0 

/ mkglf . exec 

7.5 

2 

0 

0 

/ mkglf . exec 

10 

2 

0 

0 

/mkglf . exec 

12.5 

2 

0 

0 

/mkglf . exec 

15 

2 

0 

0 

/mkglf . exec 

17.5 

2 

0 

0 

/ mkglf . exec 

20 

2 

0 

0 

/ mkglf . exec 

22.5 

2 

0 

0 

/mkglf . exec 

25 

2 

0 

0 

/mkglf . exec 

27.5 

2 

0 

0 

/mkglf . exec 

0 

2 

3 

0 

/ mkglf . exec 

2.5 

2 

3 

0 

/ mkglf . exec 

5 

2 

3 

0 

/mkglf . exec 

7.5 

2 

3 

0 

/mkglf . exec 

10 

2 

3 

0 

/mkglf . exec 

12.5 

2 

3 

0 

/ mkglf . exec 

15 

2 

3 

0 

/ mkglf . exec 

17.5 

2 

3 

0 

/mkglf . exec 

20 

2 

3 

0 

/mkglf . exec 

22.5 

2 

3 

0 

/mkglf . exec 

25 

2 

3 

0 

/ mkglf . exec 

27.5 

2 

3 

0 

/ mkglf . exec 

0 

2 

6 

0 

/mkglf . exec 

2.5 

2 

6 

0 

/mkglf . exec 

5 

2 

6 

0 

/mkglf . exec 

7.5 

2 

6 

0 

/ mkglf . exec 

10 

2 

6 

0 

/ mkglf . exec 

12.5 

2 

6 

0 

/mkglf . exec 

15 

2 

6 

0 

/mkglf . exec 

17.5 

2 

6 

0 

/mkglf . exec 

20 

2 

6 

0 

/ mkglf . exec 

22.5 

2 

6 

0 

/ mkglf . exec 

25 

2 

6 

0 

/mkglf . exec 

27.5 

2 

6 

0 

/mkglf . exec 

0 

4 

0 

0 

/mkglf . exec 

2.5 

4 

0 

0 

/ mkglf . exec 

5 

4 

0 

0 

/ mkglf . exec 

7.5 

4 

0 

0 
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. /mkgif . exec 

10 

4 

0 

0 

. /mkgif . exec 

12.5 

4 

0 

0 

. /mkgif . exec 

15 

4 

0 

0 

. /mkgif . exec 

17.5 

4 

0 

0 

. / mkgif . exec 

20 

4 

0 

0 

. / mkgif . exec 

22.5 

4 

0 

0 

. /mkgif . exec 

25 

4 

0 

0 

. /mkgif . exec 

27.5 

4 

0 

0 

. /mkgif . exec 

0 

4 

3 

0 

. / mkgif . exec 

2.5 

4 

3 

0 

. / mkgif . exec 

5 

4 

3 

0 

. /mkgif . exec 

7.5 

4 

3 

0 

. /mkgif . exec 

10 

4 

3 

0 

. /mkgif . exec 

12.5 

4 

3 

0 

. / mkgif . exec 

15 

4 

3 

0 

. / mkgif . exec 

17.5 

4 

3 

0 

. /mkgif . exec 

20 

4 

3 

0 

. /mkgif . exec 

22.5 

4 

3 

0 

. /mkgif . exec 

25 

4 

3 

0 

. / mkgif . exec 

27.5 

4 

3 

0 

. / mkgif . exec 

0 

4 

6 

0 

. /mkgif . exec 

2.5 

4 

6 

0 

. /mkgif . exec 

5 

4 

6 

0 

. /mkgif . exec 

7.5 

4 

6 

0 

. / mkgif . exec 

10 

4 

6 

0 

. / mkgif . exec 

12.5 

4 

6 

0 

. /mkgif . exec 

15 

4 

6 

0 

. /mkgif . exec 

17.5 

4 

6 

0 

. /mkgif . exec 

20 

4 

6 

0 

. / mkgif . exec 

22.5 

4 

6 

0 

. / mkgif . exec 

25 

4 

6 

0 

. /mkgif . exec 

FKS  testinR 

27.5 

4 

6 

0 

The  first  step  in  running  FKS  was  to  create  the  control  files  necessary  for  input.  This  script 
generated  the  required  file  based  on  the  name  of  the  grid  to  be  used  and  the  Froude  number  to  be 
set.  The  script  also  ran  FKS  and  created  the  appropriate  output  files. 

==>  f 2. exec  <== 

# ! /bin/bash 

rm  -f  f ks . in 

rm  error. be 

cat  <<  EOF  >  f ks . in 

BASIC  INPUT  FILE  FOR  CODE  FKS 

1)  Enter  name  of  input  file  that  defines  the  ship  hull  surface 
grid . in 

2)  Enter  Froude  number 
EOF 

echo  $2  >>  f ks . in 
cat  <<EOF  >>  f ks . in 
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3)  Set  kwave  =  1  or  0  to  compute  or  ignore  wave  component 

1 

4)  Set  klocal  =  0,1, 2, 3  or  4  to  ignore, compute, inf  fluid, double  body,  or  inf 
speed  local  component 

1 

5)  Set  kdrag  =  1  or  0  to  compute  or  ignore  wave  drag 

1 

6)  Set  kext  equal  to  -1  ,  0  ,  1  ,  or  2 

kext  =  -1  :  compute  flow  along  ship  waterline 

kext  =  0  :  compute  flow  at  ship  hull  surface 

kext  =  1  :  compute  flow  at  free-surface 

kext  =  2  :  compute  flow  within  flow  domain 

0 

7)  If  kext  >  0  :  enter  name  of  input  file  defining  flow  field  points 
(  any  string  of  characters  is  acceptable  if  kext  =  0  or  -1  ) 

EOF 

echo  wavcutxyz_$l_$2 . dat  >>  f ks . in 
cat  <<EOF  >>  f ks . in 

8)  Enter  nrowref  =  number  of  longitudinal  rows  of  patch  ref  points 

2 

9)  For  i  between  1  &  nrowref  ,  enter 

yref  (i)  ,  zref  (i)  ,  xrefmin(i)  ,  xrefmax(i)  : 

0.0  ,  0.0  ,0.0  ,1 
0.0  ,  -0.2  ,  0.0  ,  1 

10)  Enter  name  of  basic  output  file  summarizing  main  io  information 
EOF 

echo  f kssum_$l_$2 . dat  >>  f ks . in 
cat  <<EOF  >>  f ks . in 

11)  If  kdrag  =  1  :  enter  name  of  output  file  to  plot  spectrum  function 
(  any  string  of  characters  is  acceptable  if  kdrag  =  0  ) 

EOF 

echo  spec_$l_$2 . dat  >>  f ks . in 
cat  <<EOF  >>  f ks . in 

12)  If  kext  >  -1  :  enter  name  of  output  file  to  PLOT  velocity  and 

pressure  computed  at  flow  field  points 
(  any  string  of  characters  is  acceptable  if  kext  =  -1  ) 

EOF 

echo  plotuvw_$l_$2 . dat  >>  f ks . in 
cat  <<EOF  >>  f ks . in 

13)  If  kext  =  -1  or  if  kext  =  0  with  nseg  >  0  :  enter  name  of  output 
file  to  PLOT  wave  profile  (  any  string  of  characters  is 
acceptable  if  kext  >  0  or  if  kext  =  0  with  nseg  =  0  ) 

EOF 

echo  wavprof_$l_$2 . dat  >>  f ks . in 
cat  <<EOF  >>  f ks . in 
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14)  If  kext  =  0  :  enter  name  of  output  file  giving  velocity  @  centroids 
of  hull  panels  &  waterline  segments  FOR  INPUT  TO  FKNK  &  FKX 

(  any  string  of  characters  is  acceptable  if  kext  not  equal  to  0  ) 

EOF 

echo  uvwpan_$l_$2 . dat  >>  f ks . in 
cat  <<EOF  >>  f ks . in 

15)  Enter  name  of  output  file  that  defines  patches/panels/segments 
EOF 

echo  pnseg_$l_$2 . dat  >>  f ks . in 
f ks . exec 
rm  -f  pan’^ .  dat 


Another  script  file  was  made  for  use  only  in  sailing  conditions.  For  multiple  hulls,  multiple 

pre-processings  had  to  be  made  using  the  program  “g2fks.”  This  program  pre-processed  each 

condition  and  then  ran  the  condition  using  the  above  code,  “fZ.exec.” 

==>  f.exec  <== 

# ! /bin/bash 
g2fks.exec  <<EOF 
2 

$1 . grd 

0 

0 

0 

0 

EOF 

hull=$l 

speed=$2 

./f2.exec  $hull  $speed 


The  last  two  codes  show  examples  of  how  a  set  of  upright  tests  and  sailing  tests  would  be 


performed. 


==>  Upright. sh  <== 
# ! /bin/bash 


/f 2 . exec 

HOYORO 

.10 

/f 2 . exec 

HOYORO 

.  12 

/f 2 . exec 

HOYORO 

.14 

/f 2 . exec 

HOYORO 

.16 

/f 2 . exec 

HOYORO 

.18 

/f 2 . exec 

HOYORO 

.20 

/f 2 . exec 

HOYORO 

.22 

/f 2 . exec 

HOYORO 

.24 

/f 2 . exec 

HOYORO 

.26 

/f 2 . exec 

HOYORO 

.28 

/f 2 . exec 

HOYORO 

.30 

/f 2 . exec 

HOYORO 

.32 

/f 2 . exec 

HOYORO 

.34 

/f 2 . exec 

HOYORO 

.36 

./f2.exec  HOYORO  .38 
./f2.exec  HOYORO  .40 
./f2.exec  HOYORO  .42 
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==>  sail2run.sh  <== 


.  /f . exec 
.It. exec 
.  /f . exec 
.  /f . exec 
.  /f . exec 
.  /f . exec 
.It.  exec 
. /f . exec 
. /f . exec 
.  /f . exec 
.  /f . exec 
.  /f . exec 
.  /f . exec 
. /f . exec 
.  /f . exec 
.  /f . exec 
.  /f . exec 
.  /f . exec 
. /f . exec 
.  /f . exec 
.  /f . exec 
.  /f . exec 
.  /f . exec 
. /f . exec 
.  /f . exec 
.  /f . exec 
.  /f . exec 
.  /f . exec 
. /f . exec 


HOYORO  0.1 
H2.5Y0R0  0.1 
H5Y0R0  0.1 
H7.5Y0R0  0.1 
HlOYORO  0.1 
H12.5Y0R0  0.1 
H15Y0R0  0.1 
H2.5Y0R3  0.1 
H5Y0R3  0.1 
H7.5Y0R3  0.1 
H15Y4R3  0.2 
H17.5Y4R3  0.2 
H20Y4R3  0.2 
H22.5Y4R3  0.2 
H25Y4R3  0.2 
H27.5Y4R3  0.2 
H0Y4R6  0.2 
H2.5Y2R0  0.3 
H5Y2R0  0.3 
H7.5Y2R0  0.3 
H10Y2R0  0.3 
H12.5Y2R0  0.3 
H15Y2R0  0.3 
H0Y4R0  0.3 
H2.5Y4R0  0.3 
H5Y4R0  0.3 
H7.5Y4R0  0.3 
H10Y4R0  0.3 
H15Y0R0  0.4 


etc 


Gridgen  Script 

The  last  code  is  the  TCL/CL  script  which  was  used  to  import  an  iges  surface  of  the  Navy  44 
into  Gridgen  and  get  out  a  nondimensional  plotSd  grid.  This  code  has  already  been  generated  by 
the  script  “mkglf  exec” 


gg::varSet  name  HOYORO. grd 
gg::varSet  heelalpha  -0 
gg::varSet  yawalpha  0 
gg::varSet  rudderalpha  -0 
gg::varSet  pitchalpha  0 
package  require  PWI_Glyph  1 . 4 
gg : : memClear 

gg::aswSet  GENERIC  -dim  3 
gg : : defReset 
gg : : tolReset 

gg: lupdatePolicy  DISPLAY  AND  INPUT 
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#Get  the  Mk  II  Navy  44  database  file 

gg :: db Import  "/home/c03/m036330/dimenNavy44mv. dba"  -type  DBA 
set  _DB(53)  [gg : : dbGetByName  --  {Navy44a-F2-2 } ] 

NOTE:  use  set  _DB(#)  to  all  iges  surface  names. 

The  rest  of  the  set  command  have  been  omitted. 

gg : : dbTransf ormBegin  [list  \ 

$_DB(53)  \ 

$_DB(55)  \ 

NOTE:  list  ail  DB  names  required  for  transformation. 

The  rest  of  the  DBs  have  been  omitted. 

]  -maintain_Iinkage 
gg : : xf ormTransIate  ".1  .1  .1" 
gg : : xf ormTransIate  "-.1  -.1  -.1" 
gg: : dbTransf ormEnd 

set  _DB(54)  [gg :: dbGetByName  --  {Navy44a-F3-3 } ] 

set  _DB(68)  [gg :: dbGetByName  --  {Navy44a-F3-3-copyl } ] 

set  _DB(4)  [gg :: dbGetByName  --  {Navy44a-E9-9 } ] 

set  _DB(5)  [gg :: dbGetByName  --  {Navy44a-E10-10 } ] 

set  _DB(6)  [gg :: dbGetByName  --  {Navy44a-Ell-ll } ] 

set  _DB(7)  [gg :: dbGetByName  --  {Navy44a-E12-12 } ] 

set  _DB(43)  [gg :: dbGetByName  --  {Navy44a-F3-21 } ] 

set  _DB(69)  [gg :: dbGetByName  --  {Navy44a-F3-21-copyl } ] 

set  _DB(70)  [gg :: dbGetByName  --  {Navy44a-E12-12-copyl } ] 

set  _DB(71)  [gg :: dbGetByName  --  {Navy44a-Ell-ll-copyl } ] 

set  _DB(72)  [gg :: dbGetByName  --  {Navy44a-E10-10-copyl } ] 

set  _DB(73)  [gg :: dbGetByName  --  {Navy44a-E9-9-copyl } ] 

#  apply  a  rudder  angle 

set  _ggTemp_(l)  $rudderalpha 
gg :: dbTransf ormBegin  [list  \ 

$_DB(54)  \ 

$_DB(68)  \ 

$_DB(4)  \ 

$_DB(5)  \ 

$_DB(6)  \ 

$_DB(7)  \ 

$_DB(43)  \ 

$_DB(69)  \ 

$_DB(70)  \ 

$_DB(71)  \ 

$_DB(72)  \ 

$_DB(73)  \ 

]  -maintain_Iinkage 

gg : : xf ormRotate  ".908104  0  0"  ".908104  0  -10"  $_ggTemp_(l) 
gg: : dbTransf ormEnd 
unset  _ggTemp_(l) 

#  apply  a  heel  angle 

set  _ggTemp_(l)  $heelalpha 
gg :: dbTransf ormBegin  [list  \ 

$_DB(53)  \ 

$_DB(54)  \ 

$  DB(55)  \ 
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NOTE:  list  all  DB  names  required  for  transformation. 

The  rest  of  the  DBs  have  been  omitted. 

]  -maintain_linkage 

gg : : xf ormRotate  [list  000]  [list  0.500024353  0  0]  $_ggTemp_(l) 
gg: : dbTransf ormEnd 
unset  _ggTemp_(l) 

#  apply  a  trim  angle 

set  _ggTemp_(l)  $pitchalpha 
gg :: dbTransf ormBegin  [list  \ 

$_DB(53)  \ 

$_DB(54)  \ 

$_DB(55)  \ 

NOTE:  list  all  DB  names  required  for  transformation. 

The  rest  of  the  DBs  have  been  omitted. 

]  -maintain_linkage 

gg :: xf ormRotate  ".5  0  0"  ".5  10  0"  $_ggTemp_(l) 
gg: : dbTransf ormEnd 
unset  _ggTemp_(l) 

#  apply  a  yaw  angle 

set  _ggTemp_(l)  $yawalpha 
gg :: dbTransf ormBegin  [list  \ 

$_DB(53)  \ 

$_DB(54)  \ 

$_DB(55)  \ 

NOTE:  list  all  DB  names  required  for  transformation. 

The  rest  of  the  DBs  have  been  omitted. 

]  -maintain_linkage 

gg :: xf ormRotate  ".5  0  0"  ".5  0  10"  $_ggTemp_(l) 
gg: : dbTransf ormEnd 
unset  _ggTemp_(l) 

set  _DB(142)  [gg::dbPlane  -z  [lindex  "0  0  0"  2]] 

#  Intersect  the  body  with  a  z  =  0  waterplane 
gg : : dblntersect  $_DB(142)  [list  \ 

$_DB(53)  \ 

$_DB(63)  \ 

$_DB(8  6)  \ 

$_DB(103)  \ 

]  -los  [list  001] 

set  _DB(144)  [gg : : dbGetByName  --  {Navy44a-Fl-211-isectl } ] 
set  _DB(143)  [gg :: dbGetByName  --  {Navy44a-F2-3-isectl } ] 
gg : : conBegin 

gg::segBegin  -type  DB_CUB1C 

gg : : segAddControlPt  -db  [list  1  0  $_DB(144)] 
gg: :segAddControlPt  -db  [list  0.0031643302  0  $_DB(144)] 
gg :: segAddControlPt  -db  [list  0.1243136332  0  $_DB(143)] 
gg :: segAddControlPt  -db  [list  1  0  $_DB(143)] 

gg : : segEnd 

set  _CN(1)  [gg::conEnd] 

set  _DB(146)  [gg :: dbGetByName  --  {Navy44a-Fl-211-copyl-isectl } ] 


set  _DB(145)  [gg : : dbGetByName  --  {Navy44a-F2-3-copyl-isectl } ] 
gg : : conBegin 

gg::segBegin  -type  DB_CUBIC 

gg : : segAddControlPt  -db  [list  1  0  $_DB(146)] 
gg: isegAddControlPt  -db  [list  0.006846838817  0  $_DB(146)] 
gg :: segAddControlPt  -db  [list  0.1173754558  0  $_DB(145)] 
gg :: segAddControlPt  -db  [list  1  0  $_DB(145)] 
gg : : segEnd 

set  _CN(2)  [gg::conEnd] 

set  _ggTemp_ ( 1 )  [gg : : conCetPt  $_CN(1)  -db  -arc  1] 
unset  _ggTemp_(l) 

set  _ggTemp_ ( 1 )  [gg : : conCetPt  $_CN(1)  -db  -arc  1] 
unset  _ggTemp_ ( 1 ) 
gg : : conBegin 

gg::segBegin  -type  DB_CUB1C 


gg 

: segAddControlPt 

-db 

[list 

1  0 

$ 

_DB  (143) ] 

gg 

: segAddControlPt 

-db 

[list 

.51 

0 

$ 

DB(3)  ] 

gg 

: segAddControlPt 

-db 

[list 

.  60 

0 

$ 

DB(3)  ] 

gg 

: segAddControlPt 

-db 

[list 

.  99 

0 

$ 

DB(3)  ] 

gg 

: segAddControlPt 

-db 

[list 

1  0 

$ 

_DB  ( 4 1 )  ] 

gg 

: segAddControlPt 

-db 

[list 

.7316 

0 

$_DB (41) ] 

gg 

: segAddControlPt 

-db 

[list 

0  0 

$. 

_DB  (36)  ] 

gg : : segEnd 

set  _CN(3)  [gg::conEnd] 
gg : : conBegin 

gg::segBegin  -type  DB_CUB1C 

gg :: segAddControlPt  -db  [list  1  0  $_DB(144)] 
gg :: segAddControlPt  -db  [list  0.2865071407  0  $_DB(41)] 
gg :: segAddControlPt  -db  [list  0.4740180149  0  $_DB(41)] 
gg :: segAddControlPt  -db  [list  1  0  $_DB(36)] 


gg : : segEnd 

set  _CN(4)  [gg::conEnd] 
gg : : conOnDBEnt  [list  \ 


$  DB ( 3  6 ) 

\ 

I 

$_DB (82) 

\ 

J 

gg: 

conDlm 

$ 

CN  (3) 

20 

set 

_CN ( 5 ) 

[ lindex 

[gg : : conCetAll 

gg: 

conDim 

$ 

_CN  ( 5 ) 

25 

gg: 

conDim 

$ 

_CN  ( 4 ) 

20 

gg: 

:  conDim 

$ 

_CN  ( 1 ) 

25 

gg: 

domBegin 

-type 

STRUCTURED 

gg : : edgeBegin 

gg : : edgeAddCon  $_CN(3) 


gg : : edgeEnd 
gg : : edgeBegin 

gg :: edgeAddCon  $_CN(5) 
gg : : edgeEnd 
gg : : edgeBegin 

gg : : edgeAddCon  $_CN ( 4 ) 
gg : : edgeEnd 
gg : : edgeBegin 

gg :: edgeAddCon  $_CN(1) 
gg : : edgeEnd 

set  _DM(1)  [gg : : domEnd] 
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gg::conDim  $_CN(2)  25 

set  _CN(6)  [lindex  [gg : : conGetAll ]  5] 

gg::conDim  $_CN(6)  25 
gg : : domBegin  -type  STRUCTURED 
gg : : edgeBegin 

gg : : edgeAddCon  $_CN ( 2 ) 
gg : : edgeEnd 
gg : : edgeBegin 

gg : : edgeAddCon  $_CN ( 4 ) 
gg : : edgeEnd 
gg : : edgeBegin 

gg :: edgeAddCon  $_CN(6) 
gg : : edgeEnd 
gg : : edgeBegin 

gg :: edgeAddCon  $_CN(3) 
gg : : edgeEnd 

set  _DM(2)  [gg : : domEnd] 

set  ~ggTemp_(l)  [list  $_DM(1)  $_DM(2)] 

gg: : domEllSolverBegin  $_ggTemp_ ( 1 ) 
gg : : domTFISolverRun  $_ggTemp_ ( 1 ) 
gg : : domEllSolverStep  -iterations  4  -nodisplay 
gg: : domEllSolverEnd 
unset  _ggTemp_ ( 1 ) 
gg : : conOnDBEnt  $_DB(35) 
gg : : conOnDBEnt  $_DB(8) 
gg: : conOnDBEnt  $_DB(37) 
gg :: conOnDBEnt  $_DB(83) 
gg: : conOnDBEnt  $~DB(78) 
gg: : conOnDBEnt  $_DB(81) 

set  _CN(12)  [lindex  [gg :: conGetAll ]  11] 

gg::conDim  $_CN(12)  25 

set  _CN(11)  [lindex  [gg :: conGetAll ]  10] 

gg::conDim  $_CN(11)  25 

set  _CN(10)  [lindex  [gg :: conGetAll ]  9] 

gg::conDim  $_CN(10)  25 
gg :: domBegin  -type  STRUCTURED 
gg : : edgeBegin 

gg :: edgeAddCon  $_CN(6) 
gg : : edgeEnd 
gg : : edgeBegin 

gg : : edgeAddCon  $_CN (12) 
gg : : edgeEnd 
gg : : edgeBegin 

gg :: edgeAddCon  $_CN(11) 
gg : : edgeEnd 
gg : : edgeBegin 

gg :: edgeAddCon  $_CN(10) 
gg : : edgeEnd 

set  _DM(3)  [gg::domEnd] 

set  _CN(7)  [lindex  [gg :: conGetAll ]  6] 

gg::conDim  $_CN(7)  25 

set  _CN(8)  [lindex  [gg :: conGetAll ]  7] 

gg::conDim  $_CN(8)  25 

set  _CN(9)  [lindex  [gg :: conGetAll ]  8] 

gg::conDim  $_CN(9)  25 

gg:: domBegin  -type  STRUCTURED 
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gg : : edgeBegin 

gg : : edgeAddCon  $_CN(7) 
gg : : edgeEnd 
gg : : edgeBegin 

gg : : edgeAddCon  $_CN ( 8 ) 
gg : : edgeEnd 
gg : : edgeBegin 

gg :: edgeAddCon  $_CN(9) 
gg : : edgeEnd 
gg : : edgeBegin 

gg : : edgeAddCon  $_CN ( 5 ) 
gg : : edgeEnd 

set  _DM(4)  [gg : : domEnd] 

gg : : dispViewDoms  FALSE 

gg : : dbintersect  $_DB(142)  [  list  \ 

$_DB(54)  \ 

$_DB(68)  \ 

]  -los  [list  001] 

gg : : dispViewCenter  [gg : : dbUVToXYZ  [gg : : conGetPt  $_CN(1)  -db  -arc  0]] 


set 

_DB (147) 

[gg: 

:  dbCetByName  -- 

{Navy4  4a-F3-21-isectl }  ] 

set 

DB ( 1 4  8 ) 

[gg: 

:  dbCetByName  -- 

{Navy4  4a-F3-2 1-copyl-isectl 

gg:  : 

:  conBegin 

gg::segBegin  -type  DB_CUBIC 

gg : : segAddControlPt  -db  [list  0  0  $_DB(4)] 
gg : : segAddControlPt  -db  [list  0.66  0  $_DB(4)] 
gg :: segAddControlPt  -db  [list  1  0  $_DB(147)] 
gg : : segEnd 

set  _CN(13)  [gg::conEnd] 
gg : : conBegin 

gg::segBegin  -type  DB_CUBIC 

gg :: segAddControlPt  -db  [list  0  0  $_DB(73)] 
gg :: segAddControlPt  -db  [list  0.66  0  $_DB(73)] 
gg :: segAddControlPt  -db  [list  1  0  $_DB(148)] 
gg : : segEnd 

set  _CN(14)  [gg::conEnd] 
gg: : conOnDBEnt  $_DB(70) 
gg: : conOnDBEnt  $_DB(71) 
gg : : conBegin 

gg::segBegin  -type  DB_CUBIC 

gg :: segAddControlPt  -db  [list  1  0  $_DB(72)] 
gg :: segAddControlPt  -db  [list  0.1  0  $_DB(72)] 
gg :: segAddControlPt  -db  [list  1  0  $_DB(147)] 
gg : : segEnd 

set  _CN(17)  [gg::conEnd] 
gg :: conOnDBEnt  $_DB(7) 
gg :: conOnDBEnt  $_DB(6) 
gg : : conBegin 

gg::segBegin  -type  DB_CUBIC 

gg :: segAddControlPt  -db  [list  1  0  $_DB(5)] 
gg :: segAddControlPt  -db  [list  0.1  0  $_DB(5)] 
gg :: segAddControlPt  -db  [list  1  0  $_DB(148)] 
gg : : segEnd 

set  _CN(20)  [gg::conEnd] 
gg::conDim  $_CN(14)  10 
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gg::conDim  $_CN(20)  15 

set  _CN(19)  [lindex  [gg : : conGetAll ]  15] 

gg::conDim  $_CN(19)  10 

set  _CN(18)  [lindex  [gg :: conGetAll ]  14] 

gg::conDim  $_CN(18)  15 

gg : : domBegin  -type  STRUCTURED 
gg : : edgeBegin 

gg : : edgeAddCon  $_CN (14) 
gg : : edgeEnd 
gg : : edgeBegin 

gg : : edgeAddCon  $_CN (20) 
gg : : edgeEnd 
gg : : edgeBegin 

gg :: edgeAddCon  $_CN(19) 
gg : : edgeEnd 
gg : : edgeBegin 

gg :: edgeAddCon  $_CN(18) 
gg : : edgeEnd 

set  _DM(5)  [gg : : domEnd] 

set  _CN(15)  [lindex  [gg :: conGetAll ]  11] 

gg::conDim  $_CN(15)  15 

set  _CN(16)  [lindex  [gg :: conGetAll ]  12] 


gg: 

:  conDim 

$ 

CN (16) 

10 

gg: 

:  conDim 

$ 

CN (17) 

15 

gg: 

:  conDim 

$ 

_CN (13) 

10 

gg: 

: domBegin 

-type 

STRUCTURED 

gg : : edgeBegin 

gg :: edgeAddCon  $_CN(15) 
gg : : edgeEnd 
gg : : edgeBegin 

gg : : edgeAddCon  $_CN (16) 
gg : : edgeEnd 
gg : : edgeBegin 

gg :: edgeAddCon  $_CN(17) 
gg : : edgeEnd 
gg : : edgeBegin 

gg :: edgeAddCon  $_CN(13) 
gg : : edgeEnd 

set  _DM(6)  [gg::domEnd] 
gg : : dispViewDoms  TRUE 
gg: : dispViewReset 
set  _ggTemp_ ( 1 )  $name 

set  _ggTemp_ (2 )  [file  join  $_ggTemp_ ( 1 ) ] 

unset  _ggTemp_ ( 1 ) 

#  Output  as  PlotSd  format  in  local  directory 
gg : : domExport  ALL  $_ggTemp_(2)  \ 

-style  PLOTSD  -form  ASCII  -precision  SINGLE 
unset  _ggTemp_(2) 


Appendix  C  -  Data  Analysis  Codes 
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MATLAB  Script 

This  MATLAB  script  was  a  generic  script  which  would  input  a  data  file  and  tranform  a 
sprarse  grid  into  a  dense,  square  grid.  This  particular  script  was  used  on  the  tank  testing  data. 
Other  scripts  were  modified  from  this  script  to  analyze  FKS  data  or  SPLASH  data. 


==>  gridsparsedata . m  <== 

for  yc  =  0:1 

for  rc  =  0:2 

for  wc  =  3:5 

yaw=yc’^4  ; 

rudder=rc’^3 ; 

whichone=wc;  %whichone  refers  to  drag,  lift  or  yaw  in  columns  3, 4, or  5 

D=load ( ' data . txt ' ) ;  %data  is  loaded  from  the  file  data.txt 

O=load ( ' initzero . txt ' ) ; 

start=l+yaw’^5  +  rudder'^4  0/  3  ; 

if  whichone==3 

type= ' Drag ' ; 

AZ=-37 . 5; 

elseif  whichone==4 

type= ' Forward  Lift'; 

AZ=60; 

elseif  whichone==5 

type='Aft  Lift'; 

AZ=-170; 

end 

%heel  is  in  column  1,  which  is  the  y  value 
%speed  is  in  column  2,  which  is  the  x  value 

x= [0 ( : , 2) ; D (start : start+19, 2) ] ;  y= [0 ( : , 1) ; D (start : start+19, 1)]; 
z= [0 ( : , whichone) ;D(start:start+19, whichone) ] ; 
xm  =  max (x) ;  ym=max(y); 
x=x/xm;  y=y/ym; 
ti  =  0: .05:1; 

[XI, YI]  =  meshgrid ( ti, ti) ; 

ZI  =  griddata (x, y, z , XI , YI ,' cubic ') ; 

YI  =  YI  *ym; 
y=  y'^'ym; 

XI  =  XI  ’^xm; 

X  =  X  xm  ; 

%mesh (XI, YI, ZI) ,  hold 
%pIot3 (x, y, z , ' o ' ) ,  hold  off 

%ylabel ( ' Heel  Angle  [degrees ]'), xiabel (' Speed  [ ft/s ]'), zlabel ([ type  ' 

[lbs] ']); 

%title([type  '  with  '  int2str(yaw)  '  degrees  of  yaw  and  '  int2str (rudder )  ' 

degrees  of  rudder']); 

%view (AZ , 45 ) 

%  saveas (1 , [ type  ''  int2str(yaw)  ''  int2str (rudder )  ' . bmp ' ] , ' bmp ' ) 

%  save  type  XI  YI  ZI  -ASCII  -TABS 
%true  velocity  is  in  column  8 
%true  heel  is  in  column  9 
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w  =  interp2 (XI, YI, ZI, D (start : start+19, 8) , D (start : start+19, 9) , ' cubic ' ) 

%name  =  [type  ''  int2str(yaw)  ''  int2str (rudder )  '.interp'] 

%save (name, ' w' , ' -ASCII ' , ' -TABS ' ) 
if  whichone==3 
wl  =  w 

elseif  whichone==4 
w2  =  w 
ZIf  =  ZI 
zf  =  z 

elseif  whichone==5 
w3  =  w 
Zla  =  ZI 
za  =  z 

end 

end 

ZIn  =  ZIf  +  Zla 

zn  =  zf  +  za 

mesh (XI, YI, ZIn) ,  hold 

contour 3 (XI , YI , ZIn, 15 ) 

plot3 (x, y, zn, ' o ' ) ,  hold  off 

ylabel('Heel  Angle  [degrees ]'), xlabel (' Speed  [ ft/s ]'), zlabel ([' Side  Force 
[lbs] ']); 

title (['Side  force  with  '  int2str(yaw)  '  degrees  of  yaw  and  '  int2str (rudder ) 
'  degrees  of  rudder']); 
view (AZ , 45 ) 

saveas ( 1 , [ ' SF '  int2str(yaw)  ''  int2str (rudder )  ' . bmp ' ] , ' bmp ' ) 

ww  =  [wl,w2,w3] 

name  =  ['i'  int2str(yaw)  ''  int2str (rudder )  '.interp'] 

save (name, ' ww ' , ' -ASCII ' , ' -TABS ' ) 

end 

end 


Maple  6  Analysis 

Below  is  the  script  used  to  develop  lift  and  drag  for  any  yaw  angle  and  rudder  angle. 

>  restart; 

>L:  =  (AH5,AH3)  - 

>(alphal*AH5+alpha2) *AH3^2+ (betal*AH5+beta2) *AH3+ (gammal*AH5+gam 
ma2)  ; 

L  :=  (AH5,  AH3  )  ^  ( al  AH5  +  al )  AH3^  +  ( p  1  AH5  +  p2  )  AH3  +  yl  AH5  +  y2 

>  L ( 0 , 0 ) =AG2 ;  assign ( % ) ; 

y2=AG2 

>  L  ( 4 , 0 )  =AG3 ;  gainmal=solve  ( % ,  gammal )  ;  assign  ( % )  ; 

dyl  +AG2=AG3 

yl  =-\aG2  +\aG3 
4  4 

>L(0,3)=AG4;  alpha2=solve  (%  ,alpha2)  ;  assign(%); 
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9a2  +  3  |32  +AG2=AG4 
a2  =  -1^2 -  I AG2  +IaG4 

3  y  y 

>L(4,3)=AG5;  alphal=solve  (%  ,alphal)  ;  assign(%); 

36  al  -AG2  +AG4  +  12  pi  +AG3  =AG5 

>L(0,6)=AG6;  beta2=solve  (%  ,beta2)  ;  assign(%); 

-6  P2  -  3  AG2  +  4  AG4  =AG6 

Ii2  =-IaG2  +  ^AG4-^AG6 
2  3  6 

> L (4 , 6) =AG7 ;  betal=solve (% ,betal) ;  assign(%); 

-3  AGS  +  3  AG2  -  4  AG4  -  24  pi  +  4  AG5  +AG6  =AG7 

pi  =-IaG3+IaG2-^AG4+^AG5+^AG6-^AG7 
^  8  8  6  6  24  24 


>L(AH5,AH3)  ; 


1 


1 


1 


1 


1 


1 


1 


—  AG3-  —  AG2+  —  AG4-  —  AG5-  —  AG6  +  —AG7  \AH5+  —  AG2--AG4 


72 


72 


1 


36 

/ 


36 


72 


72 


18 


+  ^AG6\aH3^  + 


-\aG3+\aG2-\aG4  +  \aG5  +  ^AG6-^AG7\aH5 
o  o  o  6  24  24  i 


\aG2+^AG4-\aG6  \AH3  + 
2  3  6 


\aG2+\aG3^AH5  +AG2 
4  4  j 


>R:  =  (AH5,AH3)  - 

>subs (AG6=AD6  r  subs (AG5=AD5 , subs (AG4=AD4 , subs (AG3=AD3 , subs (AG2=AD 
2 , subs (AG7=AD7 ,L (AH5 ,AH3) )))))); 

R  :=  (AH5,  AH3  )  s\ihs{AG6  =  AD6,  subs(^G5  =AD5,  subs(^G4  =AD4, 

subs(^Gi  =AD3,  subs(^G2  =AD2,  subs(^G7  =^Z)7,  GiAHS ,  AH3))))))) 


>R(AH5,AH3)  :=R(AH5,AH3)  ; 

K{AH5,  AH3):=((^AD3-^AD2+j^AD4-j^  ADS  -^^06  +  ^  ADt]  AH5 


1 


1 


1 


^-AD2--AD4^-AD6 


AH3^  + 


1 


1 


1 


1 


1 


1 


1 


-  -  ADS  +  -  AD2  -  -  AD4  +  -  ADS  +  —  AD6  -  —  AD7\  AHS  -  -  AD2  +  -  AD4 


8 


24 


24' 


1 


1 


-AD6  \AH3  +  \--AD2  +  -AD3  \AHS+AD2 


>  xeqnL  :=solve  (L  (AHS  ,AH3)  =AY2  ,AH5)  ; 
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xeqnL  :=  4  (-AH3^  AG2  +  2  AH3^  AG4  -  AH3^  AG6  -  1 8  ^G2  +  9  AH3  AG2 

-  12  AH3  AG4  +  3  AH3  AG6  +  18  ^72)  /  (-9  AH3  AG3  +  9  AH3  AG2 

-  12  AH3  AG4  +  12  AH3  AGS  +  3  AH3  AG6  -  3  AH3  AG7  -  18  AG2  +  18  AG3 

+  AH3^AG3-AH3^AG2  +  2AH3^AG4-2AH3^AG5-AH3^AG6  +  AH3^AG7) 


> Rf y : =subs ( AH5=xeqnL , R ( AH5 , AH3 ) ) ; 

Rfy  ■=  (^4  {^AD3  -^AD2  +j^AD4  -j^ADS  -^AD6  +  ^AD7']^  {-AH3^  AG2 
+  2  AH3^AG4  -AH3^  AG6  -  18  ^G2  +  9  AH3  AG2  -  12  AH3  AG4  +  3  AH3  AG6 
+  \SAY2)^  /  (-9  AH3  AG3  +  9  AH3  AG2  -  12  AH3  AG4  +  12  AH3  AGS 
+  3  AH3  AG6  -  3  AH3  AG7  -  18  AG2  +  18  AGS  +AH3^AG3  -AH3^AG2 


1 


1 


+  2  AH3^AG4  -  2  AH3^  AGS  -  AH3^AG6+AH3^  AG7)  +jwAD2  -^AD4 

18  9 

+  ^  ads]  AH3^  +  f4  f  f- 1  AD  3  +  \  AD2  - 1 AD4  +  \  ADS  +  ^  ADS  -^AD7\ 


18 


8 


8 


24' 


24 


-AH3^AG2  +  2  AH3^  AG4  -AH3^  AGS  -  18  AG2  +  9  AH3  AG2  -  12  AH3  AG4 

+  3  AH3  AGS  +  18  yl72)1  /  (-9  AH3  AGS  +  9  AH3  AG2  -  12  AH3  AG4 

J 

+  nAH3AGS  +  3AH3AGS-3AH3AG7-\SAG2  +  \SAG3+AH3^AG3 
-AH3^AG2  +  2  AH3^  AG4  -  2  AH3^  AGS  -AH3^  AGS +  AH3^  AG7)  -^AD2 

+ 1 AD4  - 1  ADs')^  AH3  +  4  ^  AD2  + 1  ADs')^  ( -AH3^  AG2 +2  AH3^  AG4 

-AH3^  AGS  -  18  ylG2  +  9  AH3  AG2  -  12  AH3  AG4  +  3  AH3  AGS  +  18  ^72)  /  ( 
-9  AH3  AGS  +  9  AH3  AG2  -  12  AH3  AG4  +  12  AH3  AGS  +  3  AH3  AGS 

-  3  AH3  AG7  -  18  ylG2  +  18  AGS  +AH3^AG3  -AHS^AG2  +  2  AHS^  AG4 
-2AH3^AGS-AHS^AGS+AH3^AG7)+AD2 


> isolate (Rfy=res ,AH3) ; 

(-54  AD7AHS  AG2  +  9  AD7  AHS^  AGS  -  36  AD7  AHS^  AG4  +  45  AD7  AHS^  AG2 
+  54  AD2  AHS  AGS  +  117  ADS  AHS^  AG2  -  324  ADS  AHS  AG2 
+  45  ADS  AHS^AGS  -  144  ADS  AHS^  AG4  -  54  ADS  AHS  AGS 

-  1 8  AD2  AHS^  AG2  +  36  AD2  AHS^  AG4  +  54  AHS  AD7  AY2  +  162  AD2  AHS  AG2 

-  18  AD2  AHS^  AGS  -  54  AHS  ADS  AY2  -  216  AHS  ADS  AY2  +216  AHS  AD4  AY2 

-  162  AHS  AD2  AY2  +  162  AHS  ADS  AY2  +  216  ADS  AHS  AG2 

-  36  ADS  AHS^  AGS  +  144  ADS  AHS^  AG4  -  144  ADS  AHS^  AG2 

-  216  AD2  AHS  AG4  -  9  ADS  AHS^  AG7  +  36  ADS  AHS^  AGS  -  45  ADS  AHS^  AGS 
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+  54  AD6AH3  AG3  +  36  AD4  AH3^  AG7  -  144  AD4  AH3^  AGS 
+  144  AD4  AH3^  AG3  -  216  AD4  AH3  AG3  -  27  AD2  AH3^  AG7 
+  108  AD2  AH3^  AGS  -  99  AD2  AH3^  AG3  +  162  AD2  AH3  AGS 
+  6  AH3 ^  AD6  ^ G7  -  1 8  AH 3 ^  AD6  AGS +\2  AH 3 ^  AD6  AGS +  2  AH3 ^  AD4  AG7 

-  4  AH3^  AD4  AGS  +  2  AH3^  AD4  AGS  -  18  AHS^  AD4  AG7  +  48  AHS^  AD4  AGS 

-  30  AHS'^  AD4  AGS-  AHS^  AD2  AG7 +  2  AH3^  AD2  AGS-  AH3^  AD2  AGS 
+  12  AHS^  AD2  AG7  -  30  AHS^  AD2  AGS  +  18  AHS^  AD2  AGS 

+  18  AHS^  AD6AY2  +  36  AHS^  ADS  AY2 -36  AHS^  AD4  AY2  +  \S  AHS^  AD2  AY2 

-  IS  AHS^  ADS  AY2  -  \2  AHS^  ADS  AG6  -  18  AHS^  ADS  AG2  +  AHS^  ADS  AG6 

-  2  AHS^  ADS  AG4+  AHS^  ADS  AG2  -  \S  AHS^  AD7AY2-6  AH3^  AD7  AG6 
+  \SAHS^  AD7AG4-\2AHS^AD7  AG2+AHS^AD7  AG6 -2  AH3^AD7AG4 
+  AHS^AD7AG2  +  18  AH3^  ADS  AG6  +  30  AH3^  ADS  AG4  -  48  AH3^  ADS  AG4 
+  30  AH3^  ADS  AG2 -2  AHS^  ADS  AG6  +  4  AHS^  ADS  AG4-2  AHS^  ADS  AG2 

-  AHS^  AD6  AG7 +  2  AH 3^  AD6  AGS-  AH 3^  AD6  AGS +2\6  ADS  AH 3  A  G4 

+  324  ADS  AG2  -  324  ADS  AY2  -  324  AD2  AG2  +  324  AD2  AY2)  /  (-9  AHS  AGS 
+  9  AHS  AG2  -  12  AHS  AG4  +  12  AHS  AGS  +  3  AHS  AG6  -  3  AHS  AG7  -  18  ylG2 
+  18  +AH3^AG3  -AH3^AG2  +  2AH3^AG4  -  2  AHS^  AGS  -  AHS^  AG6 
+  AHS^  AG7)  = -IS  res +18  AD2 


Before  data  analysis  began,  Maple  was  also  used  to  determine  the  rudder  angle  from  the 
reading  of  the  caliper  in  tank  testing. 


>  restart: 

>  el :  =  (C3) ^2=A^2+C^2-2*A*C*cos (beta) ; 

el  :=C3^=A^  +  C^-2ACcos(P) 


> betarad:=solve (el ,beta) ; 


f 

betarad  :=n-  arccos 

V 


1  CS^-A^-C 

2  AC 


7 


> A: =10. 31;  C:=10; 


A  :=  10.31 


C:=10 


>  betarad; 

71  -  arccos(  .004849660524  CS^  -  1.000466052  ) 

> beta :=evalf (betarad*180/Pi) ; 

p  :=  180.0000000  -57.29577950  arccos(  .004849660524  -  1 .000466052  ) 


>  beta; 
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180.0000000  -57.29577950  arccos(  .004849660524  0^-  1.000466052) 


>  solve (beta=ang,C3) ; 

12124?5131  ^3031037827500000000  + 


10 

1212415131 


303245044874 1582030  +3031037827500000000  cos| 


360000000 
114591559  ^ 


>  e2 :=angle=180-57 . 2957795*arccos (0 . 004849660524* (C3+2) ^2- 
1.000466052) -11.16683841; 

e2  :=  angle  = 

168.8331616  -57.2957795  arccos(. 004849660524  (0+2)^-  1.000466052) 


>b:=solve  (e2  ,C3)  [1]  ; 
10 


b\=-2  + 


1212415131 


3032450448741582030  +3031037827500000000  cos 


200000 

1145915 


>  evalf (b, 4) ; 

-2.  +  .8248  10-*  V-3032  10 ‘‘^  +  .3031  10**^  cos(. 01 745329252  angle -2. 94669 A551 ) 


Appendix  D  -  Velocity  Prediction  Program 
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The  first  visual  basic  module  was  included  in  the  Mk  11  Navy  44  VPP.  This  module 
interpolated  through  the  data,  and  also  solved  for  the  righting  moment. 

Attribute  VB_Name  =  "Modulel" 

Function  rarm(heel) 

rarm  =  heel  ^  3  0.01831  +  heel  ^  2  5.249  +  heel  -2355.7 

End  Function 

Function  shull (speed,  rud,  yaw,  heel) 

Application . Volatile 
Const  dOO  =  1010 

Const  d03  =  dOO  +  24 

Const  d06  =  d03  +  24 

Const  d40  =  d06  +  24 

Const  d43  =  d40  +  24 

Const  d46  =  d43  +  24 

D1  =  interp(d00,  speed,  heel) 

D2  =  interp(d03,  speed,  heel) 

D3  =  interp(d06,  speed,  heel) 

D4  =  interp(d40,  speed,  heel) 

D5  =  interp(d03,  speed,  heel) 

D6  =  interp(d46,  speed,  heel) 

Shull  =  1  /  72  rud  ^  2  yaw  D4  -  1  /  72  rud  ^  2  yaw  D1  +  1  /  36 
rud  2  yaw  ’^D2-l/36’^  rud  2  yaw  ’^D5-l/72’^  rud  2  yaw 
D3  +  1  /  72  rud  ^  2  yaw  D6  +  1  /  18  rud  ^2’'Dl-l/9’'  rud  ^  2 

D2  +  1  /  18  rud  '^2’^D3-1/  8^^  rud  yaw  D4  +  1  /  8  rud  yaw  D1 

-  1  /  6  rud  yaw  D2  +  1  /  6  rud  yaw  ’^D5  +  l/  24’^  rud  yaw  D3  - 
1  /  24  rud  yaw  D6  -  1  /  2  rud  D1  +  2  /  3  rud  D2  -  1  /  6  rud 
D3  -  1  /  4  yaw  D1  +  1  /  4  yaw  ^^04  +  01 

End  Function 


Function  liftk (speed,  rud,  yaw,  heel) 

Application .Volatile 
Const  dOO  =  290 

Const  d03  =  dOO  +  24 

Const  d06  =  d03  +  24 

Const  d40  =  d06  +  24 

Const  d43  =  d40  +  24 

Const  d46  =  d43  +  24 

D1  =  interp(d00,  speed,  heel) 

D2  =  interp(d03,  speed,  heel) 

D3  =  interp(d06,  speed,  heel) 

D4  =  interp(d40,  speed,  heel) 

D5  =  interp(d43,  speed,  heel) 

D6  =  interp(d46,  speed,  heel) 

liftk  =  1  /  72  rud  ^  2  yaw  D4  -  1  /  72  rud  ^  2  yaw  D1  +  1  /  36 
rud  2  yaw  ’^02-1/36’^  rud  2  yaw  ’^05-1/72’^  rud  2  yaw 

D3  +  1  /  72  rud  ^  2  yaw  D6  +  1  /  18  rud  "'2’'Dl-l/9’'  rud  ^  2 

D2  +  1  /  18  rud  '^2’^D3-1/  8^^  rud  yaw  D4  +  1  /  8  rud  yaw  D1 

-  1  /  6  rud  yaw  D2  +  1  /  6  rud  yaw  ’^D5  +  l/  24’^  rud  yaw  ^^03- 
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1  /  24  ^  rud  ^  yaw  ^  D6  -  1  /  2  ^  rud  ^  D1  +  2  /  3  ^  rud  ^  D2  -  1  /  6  ^  rud  ^ 

D3  -  1  /  4  ^  yaw  D1  +  1  /  4  ^  yaw  ^^04  +  01 

End  Function 

Function  liftr (speed,  rud,  yaw,  heel) 

Application . Volatile 
Const  dOO  =  1586 

Const  d03  =  dOO  +  24 

Const  d06  =  d03  +  24 

Const  d40  =  d06  +  24 

Const  d43  =  d40  +  24 

Const  d46  =  d43  +  24 

D1  =  interp(d00,  speed,  heel) 

D2  =  interp(d03,  speed,  heel) 

D3  =  interp(d06,  speed,  heel) 

D4  =  interp(d40,  speed,  heel) 

D5  =  interp(d43,  speed,  heel) 

D6  =  interp(d46,  speed,  heel) 

liftr  =  1  /  72  rud  ^  2  yaw  D4  -  1  /  72  rud  ^  2  yaw  D1  +  1  /  36 

rud  ^  2  yaw  ’^D2-l/36’^  rud  ^  2  yaw  ’^D5-l/72’^  rud  ^  2  yaw 

D3  +  1  /  72  rud  2  yaw  D6  +  1  /  18  rud  '^2’^Dl-l/9’^  rud  2 

D2  +  1  /  18  rud  "^2'*'D3-l/8’^  rud  yaw  D4  +  1  /  8  rud  yaw  D1 

-  1  /  6  rud  yaw  D2  +  1  /  6  rud  yaw  ’^D5  +  l/  24’^  rud  yaw  D3  - 

1  /  24  rud  yaw  D6  -  1  /  2  rud  D1  +  2  /  3  rud  D2  -  1  /  6  rud 

D3  -  1  /  4  yaw  D1  +  1  /  4  yaw  ^^04  +  01 

End  Function 

Function  drag (speed,  rud,  yaw,  heel) 

Application . Volatile 
Const  dOO  =  1154 

Const  d03  =  dOO  +  24 

Const  d06  =  d03  +  24 

Const  d40  =  d06  +  24 

Const  d43  =  d40  +  24 

Const  d46  =  d43  +  24 

D1  =  interp(d00,  speed,  heel) 

D2  =  interp(d03,  speed,  heel) 

D3  =  interp(d06,  speed,  heel) 

D4  =  interp(d40,  speed,  heel) 

D5  =  interp(d43,  speed,  heel) 

D6  =  interp(d46,  speed,  heel) 

drag  =  1  /  72  rud  ^  2  yaw  ’^D4-l/72’^  rud  ^  2  yaw  ’^Dl  +  1/  36’^ 

rud  2  yaw  ’^02-1/36’^  rud  2  yaw  ’^05-1/72’^  rud  2  yaw 

D3  +  1  /  72  rud  ^  2  yaw  D6  +  1  /  18  rud  ^2’'Dl-l/9’'  rud  ^  2 

D2  +  1  /  18  rud  '^2’^D3-l/8’^  rud  yaw  D4  +  1  /  8  rud  yaw  D1 

-  1  /  6  rud  yaw  D2  +  1  /  6  rud  yaw  ’^D5  +  l/  24’^  rud  yaw  D3  - 

1  /  24  rud  yaw  D6  -  1  /  2  rud  D1  +  2  /  3  rud  D2  -  1  /  6  rud 

D3  -  1  /  4  yaw  D1  +  1  /  4  yaw  "^04  +  01 

If  drag  <  0  Then 
drag  =  0 
End  If 

End  Function 


Function  dragw (speed,  rud,  yaw,  heel) 
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Application . Volatile 
Const  dOO  =  1442 
Const  d03  =  dOO  +  24 

Const  d06  =  d03  +  24 

Const  d40  =  d06  +  24 

Const  d43  =  d40  +  24 

Const  d46  =  d43  +  24 

D1  =  interp(d00,  speed,  heel) 

D2  =  interp(d03,  speed,  heel) 

D3  =  interp(d06,  speed,  heel) 

D4  =  interp(d40,  speed,  heel) 

D5  =  interp(d43,  speed,  heel) 

D6  =  interp(d46,  speed,  heel) 

dragw  =  1  /  72  rud  ^  2  yaw  ’^D4-l/72’^  rud  ^  2  yaw  ’^Dl  +  1/36’^ 
rud  ^  2  yaw  ’^D2-l/36’^  rud  ^  2  yaw  ’^D5-l/72’^  rud  ^  2  yaw 
D3  +  1  /  72  rud  ^  2  yaw  D6  +  1  /  18  rud  ^2’'Dl-l/9’'  rud  ^  2 

D2  +  1  /  18  rud  '^2'*'D3-l/8'*'  rud  yaw  D4  +  1  /  8  rud  yaw  D1 

-  1  /  6  rud  yaw  D2  +  1  /  6  rud  yaw  ’^D5  +  l/  24’^  rud  yaw  D3  - 

I  /  24  rud  yaw  D6  -  1  /  2  rud  D1  +  2  /  3  rud  D2  -  1  /  6  rud 

D3  -  1  /  4  yaw  D1  +  1  /  4  yaw  ^^04  +  01 

End  Function 

Function  interp(cond,  speed,  heel) 

Dim  1  As  Integer 
Dim  j  As  Integer 
Dim  heels  As  Variant 
Dim  heelrat  As  Double 
Dim  speeds  As  Variant 
Dim  speedrat  As  Double 
Dim  drags  As  Variant 
Dim  si  As  Double 
Dim  s2  As  Double 

Dim  11  As  String:  Dim  12  As  String:  Dim  13  As  String 

II  =  "A"  &  Trim (Str (cond  +  1)) 

12  =  "B"  &  Trim ( Str ( cond) ) 

13  =  "B"  &  Trim ( Str ( cond  +  1)) 

'if  speed  <  heel  ^  3  0.0002854  -  0.01427  heel  ^  2  +  0.2817  heel  +  1 

Then 

'  speed  =  heel  ^  3  0.0002854  -  0.01427  heel  ^  2  +  0.2817  heel  +  1 

'  Sheets ( "Calc" ) .Range ("Al") .Value  =  0 

'Else 

'  Sheets ( "Calc" ). Range ( "Al "). Value  =  1 
'End  If 

heels  =  Sheets ( "DATA" ). Range ( 11 ). Resize (2 1 ,  1) 

If  heel  >  heels (21,  1)  Then 
1  =  21 
heelrat  =  1 
Elself  heel  <=  0  Then 
1  =  2 

heelrat  =  0 

Else 


1 


1 
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Do  While  heel  >  heels  (1,  1) 
i  =  i  +  1 

Loop 

heelrat  =  (heel  -  heels (i  -  1,  1))  /  (heels (1,  1)  -  heels (i  -  1,  1) ) 

End  If 

speeds  =  Sheets ( "DATA" ). Range (12 ). Resize ( 1 ,  26) 

If  speed  >  speeds (1,  26)  Then 
j  =  26 

speedrat  =  1 
Elself  speed  <=  0  Then 

j  =  2 

speedrat  =  0 

Else 

j  =  1 

Do  While  speed  >  speeds (1,  j) 
j  =  j  +  1 

Loop 


speedrat  =  (speed  -  speeds (1,  j  -  1) )  /  (speeds (1,  j)  -  speeds (1,  j  - 

1)  ) 

End  If 

drags  =  Sheets ( "DATA" ) .Range (13) .Resize (i,  j) 

si  =  - (drags (i  -  1,  j  -  1)  -  drags (i,  j  -  1) )  ^  heelrat  +  drags (i  -  1,  j  -  1) 
s2  =  -  (drags  (i  -  1,  j)  -  drags  (1,  j))  heelrat  +  drags  (i  -  1,  j) 
interp  =  -(si  -  s2)  ^  speedrat  +  si 

End  Function 

Function  interptwo ( cond,  speed,  heel) 

Dim  i  As  Integer 
Dim  j  As  Integer 
Dim  heels  As  Variant 
Dim  heelrat  As  Double 
Dim  speeds  As  Variant 
Dim  speedrat  As  Double 
Dim  drags  As  Variant 
Dim  si  As  Double 
Dim  s2  As  Double 

Dim  11  As  String:  Dim  12  As  String:  Dim  13  As  String 

11  =  "AC"  &  Trim (Str (cond  +  1)) 

12  =  "AD"  Sc  Trim  (Str  (cond)  ) 

13  =  "AD"  Sc  Trim  ( Str  (cond  +  1)) 

heels  =  Sheets ( "DATA" ). Range ( 11 ). Resize ( 2 1 ,  1) 

If  heel  >  heels  (21,  1)  Then 
i  =  21 
heelrat  =  1 
Elself  heel  <=  0  Then 
i  =  2 

heelrat  =  0 

Else 

i  =  1 

Do  While  heel  >  heels  (i,  1) 
i  =  i  +  1 

Loop 


1,  1))  /  (heels (i,  1)  -  heels (i  -  1,  1) ) 


heelrat  =  (heel  -  heels  (1  - 
End  If 

speeds  =  Sheets ( "DATA" ). Range ( 12 ). Resize ( 1 ,  21) 
If  speed  >  speeds (1,  21)  Then 

j  =  21 

speedrat  =  1 
Elself  speed  <=  0  Then 
j  =  2 

speedrat  =  0 

Else 

j  =  1 

Do  While  speed  >  speeds (1,  j) 
j  =  j  +  1 

Loop 


speedrat  =  (speed  -  speeds (1,  j  -  1) )  /  (speeds (1,  j)  -  speeds (1, 

1)  ) 

End  If 

drags  =  Sheets ( "DATA" ) .Range (13) .Resize (1,  j) 

si  =  -  (drags  (1  -  1,  j  -  1)  -  drags  (1,  j  -  1)  )  heelrat  +  drags  (1  -  1,  j 

s2  =  -  (drags  (1  -  1,  j)  -  drags  (1,  j))  heelrat  +  drags  (1  -  1,  j) 
interptwo  =  -  (si  -  s2)  speedrat  +  si 
End  Function 

Function  Linlnt(cond,  speed,  heel,  starty,  startx,  lengths,  lengthh) 

Dim  1  As  Integer 
Dim  j  As  Integer 
Dim  heels  As  Variant 
Dim  heelrat  As  Double 
Dim  speeds  As  Variant 
Dim  speedrat  As  Double 
Dim  drags  As  Variant 
Dim  si  As  Double 
Dim  s2  As  Double 

Dim  11  As  String:  Dim  12  As  String:  Dim  13  As  String 

11  =  "starty"  &  Trim ( Str ( cond  +  1)) 

12  =  "startx"  &  Trim ( Str ( cond) ) 

13  =  "startx"  &  Trim (Str (cond  +  1)) 

heels  =  Sheets ( "DATA" ). Range (11 ). Resize (lengthh,  1) 

if  heel  >  heels (lengthh,  1)  Then 
1  =  lengthh 
heelrat  =  1 
Elself  heel  <=  0  Then 
1  =  2 

heelrat  =  0 

Else 

1  =  1 

Do  While  heel  >  heels  (1,  1) 

1  =  1  +  1 

Loop 

heelrat  =  (heel  -  heels (1  -  1,  1))  /  (heels (1,  1)  -  heels (1  -  1,  1)) 

End  If 

speeds  =  Sheets ( "DATA" ). Range (12 ). Resize ( 1 ,  lengths) 
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If  Speed  >  speeds (1,  lengths)  Then 
j  =  lengths 
speedrat  =  1 
Elself  speed  <=  0  Then 

j  =  2 

Speedrat  =  0 

Else 

j  =  1 

Do  While  speed  >  speeds (1,  j) 

j  =  j  +  1 

Loop 


speedrat  =  (speed  -  speeds (1,  j  -  1) )  /  (speeds (1,  j)  -  speeds (1 

1)  ) 

End  If 


drags  =  Sheets ( "DATA" ). Range ( 13 ). Resize ( i ,  j) 
si  =  -  (drags  (i  -  1,  j  -  1)  -  drags  (i,  j  -  1)  ) 
s2  =  - (drags (i  -  1,  j)  -  drags (i,  j))  ^  heelrat 
Linint  =  -  (si  -  s2)  speedrat  +  si 
End  Function 


heelrat  +  drags (i 
+  drags  (i  -  1 ,  j ) 


j 


j  - 


-  1) 


The  second  module  solved  for  the  sail  force  coefficients. 

Attribute  VB_Name  =  "Module2" 

Function  clmain (beta) 

If  beta  <  36  Then 

clmain  =  beta  -  4  ^  -0.0000033616  +  beta  "  3  0.00033299  +  beta  -  2  - 

0.01251  +  beta  0.2182 

Else 

clmain  =  0.000000186041  beta  ^  3  -  0.0000469267  beta  2  -  0.0080303 
beta  +  1.875128288 
End  If 

End  Function 

Function  cljib(beta) 

If  beta  <  7  Then 
cl jib  =  0 

Elself  beta  <  20  Then 

cljib  =  1.401482098  Log(beta)  -  2.720302423 
Elself  beta  <  50  Then 

cljib  =  -0.000231061  beta  "  2  +  0.014840909  beta  +  1.255606061 
Elself  beta  <  101  Then 

cljib  =  -0.00006  beta  "  2  -  0.0114  beta  +  2.14 
Else 

cljib  =  0 
End  If 

End  Function 
Function  cdmain (beta) 

If  beta  <  80  Then 

cdmain  =  0.0000406693  beta  2  -  0.000370797  beta  +  0.00593405 
Elself  beta  <  148  Then 

cdmain  =  -0.0000000219664  beta  "  4  +  0.00000858857  beta  "  3  -  0.001114913 
beta  ^  2  +  0.064256605  beta  -  1.272158531 


Else 

cdmain  =  0.0025  beta  +  0.75 
End  If 

End  Function 
Function  cdjib(beta) 

If  beta  <  28  Then 

cdjib  =  0.000225322  beta  +  0.004806867 
Else 

cdjib  =  -0.0000164176  beta  "  2  +  0.011980652  beta  - 
End  If 

End  Function 
Function  clspin  (beta) 

If  beta  <  27  Then 
clspin  =  0 

Elself  beta  <  59  Then 

clspin  =  0.0000269582  beta  3  -  0.005485895  beta 
beta  -  6.497055512 
Elself  beta  <  128  Then 

clspin  =  0.00000119265  beta  "  3  -  0.000518634  beta  " 

beta  +  0.248580029 

Else 

clspin  =  -0.00000177111  beta  "  3  +  0.000770699  beta 
beta  +  8.061621883 
End  If 

End  Function 

Function  cdspin (beta) 

If  beta  <  27  Then 
cdspin  =  0 

Elself  beta  <  59  Then 

cdspin  =  0.000357143  beta  ^  2  -  0.020714286  beta  +  0 
Elself  beta  <  99  Then 

cdspin  =  -0.000222527  beta  ^  2  +  0.048659341  beta  - 
Else 

cdspin  =  0.000000932728  beta  ^  3  -  0.000427881  beta 
beta  -  2.087528817 
End  If 

End  Function 


The  final  module  eontained  the  actual  VPP  solver  itself 

Attribute  VB_Name  =  "Module3" 

Public  rs(200,  17)  As  Double 
Public  bs(200,  17)  As  Double 
Public  cTWA  As  Range:  ' cTWA  =  "B4" 

Public  cSails  As  Range:  '  cSails  =  "M4" 

Public  cTWS  As  Range:  ' cTWS  =  "C4" 

Public  cHeel  As  Range:  ' cHeel  =  "L4" 

Sub  Runf latreef ( ) 

Dim  bs(700,  17)  As  Double 
Dim  flat  As  Double 
Dim  reef  As  Double 

Range ("J33") .Formula  =  0  'Set  rudder  matrix  to  zero 
Sheets ( "Calc" ). Select  'goto  calc  sheet 


0.313550059 

2  +  0.369098301 

2  +  0.051387093 

^  2  -  0.126128629 

.448928571 

1 . 626282967 

^  2  +  0.064506767 


Set  cellrange  =  Range (Cells ( 50 ,  1)  ,  Cells (750,  18)) 
counter  =  1 

For  reef  =  0.5  To  1.0001  Step  0.02 
Range ( "N5" ) -Value  =  reef 
For  flat  =  0.5  To  1.0001  Step  0.02 
Range ("N4") .Value  =  flat 

Iterate 

If  Range ("P17") .Value  =  True  Then 
Else 

Call  record (bs,  counter) 
counter  =  counter  +  1 
cellrange . Value  =  bs 
End  If 

Next  flat 
Next  reef 

End  Sub 


Sub  IterateWSO 
Sheets ( "Calc" ) .Select 

For  WA  =  30  To  171  Step  10 
Range ( "B4 " ) . Select  ' TWA 
ActiveCell.FormulaRlCl  =  WA 
If  WA  <  46  Then 

Range ("F40") .Formula  =  "=ABS(L26)"  'Max  this 
Elself  WA  <  165  Then 

Range ( "F4 0 "). Formula  =  "=ABS(M26)"  'Max  this 

Else 

Range ("F40") .Formula  =  "=ABS(L26)"  'Max  this 
End  If 

If  WA  <  89  Then 


Range ("M4") 

•Value  =  "MJ" 

'Sail 

Comb 

Range ( "M4 " ) 

•Value  =  "MS" 

'Sail 

Comb 

End  If 

Range  ("L4")  .Value  =  1.5 
Range ("C4") .Select  ' TWS 
'For  WS  =  6  To  20  Step  1 
WS  =  6 

ActiveCell.FormulaRlCl  =  WS 

ActiveCell . Of f set ( 0 ,  1) .Range ("Al") .Select 

Call  Iteratemid 

ActiveCell . Of f set (36,  -3) . Range ( "Al : Q1 " ) .Select 
ActiveCell . Of f set ( 0 ,  6) .Range ("Al") .Activate 

Selection . Copy 
Sheets ( "Results" ) .Select 

Selection . PasteSpecial  Paste : =xlValues ,  Operation : =xlNone, 
False,  Transpose : =False 
ActiveCell . Of f set ( 1 ,  0) .Range ("Al") .Select 
Sheets  ( "Calc" )  .Select 


SkipBlanks 


ActiveCell . Of f set (-36,  -4) .Range ("Al") .Select 
Application . CutCopyMode  =  False 

WS  =  12 

ActiveCell.FormulaRlCl  =  WS 

ActiveCell . Of f set ( 0 ,  1) .Range ("Al") .Select 

Call  Iteratemid 

ActiveCell . Of f set (36,  -3) . Range ( "Al : Q1 " ) .Select 
ActiveCell . Of f set ( 0 ,  6) .Range ("Al") .Activate 

Selection . Copy 
Sheets ( "Results" ) .Select 

Selection . PasteSpecial  Paste : =xlValues ,  Operation : =xlNone,  SkipBlanks 
False,  Transpose : =False 
ActiveCell . Of f set ( 1 ,  0) .Range ("Al") .Select 
Sheets  ( "Calc" )  .Select 

ActiveCell . Of f set (-36,  -4) .Range ("Al") .Select 
Application . CutCopyMode  =  False 

WS  =  16 

ActiveCell.FormulaRlCl  =  WS 

ActiveCell . Of f set ( 0 ,  1) .Range ("Al") .Select 

Call  Iteratemid 

ActiveCell . Of f set (36,  -3) . Range ( "Al : Q1 " ) .Select 
ActiveCell . Of f set ( 0 ,  6) .Range ("Al") .Activate 

Selection . Copy 
Sheets ( "Results" ) .Select 

Selection . PasteSpecial  Paste : =xlValues ,  Operation : =xlNone,  SkipBlanks 
False,  Transpose : =False 
ActiveCell . Of f set ( 1 ,  0) .Range ("Al") .Select 
Sheets ( "Calc" ) .Select 

ActiveCell . Of f set (-36,  -4) .Range ("Al") .Select 
Application . CutCopyMode  =  False 

WS  =  20 

ActiveCell.FormulaRlCl  =  WS 

ActiveCell . Of f set ( 0 ,  1) .Range ("Al") .Select 

Call  Iteratemid 

ActiveCell . Of f set (36,  -3) . Range ( "Al : Q1 " ) .Select 
ActiveCell . Of f set ( 0 ,  6) .Range ("Al") .Activate 

Selection . Copy 
Sheets ( "Results" ) .Select 

Selection . PasteSpecial  Paste : =xlValues ,  Operation : =xlNone,  SkipBlanks 
False,  Transpose : =False 
ActiveCell . Of f set ( 1 ,  0) .Range ("Al") .Select 
Sheets  ( "Calc" )  .Select 

ActiveCell . Of f set (-36,  -4) .Range ("Al") .Select 
Application . CutCopyMode  =  False 

'Next  WS 

Next  WA 
End  Sub 


Sub  Iteratemid 0 

Dim  bs(400,  17)  As  Double 

Dim  rs(400,  17)  As  Double 


Dim  flat  As  Double 
Dim  reef  As  Double 
Dim  counter  As  Integer 
Dim  cellrange  As  Range 
Dim  rudder  As  Double 

Sheets  ( "Calc" ). Select  'goto  calc  sheet 

Set  cellrange  =  Range (Cells ( 50 ,  1),  Cells (450,  18)) 

counter  =  1 

Range ("J33") .Formula  =  0  'Set  rudder  matrix  to  zero 
Range ("P14") .Value  =  1  'Turn  calc  off 

rudder  =  0:  Range ("K4") .Value  =  rudder  'Get  rudder  angle 

If  Range ("C4") .Value  <  13  Then  'Get  windspeed 
flat  =  1:  Range ("N4") .Value  =  flat  'flat 
reef  =  1:  Range ("N5") .Value  =  reef  'reef 
Elself  Range ("C4") .Value  <  17  Then  'Get  windspeed 
flat  =  0.94:  Range ("N4") .Value  =  flat  'flat 
reef  =  0.94:  Range ("N5") .Value  =  reef  'reef 

Else 

flat  =  0.9:  Range ("N4") .Value  =  flat  'set  initial  flat 
reef  =  0.9:  Range ("N5") .Value  =  reef  'get  initial  reef 
End  If 

Calculate 

Range ("P14") .Value  =  0  'Turn  calc  on 
movement  =  True 
rfd  =0.02 

If  flat  +  reef  =  2  Then 
movement  =  False 
End  If 

Iterate 

Call  record (bs,  counter) 
counter  =  counter  +  1 
stopit  =  0 

While  movement  =  True 
movement  =  False 

vvg  =  already (bs,  flat,  reef,  cellrange) 

If  vvg  =  0  Then 
Iterate 

Call  record (bs,  counter) 
counter  =  counter  +  1 
first  =  Range ("F40") .Value 

Else 

first  =  vvg 
If  counter  >  3  Then 
stopit  =  1 
End  If 
End  If 

flat  =  flat  +  rfd:  Range ("N4") .Value  =  flat 
If  flat  <  1  Then 

vvg  =  already (bs,  flat,  reef,  cellrange) 


If  vvg  =  0  Then 
Iterate 

Call  record (bs,  counter) 
counter  =  counter  +  1 
flatup  =  Range ("F40") .Value 

Else 

flatup  =  vvg 
End  If 

flat  =  flat  -  rfd:  Range ("N4") .Value  =  flat 
Else:  flatup  =  0 
End  If 

If  reef  <  1  Then 

reef  =  reef  +  rfd:  Range ("N5") .Value  =  reef 
vvg  =  already (bs,  flat,  reef,  cellrange) 

If  vvg  =  0  Then 
Iterate 

Call  record (bs,  counter) 
counter  =  counter  +  1 
reefup  =  Range ("F40") .Value 

Else 

reefup  =  vvg 
End  If 

reef  =  reef  -  rfd:  Range ("N5") .Value  =  reef 
Else:  reefup  =  0 
End  If 

flat  =  flat  -  rfd:  Range ("N4") .Value  =  flat 
vvg  =  already (bs,  flat,  reef,  cellrange) 

If  vvg  =  0  Then 
Iterate 

Call  record (bs,  counter) 
counter  =  counter  +  1 
flatdn  =  Range ("F40") .Value 

Else 

flatdn  =  vvg 
End  If 

flat  =  flat  +  rfd:  Range ("N4") .Value  =  flat 

reef  =  reef  -  rfd:  Range ("N5") .Value  =  reef 
vvg  =  already (bs,  flat,  reef,  cellrange) 

If  vvg  =  0  Then 
Iterate 

Call  record (bs,  counter) 
counter  =  counter  +  1 
reefdn  =  Range ( "F40" ) -Value 

Else 

reefdn  =  vvg 
End  If 

reef  =  reef  +  rfd:  Range ("N5") .Value  =  reef 
If  flatup  >  first  Then 

flat  =  flat  +  rfd:  Range ("N4") .Value  =  flat:  movement  =  True 
End  If 

If  reefup  >  first  Then 

reef  =  reef  +  rfd:  Range ("N5") .Value  =  reef:  movement  = 


True 
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End  If 

If  flatdn  >  first 
flat  =  flat  - 
End  If 

If  reefdn  >  first 
reef  =  reef  - 
End  If 


Then 

rfd:  Range ("N4") .Value 
Then 

rfd:  Range ("N5") .Value 


flat:  movement 


reef :  movement 


True 


True 


If  flat  >  1  Then 

flat  =  1:  Range ("N4") .Value  =  flat 
movement  =  False 
End  If 

If  reef  >  1  Then 

reef  =  1:  Range ("N5") .Value  =  reef 
movement  =  False 
End  If 

If  flat  <  0.5  Then  movement  =  False 
If  reef  <  0.5  Then  movement  =  False 


If  counter  >  40  Then  movement  =  False 

If  stopit  =  1  Then  movement  =  False 

'flat  =  flat  +  0.02:  Range ("N4") .Value  =  flat 

'reef  =  reef  +  0.02:  Range ("N5") .Value  =  reef 

'If  counter  >  3  Then 

'If  Range ("C4") .Value  >  16  Then 

'  If  Range ("F40") .Value  <  bs (counter  -  3,  6)  Then  'Max  this  variable 

'  reef  =  reef  -  0.02:  Range ("N5") .Value  =  reef 

'  flat  =  flat  -  0.02:  Range ("N4") .Value  =  flat 

'  End  If 

'End  If 

'End  If 

'If  flat  <  0.7  Then  reef  =  1.02 
'If  reef  <  0.7  Then  flat  =  1.02 
cellrange . Value  =  bs 

Wend 

If  flat  =1.02  Then 

flat  =  1:  Range ("N4") .Value  =  flat 
End  If 

If  reef  =  1.02  Then 

reef  =  1:  Range ("N5") .Value  =  reef 
End  If 
Iterate 

first  =  Range ("F40") .Value 
If  Range ("P17") .Value  =  False  Then 
Call  record (bs,  counter) 
counter  =  counter  +  1 
End  If 

cellrange . Value  =  bs 
where  =  MaxL(bs) 

flat  =  bs (where,  11) :  Range ("N4") .Value  =  flat 
reef  =  bs (where,  12) :  Range ("N5") .Value  =  reef 
Range ("F4") .Value  =  bs (where,  5)  'TBS 
Range ("L4") .Value  =  bs (where,  8)  'Heel 
Range ("J4") .Value  =  bs (where,  9)  'Yaw 
Range ( "K4 "). Value  =  bs (where,  10)  'Rudder 
Range ( "04 ") .Value  =  bs (where,  13)  'Rake 


Reset 
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dorudder : 


Dim  counterr  As  Integer 
Dim  ruddelta  As  Double 
Dim  error2  As  Double 

rudder  =  0:  Range ("K4") .Value  =  rudder 
ruddelta  =  2 
counterr  =  1 

Range  ("J33")  .Formula  =  "=Abs  ( 12 1 )  ^^2 "  'Set  rudder  matrix  to  not  zero 
While  rudder  <  6 
Iterate 

If  Range ("P15") .Value  =  True  Then  'Can  it  change  wrt  yaw? 

Range ("K4") .Value  =  rudder  -  ruddelta 
rudder  =  7 
GoTo  skip 
End  If 

Call  record (rs,  counterr) 

rudder  =  rudder  +  ruddelta:  Range ("K4") .Value  =  rudder 
counterr  =  counterr  +  1 
skip : 

' cellrange . Value  =  rs 

Wend 

ruddelta  =  1 
where  =  MaxL(rs) 


rudder 


flat  =  rs (where,  11) 

1  :  Range ( " 

N4") 

.Value  =  flat 

reef  =  rs (where,  12) 

1  :  Range ( " 

N5") 

.Value  =  reef 

Range ("F4") .Value  = 

rs (where. 

5) 

'  TBS 

Range ( "L4 ") .Value  = 

rs  (where. 

8) 

'  Heel 

Range ( "J4 "). Value  = 

rs  (where. 

9) 

'Yaw 

Range ("K4") .Value  = 

rs (where. 

10) 

' Rudder 

Range ("04") .Value  = 

rs (where. 

13) 

'  Rake 

rudder  =  Range ("K4") 

1  .Value 

Reset 

Iterate 

Call  record (rs,  counterr) 
counterr  =  counterr  +  1 
If  rudder  <  3  Then 
ruddelta  =  0.2 

Else 

ruddelta  =  -0.2 
End  If 


While  Abs (ruddelta)  >  0.05 
Iterate 

If  Range ("P15") .Value  =  True  Then  'Can  it  change  wrt  yaw? 
Range ("K4") .Value  =  rudder  -  ruddelta 
ruddelta  =  0 
GoTo  skip2 
End  If 

Call  record (rs,  counterr) 

error?  =  rs (counterr  -  1,  6)  -  rs (counterr  -  2,  6) 

If  error?  <  0  Then 
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ruddelta  =  -ruddelta  ^1/6 
End  If 

rudder  =  rudder  +  ruddelta:  Range ("K4") .Value  = 
counterr  =  counterr  +  1 


If  rudder  <  0  Then 

ruddelta  =  0:  rudder  =  0:  Range ("K4") .Value 
End  If 

If  rudder  >  6  Then 

ruddelta  =  0:  rudder  =  6:  Range ("K4") .Value 
End  If 


' cellrange . Value  =  rs 
skip2 : 

Wend 

cellrange . Value  =  rs 
where  =  MaxL(rs) 

flat  =  rs (where,  11) :  Range ("N4") 
reef  =  rs (where,  12):  Range ("N5") 
Range ("F4") .Value  =  rs (where,  5) 
Range ("L4") .Value  =  rs (where,  8) 
Range ("J4") .Value  =  rs (where,  9) 
Range ("K4") .Value  =  rs (where,  10) 
Range ("04") .Value  =  rs (where,  13) 
Reset 


.Value  = 
.Value  = 
'  TBS 
'Heel 
'  Yaw 
' Rudder 
'  Rake 


flat 

reef 


Iterate 


cellrange . Value  =  bs 
Range ("J33") .Formula  =  "0" 

Exit  Sub 

RemoveError : 

'flat  =  0.5:  Range ("N4") .Value  =  flat 

reef  =  reef  -  0.01:  Range ("N5") .Value  =  reef 

Range ("P14") .Value  =  1 

Calculate 

Range ("P14") .Value  =  0 
Calculate 


End  Sub 
Sub  Iterate  ( ) 

Dim  il  As  Long 
Dim  i2  As  Integer 

11  =  Application . Maxiterations 
Application . Maxiterations  =  10 

12  =  0 

Range ("P16") .Value  =  6 
Cagain : 

Calculate 
i2  =  i2  +  1 

If  Int(i2  /  40)  =12/40  Then 


rudder 


=  rudder 


=  rudder 


0.5  ^  2 


Range  ("P16")  .Value  =  (12) 

End  If 

If  Range ("P18") .Value  =  True  Then 
ResetCalc 

Range  ("PI  6")  .Value  =  Range  ( "PI  6"  )  .Value 
End  If 

If  Range ("L35") .Value  =  False  Then 
GoTo  Cagain 
End  If 

Application . Maxlterations  =  1 
End  Sub 

Sub  ResetCalc ( ) 

Dim  11  As  Long 

Dim  heel  As  Double 

11  =  Application . Maxlterations 

Application . Maxlterations  =  2 

heel  =  Range ("L4") .Value 

Range ("L4") .Value  =  Range ("L4") .Value  -  1 
Calculate 

Range ( "L4 "). Value  =  3 
Range ("P14") .Value  =  1 
Calculate 

Range ("P14") .Value  =  0 
Calculate 

Range ("L4") .Value  =  heel 
Range ("P14") .Value  =  1 
Calculate 

Range ("P14") .Value  =  0 
Calculate 

Application . Maxlterations  =  11 
End  Sub 

Function  MaxL (speeds)  As  Double 
bubba  =  0 

For  1  =  0  To  UBound ( speeds ,  1) 

If  speeds (1,  6)  >  bubba  Then 

bubba  =  speeds (1,  6) 

MaxL  =  1 
End  If 
Next  1 

'MsgBox  (speeds (MaxL,  6)) 

End  Function 

Sub  record (bs,  counter) 
bs  (counter  -  1,  0)  =  counter 


bs  (counter  - 

1, 

1) 

=  Range 

("A40") 

.  Value 

'TWA 

bs  (counter  - 

1, 

2) 

=  Range 

("B40") 

.  Value 

'TWS 

bs  (counter  - 

1, 

3) 

=  Range 

("C40") 

.  Value 

'AWA 

bs  (counter  - 

1, 

4) 

=  Range 

("D40") 

.  Value 

'AWS 

bs  (counter  - 

1, 

5) 

=  Range 

("E40") 

.  Value 

'TBS 

bs  (counter  - 

1, 

6) 

=  Range 

("F40") 

.  Value 

'VMG 

bs  (counter  - 

1, 

7) 

=  Range 

("G40") 

.  Value 

'ABS 

bs  (counter  - 

1, 

8) 

=  Range 

("H40") 

.  Value 

'  Heel 

bs  (counter  - 

1, 

9) 

=  Range 

("140") 

.  Value 

'Yaw 

145 


bs  (counter  - 

1, 

10) 

=  Range (" J40" 

) .Value 

' Rudder 

bs  (counter  - 

1, 

11) 

=  Range ("K40" 

) .Value 

'Flat 

bs  (counter  - 

1, 

12) 

=  Range ("L40" 

) .Value 

'  Reef 

bs  (counter  - 

1, 

13) 

=  Range ("M40" 

) .Value 

'  Rake 

bs  (counter  - 

1, 

14) 

=  Range ("N40" 

) .Value 

' Heelmom 

bs  (counter  - 

1, 

15) 

=  Range ("040" 

) .Value 

' Sidef ore 

bs  (counter  - 

1, 

16) 

=  Range ("P40" 

) .Value 

' Drive 

bs  (counter  - 

1, 

17) 

=  Range ("J34" 

) .Value 

'  Error 

End  Sub 

Function  already (bs,  flat,  reef,  cellrange) 

already  =  0 

cellrange . Value  =  bs 

For  i  =  0  To  UBound(bs,  1) 

If  Abs(bs(i,  11)  -  flat)  <  0.001  Then 

If  Abs (bs (i,  12)  -  reef)  <  0.001  Then 
already  =  bs  (i,  6) 

End  If 
End  If 
Next  i 

End  Function 


Velocity  Prediction  Program  Spreadsheet 

The  next  few  pages  eontain  screenshots  of  the  contents  of  the  actual  cells  in  the  tank  data 

spreadsheet. 
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From  the  “calc”  sheet  which  was  the  main  engine  were  these  cells: 


1  ^  •  Security...  ^  S  , 


.23333 

J  z  z  z  z 


r 


lo  "o  ~ 

S _ M  Q 


Ismail 


N  3  I 

1  Uff 


*-  <N  ^  F-  o  W  ' 

»-  0>  ro  00  O  PI 
iv>  rr)  in  Oi  -g  ^ 


31  I-  ^  fN  01  ^ 

JO  01  00  00  00^^  til 
N  <N  (N  (N  rvj  1+’*^ 

g  III;  I 


]  Paste  Special...  g  1] 
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File  Edit  View  Insert  Format  lools  Data  Window  Help  Solver...  Goal  Seek...  Aaobat  -  IbI  X| 
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from  the  sail  calculations: 


■  Eji  VPP25J(LS 


